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ABSTRACT 

An assessment has been made of the magnitude and significance of 
the pollution loadings from urban runoff in the province of Ontario. The 
study was conducted under the provisions of the Canada-Ontario and 
Canada-United States Agreement on Great Lakes Water Quality. 

The selected local authorities were interviewed to obtain 
detailed local data. In addition, various maps and demographic informa- 
tion were supplied by the Ministry of the Environment. Methodology used 
by the American Public Works Association and the University of Florida 
for a study of similar problems in the United States provided a basis for 
data manipulation and preparation of cost estimates, modified where 
possible to reflect conditions relevant to Ontario. 

The cost estimate performed indicated that to obtain 25 percent 
control of BOD, an annual cost of $10,861,000 would be incurred. These 
costs are exclusive of the storm flow conveyance system. For BOD control 
at 75 percent, the annual cost would be $95,471,000. 



RESUT1E 



Nous avons evalufi 1 'importance de la pollution causae par les 
effluents urbains en Ontario, dans le cadre des accords Canada-Ontario et 
Canada - Etats-Unis sur la quality des eaux des Grands lacs. 

Nous avons interrog^ des responsables locaux afin d'obtenir des 
donn(5es locales d^tailK^es. Le mlnistSre de 1 'Rnvironnement nous a, en 
outre, fourni diverses cartes et informations d^moqraphiques. Les 
methodes utilisSes par I' Arnerican Public Works Association et 
1 'university de la Floride pour 1 'etude de probl^nes sernblables aux 
Etats-Unis nous ont serv1 de module pour traiter les donn^es et estimer 
les coOts. Nous les avons adaptcs, lorsque possible, aux conditions 
propres 3 1 'Ontario. 

[Jos estimations indiquent qu'une reduction de 25 p. 100 de la 
DBO coQterait 10 861 000 de dollars annuellement, en excluant les coOts 
du systeme d'^vacuation des eaux pluviales. Une rf^ductlon de la DBO de 
75 p. 100 coQterait 95 471 000 dollars par an. 
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SUMMARY 

In 1975 the American Public Works Association, with the University 
of Florida, conducted a study for the Urban Drainage Subcommittee of the 
Canada-Ontario Agreement on Great Lakes Water Quality. 

The study was conducted by utilizing methods and procedures 
developed for a similar study of the United States. Information was 
supplied directly from the Ontario Ministry of the Environment, Environ- 
ment Canada, and field interviews conducted in ten representative cities. 

The characterization of urban storm water runoff and combined 
sewer overflows was taken from the United States study, modified slightly 
for conditions in Ontario. Information concerning pollution from 
snowmelt was obtained from Canadian studies. 

Storm water modelling using the complete U.S. EPA model STORM 
for four small cities, and similar methods for an additional 52 cities, 
were used to determine storm flows and potential pollution loads. 
Relationships were developed between population density and pollution 
control costs for separate and combined sewered areas. Based upon an 
assumed cost and availability of land by population density, optimization 
of storage versus treatment was considered. 

Control of pollution has been generally limited to evaluation of 
BOD, although other parameters are being controlled at the same time at 
various degrees of efficiency. The approach used could be adapted for 
other quality parameters. 

An important aspect of water quality planning is the tradeoff 
involved in the decision-making process when alternatives such as 
providing advanced dry-weather treatment and control of storm water or 
combined sewer overflows are under consideration. The study indicates 
that a significant portion of the wet-weather pollution should be 
controlled prior to initiation of advanced wastewater treatment. 

The study found that local officials in Ontario were ^^ery 
concerned with flood control aspects of storm flow and were less 
concerned with quality aspects. Information required for detailed 
individual model studies was not available, nor were plans being made to 
gather information on key parameters. 
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Thus, local officials should be made aware of the itnportance of 
pollution from urban storm water runoff and combined sewer overflows, and 
the significance of such polluted discharges compared to their present 
discharge of treated wastewater treatsnent plant effluents. 

The storage requirements necessary to economically treat storm 
flow may also serve to alleviate local flooding problems. Likewise, the 
gathering of data for storr.i water niodelling would oe of direct benefit in 
the planning and improvement of the drainage system, a readily perceived 
benefit. 

Although the cost calculations may not oe accurate to a high 
degree for an individual city, it is believed that they dre reliable for 
preliminary assessments when considering the total urban drainage to the 
Great Lakes from the province of Ontario. 
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CONCLUSIONS 

1. City officials perceived many problems with their storm water 
control systems. Local flooding due to hydraulic overloading, 
infiltration/inflow, deposition of solids, and untreated bypasses 
were the most common concerns. 

2. Interest in modelling storm water runoff exists, particularly with 
regard to quantity predictions. Most authorities appeared to not 
have the necessary resources to conduct studies. 

3. To facilitate modelling, extensive efforts were made to characterize 
urban land use patterns. The distribution of developed land use in 
urban areas was taken as: 



Residential 


52.5 


Commercial 


10.3 


Industrial 


14.0 


Other 


23,2 



i4.* The extent of combined sewer systems was determined for 49 of the 
56 cities, based on available data. 

5. The loading factors used to calculate pollution loadings were 
based upon a study of available applicable data. However, the 
overall receiving water quality impact of various sources has not 
been evaluated in the urban setting. Among the areas where research 
is minimal are: snowmelt, wear products from street surfaces, urban 
sediments and erosion products, tree and leaf litter, and accumula- 
tion from non-street impervious areas. 

6. Primary treatment devices using physical processes such as screen- 
ing, settling, and flotation have been developed and tested for 
application to combined sewer overflows. Their application to urban 
storm flows should be equally effective. BOD removal efficiency of 
40 percent appears reasonable. 

7. For all urban areas of over 10,000 population draining into the 
Great Lakes, the annual cost of providing 25i BOD removal to urban 
storm flow and combined sewer overflows using secondary treatment 
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and storage is estimated at $10,861,000, and for 50 percent control, 
$31,744,000. 

8. Secondary treatment devices which use biological and physical- 
chemical processes are suitable for treating both combined sewer 
overflows and urban storm flow. Contact stabilization is feasible 
only if the dry-weather flow (DWF) facility is of an activated 
sludge type. BOD removal efficiency of 85 percent appears 
reasonable. 

9. The annual cost of providing 75% BOD removal for urban storm flow 
and combined sewer overflows using secondary treatment and storage 
is estimated to be $95,471,000. 

NOTE: The above costs are in terms of 1975 dollars. Costs include 
land, engineering, and sludge disposal but do not include the cost 
of a sewer system to transport flows to the point of treatment. The 
primary purposes of the cost estimate were to: (1) develop an order 
of magnitude for the costs which might be anticipatted as control 
measures are required; and (2) develop a procedure for estimating 
costs as more refined and specific data become available. 

10. Most DWF facilities offer opportunity for treating wet-weather 
flow (WWF). The optimum mix of facilities must be determined on a 
case-by-case basis. 

11. The optimal cost for WWF facilities is a function of the size of 
the facility, the unit cost of which decreases with size, and the 
cost of conveying the flow to the facility, a cost which Increases 
with the size of the facility. Such determinations must be made for 
specific sites. For the purpose of this study a representative 
treatment cost was determined. 

12. The optimal sizing of treatment and storage facilities will vary 
for each area and for the level of control required. 

13. For BOD removals of less than or equal to about 10 percent, the 
optimal treatment strategy is to use primary treatment devices for a 
portion of the flow. Secondary treatment devices will be required 
for higher levels of control. 
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14. The relationship between level of control and number of overflow 
events can be predicted by methods which were developed. 

15. To achieve 50 percent BOD removal, intensive control of combined 
sewer overflows results in a lower marginal cost than control in 
storm sewer areas. 

16. Operation and maintenance costs will be affected by the number of 
hours that the facility is operated. The amount of storage which is 
provided allows a smaller capacity treatment facility which will 
operate for a longer period of time to achieve the same relative 
treatment. For the purpose of this study, annual operation and 
maintenance costs have been assumed to be 20 percent of the total 
costs of the treatment facility. Actual costs will vary by the type 
and size of the Individual treatment units. 

17. The assumption of a first flush, a high concentration of pollu- 
tants in the first portion of the runoff, has a significant effect 
on cost assumptions. Control costs are about one-third less If a 
first flush is assumed. 

18. The cost-benefit relationship between tertiary treatment of DWF 
and provision of WWF facilities needs to be Investigated prior to 
instituting either control measure. When analyzed under the assump- 
tions used in this report, it was found that about 16 percent of the 
storm flow should be controlled before Initiating tertiary treatment 
if additional removal of organics is the purpose of tertiary 
treatment. 

19. Rooftop and parking lot storage, surface and underground tanks, 
and storage in treatment units are effective flow attenuation 
control alternatives. The cost of providing such facilities varies, 
based upon whether or not they can be designed Into new or reconstruc- 
tion projects, and on the population density of the area in which 
they are to be located. 

20. The unit costs of pollution control are lowest In the unsewered 
areas because of relatively low storage costs. 
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RECOMMENDATIONS 

1. Local authorities should be made aware of the benefits of and need 
for urban storm flow planning for quality as well as quantity 
considerations. Local officials expressed concern with flooding 
problems but exhibited relatively minor interest in the polluting 
aspects of storm flow. 

2. The Ministry of the Environment should take the lead in initiating 
storm water modelling efforts. Local officials showed an interest 
in modelling but generally did not have the resources to undertake 
such studies. The results of model studies should be of major 
assistance in planning for storm flow quantity. 

3. Inventories of local demographic characteristics are needed, and 
compilation should be made to assist further modelling efforts. 
Such inventories should include the area served by each type of 
sewer. 

4. Reports dealing with runoff and quality predictions should be 
carefully structured. Quality parameters must be precisely defined 
along with the averaging method used. Structured demographic data 
such as population density and land use should be reported for each 
area where studies are conducted. 

5. The cost assessment methodology should be extended to cover 
additional storage and treatment combinations - cost and performance 
data for storage and treatment units, and the impact of different 
storage reservoir operating policies. The relationship between the 
number of overflow events and percent runoff controlled would be 
useful . 

5. Careful investigation of assumptions used in this study concerning 
percent imperviousness is needed. Lower unit control costs may be 
anticipated as it is believed that conservative functions were used. 

7. The tradeoff between storm flow control and tertiary treatment 

should be evaluated on the basis of control of pollutants other than 
suspended solids and BOD. 
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8. In conducting receiving water studies to detennine the effect of 
urban storm flow pollution, the water quality indicators that will 
be used for planning purposes should be identified before the start 
of data collection. The data collection system must be designed to 
obtain representative samples before, during, and after storm 
events. 

9. Modelling efforts require additional work. Such work should 
include: 

a. The response of receiving waters to urban runoff and 
dry-weather flow inputs should be characterized when storage of 
waste streams is considered in combination with treatment. 

b. Simplified techniques to approximate the complex mechanisms of 
pollutant transport in lakes and bays should be developed. 

10. Results of present studies of runoff discharges from completely 
developed urban drainage basins should be carefully evaluated. Time 
related responses of these systems to flow and concentration for a 
variety of rainfall and runoff events will assist future modelling 
efforts in the analysis of discrete samples of runoff to provide 
quality information, and provide indications of runoff characteriza- 
tion over time. 

11. The use of multipurpose units should be given consideration due to 
the potential cost savings. The following figure includes the 
estimated cost by percent pollutant (BOD) content for single purpose 
and multipurpose systems. 
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1^ INTRODUCTION 

1 «1 Purpose 

The purpose of the study was to ascertain the magnitude and 
significance of pollution loadings from urban storm water runoff. The 
three principal objectives of the study were to: 

1) prepare planning estimates of the quantity and quality of 
urban storm water runoff contributions to the Great Lakes 
Ontario watershed. 

2) develop cost estimates for implementing control and abate- 
ment practices with presently available technology. 

3) transfer the methodology employed in the analysis to 
technical personnel designated by the Urban Drainage 
Subcommittee. 

\^M Scope 

The American Public Works Association and the University of 
Florida have previously prepared an assessment of the pollutant effects 
of urban storm water runoff and combined sewer overflows [1,2,3]. Data 
so developed and information from the project officers, as well as from 
field visits to ten Ontario municipalities, were used for this study. 
Computer models were programmed and used to predict pollution loadings 
and to make overall cost assessments for various storage-treatment combina- 
tions at the downstream end of an entire urban catchment. Model outputs 
allow prediction of given levels of control which may be stated as: 

1) percent of runoff captured, 

2) percent of BOD or other pollutant removed. 

3) number of overflows per year, and/or 

4) quantity of overflows per year. 

1 .3 Background 

Canada-Ontario Agreement (COA) Research Report No. 26, "Review 
of Canadian Design Practice and Comparison of Urban Hydrologic Models," 
October 1975, by J.F. MacLaren Ltd. [64], describes lucidly the concepts 
involved in storm water management. 



"storm water management, considering drainage as a subsystem 
of the total urban system with environmental aspects and possible 
benefits, is a relatively new concept. 

"The traditional storm drainage design philosophy was to collect 
the runoff and carry it away as fast as possible out to the boundaries of 
the considered watershed. This was done by connecting all impervious 
areas such as roofs and driveways to a network of gutters and conduits 
with considerably higher velocity and density than in the natural 
drainage system. Storm water was also considered clean and there was no 
concern with regard to pollution from separate storm sewers. The design 
of the storm sewer system was carried on independently from the studies 
for flood control from rare events. Negative consequences of this 
philosophy, such as drastic increase of the peak flows at the outlet of 
the urbanized watershed, increased incidence of local flooding, depletion 
of groundwater, considerable increase in the cost of new storm sewerage 
systems and relief sewers, and the environmental damage, are now evident 
and many attempts for an Innovative approach are underway. 

"The key to the Implementation of new management methods, 
however, is the use of improved hydrologic tools. The design of storage, 
for example, which is the simplest method for reduction of flow peaks, is 
possible only through the synthesis of hydrographs. Storage in an urban 
system is not necessarily concentrated In a reservoir but may be 
distributed over different elements of the watershed such as parking 
lots, roofs, elements of the sewer network, etc. Other methods of peak 
reduction are the retardation of flow by reduction of velocity or 
increase of infiltrated volumes. The traditional design method for 
drainage systems, the Rational Method, Is aimed at providing only design 
peak flows and cannot be used for the study of management techniques. 
Even the use of the Rational Method for the derivation of design peak 
flows has been subject to numerous criticisms. 

"Therefore, an increasing number of more sophisticated models 
dealing with urban storm water runoff, some of which include quality 
considerations, are being developed." 

1 ,4 Runoff in General 

Little doubt now exists that urban storm water runoff 
represents a significant source of water pollution, and affects the 
quality of our streams, estuaries, lakes, and oceans. 

Considerable research has taken place to better understand the 
contamination of runoff In both urban and nonurban environments. Of 
particular Interest, however. Is urban surface runoff and its contribu- 
tions to the deterioration of receiving water quality. The polluting 
effects of runoff may be classified In terms of direct and indirect 
pollution. Direct pollution Includes discharges In runoff from separate 
storm sewer collection systems, or contributions within uncontained 
runoff entering the receiving water at locations other than clearly 
defined points of discharge. 



Indirect pollution involves point discharges or overflows from 
the planned or unplanned addition of storm water to other wastewater 
flows. This may include combined sewer overflows and overflows resulting 
from uncontrolled runoff inflow into sanitary sewerage systems and, in 
some cases, excessive infiltration. 

Traditionally, direct pollution from runoff has been disregarded. 
Surface runoff was generally characterized as a resource to be quantita- 
tively controlled. Drainage and flood control objectives were paramount 
and runoff pollution was considered nonexistent or at most a low priority 
problem. Although early investigative efforts in Europe [4] and the 
United States [5] began to suggest the importance of surface runoff 
pollution, serious consideration of its effects is a fairly recent 
phenomenon. It was not until a 1964 report by the U.S. Public Health 
Service that the issue of runoff quality began to assume national 
importance [6]. In the ensuing period, a number of research efforts have 
sought to characterize runoff pollution, its impact, and its control and 
abatement. 

One approach has been to empirically characterize discharges in 
various drainage basins across the country. This has often involved the 
study of drainage flows from urban or urbanizing drainage basins. On 
occasion, relationships between discharge and effluent quality data have 
been related to physical basin characteristics and given rainfall events. 
Inconsistency exists within this body of information, however, because of 
variation in research objectives being addressed, the pollutants being 
evaluated, the sampling techniques used and the measurements made. 

These studies have identified a number of contaminating consti- 
tuents. Some may demand considerable amounts of oxygen. Other 
constituents, such as pathogenic organisms, may produce the risk of 
infectious disease. Some contaminants are nutrients capable of promoting 
the growth of algae and aquatic plants within a receiving water. Others 
may be toxic to plants and animals. Still others adversely affect 
natural stream purification processes or give rise to sediment and solids 
depositions. From the subsequent water user's point of view, some 
pollutants cause a water supply to become hard, corrosive, or may render 
otherwise potable water unacceptable from the standpoint of colour. 



turbidity, odour, or appearance. In the same vein, many of these 
constituents overtax existing water treatment facilities or make their 
operation uneconomical. 

In an urban space, pollutants may be deposited for subsequent 
pickup by surface runoff or they may be directly introduced into drainage 
flows. The products of combustion and other suspended materials in the 
air (particulates and other emissions) may be scavenged from the air by 
falling rain. Deposits of airborne materials on pervious and impervious 
surfaces may be removed to contaminate runoff flows. Street paving and 
surfacing materials; debris from open areas, including erosion products, 
organic plant and animal wastes, and a variety of chemicals such as 
fertilizers, soil conditioners, and pesticides; transportation related 
material including deposits of fuel, lubricants, hydraulic fluids, 
coolants, tire, clutch and brake wear products, exhaust emission 
particulates, rust and dirt; street litter, household and commercial 
wastes; and finally snow and ice control, antiskid and corrosion 
inhibiting materials all may contribute to the contamination of runoff 
and its subsequent effects on receiving water quality. 

Sediment is perhaps the largest single source of water pollution, 
Current estimates suggest that four billion tons [7] of sediment makes 
its way to the rivers of North America annually. Sediments are soils or 
other surficial materials that are products of erosion and may be 
transported or deposited by the action of wind, water, snow, ice, or 
gravity [8]. 

Erosion and sedimentation are natural and continually occurring 
geological processes. Normally, soils are protected by vegetation and 
vegetative residue. In areas where moisture is too limited or fertility 
too low to sustain close-growing vegetation, the land is subject to 
periodic erosion from intense rains. Man's actions, including construc- 
tion and many types of urban activity, often remove vegetation in 
localized areas, which tends to increase the rate of erosion. Removal of 
the protective cover allows the forces of wind and water to act more 
directly and forcefully on the exposed soil particles. 

Nonpoint source pollutants are organic and inorganic materials 
entering storm water from unspecific or unlocalized sources in sufficient 



quantity to constitute a pollution problem. In a rural environment, they 
include sediment, plant nutrients, pesticides, and animal wastes from 
cropland, rangeland, pastures, and farm feedlots. Sediment is the major 
pollutant in tenns of volume [9], and may be a carrier of some pesticides 
and plant nutrients. In an urban enviroment, similar pollutants may be 
experienced from impervious areas as well as those materials that are 
unique to urban activities, such as transportation related pollutant 
sources, air pollution, and so on. 

A body of knowledge is now being developed through the study of 
some of the pollutant source characteristics previously described. 
Although this area of study was developed primarily for nonurban environ- 
ments and nonpoint discharges, some generalizations are now being applied 
in urban cases to estimate pollution effects. The use of the Universal 
Soil Loss Equation for the estimation of sediment contributions [10] is a 
good example of a nonurban technology used in appropriate urban applica- 
tions. 

In urbanized areas, the polluting potential of street litter 
accumulations has been studied to assess the magnitudes of the pollutants 
that are available to surface runoff. Considering the developed urban 
street as a temporary sink for the accumulation of pollutants that dre 
representative waste products of a complex urban environment, methods for 
estimating the quantity of runoff pollution have been devised under the 
assumptions that the urban street is a logical extension of the urban 
drainage system and that the runoff and polluting contributions from 
pervious areas will be negligible for most runoff events. This approach 
to the mechanism of urban runoff pollution may be construed as a special 
case of the study of contaminant source characteristics. 

All of these methods represent some of the various mechanisms 
that have been used for the assessment of the direct contributions of 
urban runoff to pollution. The priorities associated with the evalua- 
tion, abatement, and control of indirect pollution have generally been 
much higher. Indirect contributions are due to the admixture of runoff 
with other wastewater flows. Interest in uncontrolled discharges of 
combined sewer overflows has generally taken the form of sampling 
programs, and contributions to pollution have been determined through 



discharge measurements on a case-by-case basis. Similarly, the control 
and abatement of combined sewer overflows has been developed on a site- 
specific basis. 

These pollution phenomena are a fact of life. They have been 
accepted as a matter of fact in the past and their impact on receiving 
waters has had an adverse effect on the water resources of the United 
States and Canada, and on most if not all of the developed nations. 

The eventual need to control, minimize, or eliminate the pollu- 
ting effects of urban and nonurban runoff wastewaters is incontestable. 
Yet, the costs of achieving this goal are so high that investments must 
be proven necessary and essential in terms of the benefits to be derived. 
At the same time that idealistic efforts must be made to overcome runoff 
pollution, further funds will be required to upgrade dry-weather waste- 
water treatment standards to reduce or eliminate the impacts of 
inadequately treated or untreated sanitary flows on receiving waters. 
In addition, it may be necessary to control, at least in some measure, 
the polluting impact of nonurban surface drainage wastes runoff on the 
same streams, lakes, and other water sources which are affected by 
wastewaters from urbanized areas. 

Of overriding importance in establishing urban wastewater 
policies for the future is knowledge of the polluting potentials of stomi 
water discharges, combined wastewater overflows, and dry-weather flow 
spin constituents, and their comparative or relative impacts on receiv- 
ing waters, and also clarification of the "natural" cleanliness or 
pollution condition of natural streams and other receiving waters before 
any urban runoff waters reach them. 

In order that the present study be wholly practical , commentaries 
on pertinent information gained from field interviews in Burlington, 
Guelph, Kingston, Kitchener, Milton, St. Catharines, Saul t Ste. Marie*, 
Thunder Bay, Toronto, and Windsor, are contained in Section 2. 

Land use information relating to population, land area and 
location, and population density and land use distribution are examined 
in Section 3. 

Utilizing such demographic information, together with pollution 
loading information set out in Section 4, a province-wide cost assessment 



was prepared. This assessment can be improved or expanded as further 
data become available, as detailed in Section 5 of the report. 



2. THE PROBLEM IN ONTARIO 

Problems associated with receiving water pollution from urban 
storm water runoff have slowly become recognizable because sanitary waste 
discharges from local authorities have been corrected by wastewater 
treatment facilities. Combined sewer overflows and storm water discharges 
have been identified as having the potential to become a major source of 
pollution in receiving waters. 

In Ontario it was found that few local authorities were concerned 
with or had identified problems associated with storm water runoff at this 
time. Rather, attention was primarily focused upon problems associated 
with flood control aspects, such as flooded basements and overloaded 
sewers. 

2.1 Known Sewer System Problems as Perceived in the Province 

On-site interviews were conducted at the ten representative 
cities shown in Figure 1 by a trained representative of the APWA. 
Officials of the Ministry of Environment were present for many of the 
sessions. The interview outline is contained in Appendix III. A broad 
range of questions were asked to determine basic infonnation and to allow 
the Ministry to gain insight into local perceptions of the storm water 
pollution problem. The survey attempted to identify the extent of local 
problems with four general conditions: 

a. hydraulic overloading of sewers, 

b. existence of system bypasses, 

c. solids deposition in sewers, and 

d. infiltration/inflow problems. 

Hydraulic overloading is caused by either inadequately sized 
sewers or excessive flows. The result of overloading is generally 
bypassing of treatment facilities. System bypasses allow excessive flows, 
or combined sewer overflows, to escape from the system without treatment. 
Solids deposition in sewers results from the use of large or combined 
sewers, flat grades or poor interior surfaces. Solids are deposited 
during dry weather and are then flushed out of the system during stonn 
events. Infiltration and inflow represent a multitude of problems which 
allow either surface or groundwater to enter the system and contribute to 
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FIGURE 1 ONTARIO-GREAT LAKES BASIN TEST CITIES. 



bypasses or excessive flows at treabnent facilities. All of these 
problems may be interrelated; all affect receiving water quality. 

Table 1 summarizes the responses received as either major or 
minor, or as not a problem. Hydraulic overloading due to size of pipe and 
existence of infiltration/inflow conditions appears to be the major 
problem. 

TABLE 1. SEWER SYSTEM PROBLEM AREAS IDENTIFIED 



City 



Overloaded 
Sewers 



Bypasses 



Solids 




Deposited 


Infiltration/ 


In Sewer 


Inflow 


* 


X 


% 


t 





- 


-r 


X 


«• 


« 



Burlington 

Guelph 

Kingston 

Kitchener 

Milton 

Sault Ste. Marie 

St. Catharines 

Thunder Bay 

Toronto 

Windsor 















X 






Key: X identified general problem 
identified minor problem 
- not identified as problem 

Source: APWA Survey, 1975 



2.2 



General Comments From Field Interviews 



During the course of the field interviews many relevant 
comments were made by local officials which are of importance in 
evaluating the overall problems in the province. This section will 
highlight the major comments. 
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2.2.1 Windsor 

The City of Windsor encounters problems in the area known as 
Riverside. Riverside is served by a separate system with many interconnec- 
tions. In some areas the sanitary sewer is laid directly under the storm 
sewer with a common manhole and a plate at the invert of the storm sewer. 
This leads to many problems: leakage due to improper plate replacement, 
missing plates, etc. Tapping into the sanitary sewer system is also very 
difficult. 

There are a number of sewer sections that are substantially 
deficient in hydraulic capacity. 

Stonn runoff causes flooding in basements and ponding on the 
streets. Stom drainage facilities are inadequate in some parts of the 
city. Most of the sewer problems occur in early and late winter. 

Windsor's past experience in receiving grant-in-aid funds has 
been good and the city would like to receive additional grants. City 
officials prefer the Central Mortgage and Housing Corporation grants 
because of the simple grant requirements as compared to provincial grants 
from Ontario Ministry of Transportation and Communication for stonn water 
drainage projects. 

2.2.2 Burlington 

Burlington has established an allocation procedure for funding 
corrective work. However funding is considered inadequate. 

Many sewers in the older areas are overloaded, primarily due to 
underdesign. Burlington also has serious infiltration problems. 
Infiltration has resulted from open joints and poor bedding. Many 
subdivisions were built between 1955/1960 with only sanitary sewers; 
open ditches are used for storm drainage. 

Basement and backyard flooding in the older areas is attributed 
to storm sewers backing up because of high rates of inflow and infiltra- 
tion. This situation generally occurs during the January thaw. As far as 
abatement of runoff is concerned. Burlington uses neither storage nor 
detention facilities. There have been no studies of storm water quality. 
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2.2.3 St. Catharines 

The financial arrangements for stonn and sanitary wastewater 
works in St. Catharines appear good. However, programs are cut to meet 
the monies available. Grants-in-aid may act as a restraint as well as an 
encouragement. 

Hydraulic overloading shows up in areas of recurrent flooding. 
This is believed to be due to the undersizing of the combined sewers. 

The three-year cleaning frequency of the sewers appears to be 
effective. There has been no concurrent improvement of receiving waters 
although a long-term program is under study which includes designs for 
watercourses. 

There is no formal receiving water sampling program. There has 
been some effort to verify volumetric and flow gauge data obtained from 
their sampling activities. Local officials appeared Interested in runoff 
planning and hoped to learn how to apply storm water management models for 
input to their proposed study. 

Heavy infiltration/Inflow is believed to be the cause of 
widespread flooding all over the city. The watercourses were designed 
with inadequate cross-sectional areas. The problem is under study with 
redesign planned in the near future. ""~~-^ 

The practice of building extensions on the separated sewers and 
then letting these sewer extensions drain into the combined sewers has 
only aggravated the pollution problem. ____ 

The bypassing of sewage during rainstorms is due to the hydraulic 
inadequacies of the system. This will be corrected by a planned program 
of improvement. 

2.2.4 Kingston 

The financial arrangements for stonn and sanitary wastewater 
works are fair to good. The city relies generally upon the Ontario 
Ministry of Transportation and Communication for grants-in-aid. They have 
found that other grant programs involve considerable overhead and 
therefore they are avoided unless the advantages are substantial. 

Infiltration and inflow are high, amounting to five to seven 
times dry-weather flow (DWF) during peak wet-weather flow periods (WWF). 
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The causes are probably bad joints on sewers and laterals, broken sewers, 
and bad backfill practices. 

Basement flooding is only a minor problem because of four things: 
1) the construction of additional storm sewers, 2) the elimination of the 
cross connections between storm and sanitary sewers, 3) disconnection of 
roof leaders, and 4) Improved sequencing and operation of pumps. 

Kingston expressed some interest in runoff planning. They had no 
quantity /quality management investigation planned or running at the 
time. They appeared to be very interested In runoff planning if Ministry- 
level assistance were available. 

2.2.5 Thunder Bay 

Monies received from the federal Central Mortgage and Housing 
Corporation were not substantial, about 16 percent of system expenditures. 
The city had great praise for the Central Mortgage and Housing Corporation 
and the simplicity of the grants application. Thunder Bay has experienced 
little difficulty in securing grants. 

There are no real legal difficulties in requiring the removal of 
the roof leader connections. However, no attempt has yet been made to 
enforce removal . 

There are many Inadequately sewered sections of the city. Storin 
water collection systems have not always been installed. Basement 
flooding is prevalent in the Fort William area. 

The city was constructing sewer system relief bypasses with the 
hope that such construction would eliminate the need for additional storm 
sewers. 

2.2.6 Sault Ste. Marie 

Allocation of money for storm and sanitary wastewater works in 
the past appears to have been at a reasonable level. 

Sault Ste. Marie has no fixed program for maintenance of their 
sewer lines; however, the city department cleans the entire system once 
every three years. 

Planning activities for quantity /quality had not been conducted. 
However, Sault Ste. Marie was very responsive to the idea of using runoff 
planning in the future. 
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Judging from the records of the wastewater treatment facility, 
infiltration and inflow must be assumed to be high. A general infiltra- 
tion problem exists all over the city. Flooding of streets is not a 
serious problem, 

2.2.7 Toronto 

Grant experience has been satisfactory. The relief afforded by 
the road storm sewers is not sufficient. However, the city had not 
anticipated major relief from such storm sewers. 

Storm water runoff causes underpasses and basements to flood, due 
primarily to hydraulic overloading. 

Both the Metropolitan and City of Toronto governments are interest- 
ed in runoff planning and would provide input to a study sponsored by the 
Ministry. Toronto has computer capacity and good programmers, but increased 
sophistication is required. Ministry assistance is, of course, a policy 
matter, but it is believed that a cooperative program would be well received 

Toronto presently has no runoff abatement plans. Storage or 
detention is not applicable at the moment and, therefore, likely locations 
are not known. 

2.2.8 Kitchener 

The City of Kitchener has enjoyed a good relationship with the 
federal government. Kitchener's record for applying for grants and 
obtaining them has been good. 

A study by Proctor and Redfern, Consulting Engineers, on bottle- 
necks, both present and potential, in the sanitary collection systan, 
should identify and allow correction of areas with any major solid 
deposition problems. 

Consultants have recommended that an open storage reservoir be 
built in a natural catchment by damming a branch of Schneider Creek. This 
would in effect reduce hydraulic flows in that particular branch of 
Schneider Creek, and hence, act as a storage/detention basin. The basis 
for such a recommendation is the Grand River Conservation Authority report. 

Roof leader connections are troublesome. There is a history of 
allowing footing drains to be connected to the sanitary sewers because the 
storm drains are so shallow. Insofar as the disconnection of roof leaders 
is concerned, enforcement has not been attempted. 
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There was great interest expressed in learning how to apply 
runoff planning tools to the Kitchener area, 

2.2.9 Guelph 

Guelph's past experience with the receipt of grants-in-aid has 
been basically good, although paperwork requirements are considered to be 
too extensive. 

In the older areas ponding is a problem because sewers are 
inadequate. 

The water usage per capita is increasing and relief trunk sewers 
may be required. It is believed that a nonlinear rate structure would 
help discourage excess use. 

Local deficiencies are recognized In certain areas; sanitary 
storm sewers, infiltration and inflow, to name a few. Local officials are 
interested only in available runoff planning tools if the cost-benefit is 
apparent. Local officials would be able to produce local inputs by 
providing manpower for monitoring. 

2.2.10 Milton 

The Town of Milton has experienced severe problems with its 
sewerage system. The main reason is because of the incomplete system in 
the old part of town. This area is served by ditches and culverts. The 
town hopes that this problem will be eliminated in the near future. 

Basement flooding is usually caused by sanitary sewers backing up 
through the foundation drains. There are no combined sewer overflow 
regulators and there are no detention facilities. 

Local authorities are interested in runoff planning, but through 
the town's consulting engineers. As far as computer capacity is concern- 
ed, there is no available source with a sufficient memory, and they would 
happily accept the assistance of the Ministry of the Environment. 

2.3 Climate 

Climatic information was obtained which covered altitude above 
sea level, mean annual temperature, mean maximum daily temperatures, mean 
minimum daily temperatures, daily range of temperatures, extreme low and 
high temperatures, mean annual precipitation in inches, mean annual water 
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surplus in inches and mean daily temperature. Three government publica- 
tions were consulted in this connection. They were: "The Climate of 
Southern Ontario". Climatological Studies No. 5, Environment Canada. 1974 
[61]; "The Climate of Northern Ontario," Climatological Studies No. 6, 
Environment Canada, 1968 [65]; and "The Climate of the Great Lakes Basin," 
Climatological Studies No. 20, Environment Canada, 1972 [60]. 

A detailed listing of all the factors is included in the summary 
of responses from the ten cities furnished separately and is not included 
in this report. Selected factors are shown in Table 2. It was ascertained 
from a study of the climatic data that in some cases although the rainfall 
frequency for a particular year was typical, the snowfall was not. 
Because of this, one year of data from four cities, Burlington (1973), St. 
Catharines (1973), Kingston (1965), and Sault Ste. Marie (1969). was 
chosen as being typical for the computer analysis. This selection of four 
cities to represent the whole region was deemed adequate because of the 
fact that monthly precipitation below 44° latitude is very stable, 
while above 44° the trend is towards a peak in the summer months. 
Snowfall was also an important component of Ontario's precipitation total, 
ranging from a low of 4.0 in. (10.2 cm) water equivalent near Windsor to a 
high of 14.0 in. (36 cm) north of Sault Ste. Marie. 

2.3.1 Windsor 

From the book "The Climate of Southern Ontario" [61], the 
following statistics were taken for the City of Windsor. The mean annual 
temperature is 48°F (9°C). The extreme low temperature for the year 
has been recorded as -21^V (-330c) while the extreme high temperature 
has been recorded at lOe^F (41°C). The mean annual precipitation is 
30 inches {76 cm); the mean annual snowfall is 40 inches (102 cm). The 
annual water surplus in inches is somewhat lower than for other parts of 
Ontario, being in the order of 10 inches (25 cm). 

2.3.2 Burlington 

In Burlington, the climatic area is described as the Lake 
Ontario shore. The city is at altitude of about 300 feet (91 m) above 
sea level. The mean annual temperature is 450F (70c). The 
extreme low temperature is -300F (-340C) and the extreme high is 
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TABLE 2. POPULATION, ALTITUDE, AREA, AND CLIMATE OF INTERVIEWED CITIES 



City 


Pop 
1973 


Area/ 
1,000 
acres 
(hectares) 


Altitude 

feet 
(metres) 


Annual 

Precipit' n 

inches 

(cm) 


Snowf al 1 

inches 

(cm) 


Days of 

Precipit' n 


Mean 
Temp 
Op 
(OC) 


Burl ington 


91,554 


13.8 
(5.6) 


300 
(91.2) 


34 
(86) 


65 
(165) 




45 

(7.2) 


Guelph 


63.009 


11.3 
(4.6) 


1,150 
(350) 


33 
(84) 






44 
(6.67) 


Kingston 


59,289 


5.7 
(2.3) 


245-355 
(75-108) 


35 
(89) 








Kitchener 


121,441 


22.1 
(8.9) 


1,100 

(335) 


34 
(86) 






44 
(6.67) 


Milton 


15,667 




650 
(198) 


30-38 
(76-97) 


48 
(122) 




45 
(7.2) 


Sault Ste. 
















Marie 


77.501 


12.7 
(5.1) 


700 
(213) 


37 
(94) 


112 
(284) 




40 

(4.44) 


St. Catharines 


112,299 


16.3 

(6.6) 


300-800 
(91.3-243) 


31 
(79) 


45 

(114) 




40 

(8.9) 


Thunder Bay 


105,954 


141 
(57.1) 


645 

(196) 


29 

(74) 


73 

(185) 


141 


36 

m 


Toronto 


676,363 


24.0 
(9.7) 




26.7 
(68) 


59.2 
(150) 


134 


47 
(8.3) 


Windsor 


199,250 


26.3 

(10.6) 




30 
(76) 


40 
(102) 


137 


48 

(8.9) 



104OF (40^0). The mean annual precipitation is 34 inches (86 cm) 

and the mean annual snowfall is 65 inches (165 cm). The mean annual water 

surplus is about 13 inches (33 cm), 

2.3,3 St. Catharines 

St. Catharines straddles the Niagara escarpment. This deserves 
special attention because of its influence on the climate of the contiguous 
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region. The escarpment extends from Queenston, where it is a 300 foot 
(91 m) bluff along the south shore of Lake Ontario through Hamilton northward 
to Collingwood on Georgian Bay where it is almost 8,000 feet (2,440 m) 
high. The annual precipitation is 31 inches (79 cm). The mean annual 
temperature is 480f (9°C). The extreme low temperature is -160F (-27°C) 
and the extreme high temperature is 1040F (400C). The mean annual snowfall is 
45 inches (114 cm). The altitude above sea level is approximately 300 
feet (91 m). 

2.3.4 Kingston 

Located in Southern Ontario, Kingston's climate is moderated by 
the Great Lakes. Hence, Kingston's summers are cooler and winters milder 
than those in eastern Ontario or in the United States west of the Great 
Lakes. In July the mean temperature range is from 640f to 700F 
(17.70c to 21.1°C). The corresponding mean temperature in January 
Is 180F to 240F (-80C to -4^0). The effect of the lakes, particularly 
in the summer when sunny days bring cool lake breezes, tends to lower the 
maximum temperatures. Precipitation is evenly distributed throughout the 
year with annual precipitation ranging from 30 to 40 Inches (76 to 102 
cm). The altitude above sea level is about 300 feet (91 m) . 

2.3.5 Thunder Bay 

Thunder Bay has an altitude of 600-1,400 feet (183-427 m) above 
sea level. The mean annual temperature is 360f {20c). The mean 
annual minimum temperature is -40OF (-40OC). The mean annual snow- 
fall is 73 inches {185 cm) and the water surplus is 5 Inches (13 cm). 

2.3.6 Sault Ste. Marie 

The information obtained about Sault Ste. Marie is taken from 
"The Climate of Northern Ontario". The statistics are as follows: 
altitude, about 700 feet (213 m) ; mean annual temperature is 40OF 
(4.40c); mean annual minimum temperature is -30^V (-34°C); the 
mean annual precipitation in inches is 33 (84 cm); and the mean annual 
water surplus in Inches is 13 (33 cm). The mean annual snowfall is 112 
inches {284 cm). 
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2.3.7 Toronto 

Toronto, like the rest of the cities around the Great Lakes 
area, has its weather significantly altered by the Great Lakes. The 
usual general circulation over Ontario gives Toronto a temperature range 
from lOOp {-120c) (January) to 700f {2^^Z) (July) and a mean annual 
rainfall of 26.7 inches (68 cm). The winters are cold and dry with an 
annual snowfall of 59.2 inches (150 cm). The mean annual water surplus is 
about 11 inches (28 cm) . 

2.3.8 Kitchener 

Kitchener is located in the climatic area known as the Huron 
Slope. The mean annual temperature Is 44°F (6.60c). The extreme 
low temperature recorded Is -430f {-420c). The extreme high temperature 
Is 102OF (390c). The mean annual precipitation varies between 32 and 39 
inches (81 and 99 cm). The mean annual water surplus is 6 inches (15 cm). 
The altitude above sea level is about 1,100 feet (336 m) . 

2.3.9 Guelph 

The City of Guelph is also situated in the climatic region 
known as the Huron Slope and has a mean altitude above sea level of 
approximately 1,150 feet (351 m) . The mean annual temperature Is 440f 
(6.60c). The extreme high temperature has been 104Of (40Oc), 

2.3.10 Milton 

The mean annual temperature is 450f (70c). The daily range in 
temperature Is 180f (IOOc) in January and 220f (11°C) in July. The 
extreme low temperature has been recorded as -390f (-390c) the extreme high 
is 1050F (40,5Oc). The mean annual precipitation Is given as between 30 and 
38 inches (76 and 97 cm) of rainfall. The mean annual water surplus Is 
12 inches (30 cm). The mean annual snowfall is 48 inches (122 cm). The 
altitude above sea level Is approximately 650 feet (198 m) . 

2*4 Wastewater Treatment and Collection Systems 

Wastewater treatment facilities were also examined. As 
nutrient removal (phosphates) is required by Ontario regulations, most 
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plants are of the secondary type (activated sludge plus physical-chemical) 
For most of the cities examined the responsibility for wastewater treat- 
ment lies either with the provincial Ministry of Environment with costs 
borne by the municipality served, or with the appropriate regional 
government. Exceptions are Guelph, Windsor, and Kingston where all 
collection and treatment are city responsibilities. In Toronto, waste- 
water treatment is the responsibility of the Municipality of Metropolitan 
Toronto. Table 3 shows the prevailing situation. 

A summary by cities of treatment plant type, hydraulic capacities 
and percentage removal of pollutants is included in Table 4. 

Receiving waters are as follows: 



2.4.1 



Windsor 

Burlington 

Kingston 

Thunder Bay 

Sault Ste. Marie 

Toronto 

Kitchener 

Guelph 
Milton 
St. Catharines 

Windsor 



Detroit River to Lake Erie 

Lake Ontario 

Cataraqui River to Lake Ontario 

Kaministiqua River to Lake Superior 

St. Mary River to Lake Huron 

Don and Number Rivers to Lake Ontario 

Schneider Creek and 

Speed River to Lake Erie 

Eramosa River to Lake Erie 

16 Mile Creek to Lake Ontario 

Creeks to Lake Ontario 



Windsor is basically served by two sewage treatment plants. 
One is served by combined sewers; the other is served by separate sewers. 
The one located in West Windsor has a capacity of 21 Imgd (1104 L/sec) 
while the other, located in Little River, Ontario, has a capacity of 4,5 
Imgd (236 L/sec). The Little River is an activated sludge plant with 
phosphorus removal. Dry-weather flow goes directly into the interceptor 
sewer. 

The West Windsor plant has summary treatment along with 
phosphorus removal. The Little River plant is equipped with pumping 
stations, grit chambers, sedimentation basin, aeration tanks, secondary 
sedimentation, vacuum filters, and centrifugation. 
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TABLE 3. RESPONSIBILITY FOR SEWER SYSTEM 



City 


Sanitary 
Sewers 


Storm 
Sewers 


Combi ned 

Sewers 


Wastewater 

Treatment 


Windsor 


City 


City 


City 




City 


Burlington 


Region (Halton) 


City 


N/A 




Region 


St. Catharines 


Region (Niagara) & City 
Depending on size 


City 


Shared 




Region 


Kingston 


City 


City 


City 




City 


Thunder Bay 


Ministry of Environment 


City 


City 




Ministry of the 
Environment 


Sault Ste. Marie 


City 


City 


City 




Ministry of the 
Environment 


Toronto 


Shared according to size 
with Mun. of Metro Toronto 


Shared 


Shared with 
Metro Toronto 


Metro 
Toronto 


Kitchener 


City 


City 


City 




Region (Waterloo) 


Guelph 


City 


City 


N/A 




City 


Mi 1 ton 


Region (Halton) 


Town 


N/A 




Region 



TABLE 4. SUMMARY OF WASTEWATER TREATMENT FACILITIES 



City 



Type Number Total Mean Daily 
of of Flow Treated 
Treatment Plants { Imgd) (m^/s 



Efficiency of Removal 
BOD5 SS Nutrients 
[%) {%) {%) 



Burl ington 


Activated 
Sludge 


3 


1.57 
10.6 
1.55 


0.083 
0.557 
0.08 


94.5 
85.6 
55.4 


90.8 
78,3 
40.8 


55.4 
82.4 
27.1 


Windsor 


Primary 

Activated 

Sludge 


2 


22.1 
5.0 


1.16 

0.263 


48.0 
74.0 


62.0 
80.0 


70.0 
41.0 


St. Catharines 


Activated 

Sludge 

Primary 


2 


15.4 


0.81 


86.2 
57.0 


90.5 
47.6 


58.9 
12.8 


Kingston 


Primary 


1 


11.0 


0.58 


50.0 


75.0 


- 


Thunder Bay 


Primary 


2 


6.01 
6.20 


0.32 
0.33 


23.0 
22.0 


49.0 
52.0 


22.0 
28.5 


Sault Ste. Marie 


Primary 


1 


9.90 


0.52 


43.0 


66.0 


38.0 


Toronto 
Metro 


Conventional 
Secondary 


3 


69.0 

176.0 

8.1 


3.63 
9.26 
0.43 


9b. 
82.0 
87.0 


95.0 
91.0 
93.0 


87. U 
93.0 


Kitchener 


Primary 


1 


14.3 


0.75 


89. 0(?) 


90. 0(?) 


- 


Guelph 


Conventional 
Secondary 


3 


Not Operating 
3.0 
3.0 


0.16 
0.16 


77.0 
81.5 


90.5 
93.0 


- 


Milton 


Conventional 
Secondary 


1 


1.0 


0.053 


89.7 


94.7 


- 
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2.4.2 Burl jngton 

The sanitary wastewater works in Burlington were recently 
turned over to the regional government. The operation and maintenance of 
the storm sewer system is the responsibility of the city. The operation 
and maintenance of the sanitary sewer system is the responsibility of the 
region. There are no combined sewers in Burlington, but there Are 
separate storm sewers and open channels. At the time there were no 
detention or retention facilities. 

The mean daily flows in Burlington's three wastewater treatment 
plants are 1) Drury Lane, 1.S7 Imgd (82 L/sec); 2) Skyway, 1U.6 Imgd (552 
L/sec); and 3) Elizabeth Gardens, 1.55 Imgd (81 L/sec). The treatment 
unit processes that they employ are: 1) Drury Lane-conventional activated 
sludge; 2) Skyway - extended aeration, aerobic digestion, and sludge 
thickening; and 3) Elizabeth Gardens-conventional activated sludge. 

2.4.3 St. Catharines 

In the City of St. Catharines there are two types of sanitary 
sewers: regional trunk sanitary sewers and local sanitary sewers. The 
wastewater treatment plant is solely the responsibility of the region 
with charges being made on a user basis. The operation and maintenance of 
the storm sewers is the responsibility of the city. All the pumping 
stations are maintained by the region at the cost of the city. 

Regulator facilities are set to overflow when the wet-weather 
flow exceeds 2.5 times the dry-weather flow. There is no detention. 
However, the city has developed a project proposal for evaluation of the 
performance of a combined wastewater retention facility. In areas where 
there are combined sewers, restricting runoff from roofs and holding on 
open spaces is practiced. 

The Port Weller treatment plant has a mean daily flow of 15.4 Imgd 
(810 L/sec). The Port Weller treatment process is activated sludge, while 
Port Dalhousie has a primary treatment unit. 

2.4.4 Kin gston 

For Kingston, the wastewater treatment works' capacity is about 
19 Imgd (990 L/sec). Treatment unit processes consist of pumping 
stations, grit chambers, sedimentation tanks, and digesters. The effluent 
is chlorinated and discharged into Lake Ontario. 
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Effluent strength for BOD is between 5-100 mg/L. Nutrients are 
about 5-7 mgA and are mostly phosphates. The efficiency of the treatment 
process is approximately 5-50 percent for BOD and 75 percent for suspended 
solids. 

Limited size detention tanks were under consideration. Chlorina- 
tion is proposed if the detention tanks are used, but there were no 
current plans for the wastewater treatment plant. The average flow for 
the plant is 11.0 Imgd (578 L/sec). Primary treatment with separate 
sludge digestion is used. 

2.4.5 Thunder Bay 

In Thunder Bay the sanitary wastewater works are the responsi- 
bility of the Ministry of Environment of Ontario with the sanitary waste- 
water sewers being designed by the city and operated by the Ministry. The 
design, construction and operation of the sewage treatment plant is also 
the responsibility of the Ministry. Thunder Bay has two wastewater 
treatment plants. South Plant. 6.01 Imgd (314 L/sec) and North Plant, 6.20 
Imgd (324 L/sec). Expansion of the South Treatment Plant was underway and 
the North Plant was to be abandoned shortly thereafter. The South Treat- 
ment Plant contains screening, pumping, grit removal, primary sedimenta- 
tion, chlorination, and two stage separate sludge digesters. No detention 
facilities were planned at the time. 

2.4.6 Sault Ste. Marie 

Sault Ste. Marie is served by a separate storm sewer and a 
sanitary sewer running side by side and entering a common manhole. 

The wastewater treatment plant and the pumping stations are 
maintained by the Ministry of the Environment. Insomuch as dry-weather 
treatment is concerned, effluent flows can be adjusted from the lift 
stations as necessary, and through the bypass capability; these are used 
for emergency purposes less than once a year. 

There are no detention facilities in use with the exception of a 
conservation scheme known as the Fort Creek Conservation Dam. There are 
no control activities for wet-weather flow. 

The mean daily flow is 9.90 Imgd (516A sec). The plant uses 
comminution, grit removal, primary sedimentation, and chlorination as 
treatment processes. 



2.4.7 Toronto 

In Toronto, there are two classes of sewers: those that are 
owned and operated by the city, and the Metro sewers. The Metro sewers 
are all trunk sewers that serve 1,000 acres (405 ha) or more. The city 
sewers come under the jurisdiction of the Department of Public Works of 
Toronto. All the wastewater treatment works are controlled by the Metro 
Corporation and the same division prevails for operations maintenance and 
capital construction. 

Ninety-five percent of the city is served by combined sewers. There 
is a separate sewer district covering Swansea (discharges into the Humber 
plant) and Forest Hill Village (discharges into the Toronto system). 

There were no wet-weather control activities at the time. The 
mean daily flows in Toronto's three wastewater treatment plants are: 1) 
Humber, 69.0 Imgd (3,600 L/sec); 2) Main, 176.0 Imgd (9,260 L/sec); and 
3) North, 8,1 Imgd (426 L/sec); all employ conventional secondary 
treatment. 

2.4.8 Kitchener 

In Kitchener all sewers are the responsibility of the city. 
The treatment works are the responsibility of the regional government. 
The sewage treatment plants, all of which are called regional, are, for 
all intents and purposes, provincial sewage treatment plants. There are 
very few combined sewers, as there are separate sanitary and storm sewer 
systems. The treatment plant has been under continuous expansion. 

Two pumping stations have bypasses to the rivers. The sewage 
treatment plant bypass operates only when a breakdown takes place. The 
effluent is chlorinated. Storage methods are not used. None of the other 
techniques, such as tanks or detention facilities, are employed. The mean 
daily flow in the Kitchener wastewater treatment plant is 14.3 Imgd (746 
L/sec), utilizing secondary treatment. 

2.4.9 Guelph 

Guelph has the responsibility for the design and construction 
of the storm water facilities. The city is also responsible for the 
wastewater treatment works. This 1s unusual in Ontario. The plant is a 
secondary treatment plant with phosphorus removal. In addition, the 



effluent is chlorinated. Operation and maintenance of the plant, as well 
as the storm and sanitary sewer systems, is the responsibility of the 
city. 

Wastewater treatment works handle a mean daily flow of about 8,5 
Imgd (433 L/sec) and peak is somewhat higher. The flow is gravity fed to 
the plant, which is a conventional activated sludge plant. Actually, 
there are three plants in Guelph, but all are at one geographical 
location. One has a capacity of 2 Imgd (104 L/sec) and is the oldest of 
the plants; two newer plants have a capacity of 4 Imgd {208 L/sec). The 2 
Imgd (104 L/sec) capacity plant is presently held as a standby. Detention 
and retention basins were expected to be constructed during the 1975-1976 
season. 

2.4.10 Milton 

Milton has a conventional activated sludge plant, with tertiary 
treatment using lime. This treatment process is a regional responsibility. 
The rated capacity is 1.58 Imgd (82.4 L/sec); the average daily flow is 
1.0 Imgd (53 L/sec). The raw sewage has a five-day BOD of 151 mgA, 
suspended solids of 195 mgA, nitrogen of 51 mg/L, and a phosphorus total 
of 7.5 mg/L. The final effluent has a five-day BOD of 6.5 mgA, suspended 
solids of 24 mgA, nitrogen of 7.5 mgA, and phosphorus of 2.9 mgA. The 
percentages of removal are therefore 95.6, 37.7, 85,3, and 61.3, respec- 
tively, for five-day BOD, suspended solids, nitrogen, and phosphorus. The 
Ministry of the Environment is committed to allowing a connected popula- 
tion of 18,000 persons to contribute to the existing plant. Milton has 
established a subdivision priority list to increase the population to 
18,000. A further expansion to 2.84 Imgd (148 L/sec) will allow for an 
additional 10,000 persons. Treatment facilities at the plant are >^ery 
erratic due to pulsing caused by effluent pumps. 

There are no wet-weather controls and detention facilities in 
Milton at this time. The town, however, still physically operates and 
maintains its own sewers. 

2.5 Development Characteristics 

Key city development characteristics are listed in Table 5. As 
noted, most cities have considerable industry. 



TABLE 5. DEVELOPMENT CHARACTERISTICS 



City 



Major Industries 



Growth Expected 



Windsor 
Burlington 
St. Catharines 

Kingston 

Thunder Bay 

Sault Ste. Marie 
Toronto 

Kitchener 



Guelph 
Milton 



Automotive (Primary) 

Distilling 

Secondary types 
Steel 

Secondary Automotive 

Agricultural 

Electronic 

Sparse with aluminum 
and paint industries 
in adjacent township 

Grain shipping & storage 
Pulp & Paper 
Steel & Iron 

One prominent industry 
Steel 

Almost every conceivable 
type of industry {primary 
and secondary) and 
commercial endeavor 

Furniture manufacturing 

Automotive (secondary) 

Rubber 

Meat Packing 

Shoe manufacturing 

Diversified secondary 
Steel (secondary) 



Moderate 



Uniform 



Moderate 



Small 



Moderate 

Uniform 

Uncertain due 
to political 
development disputes 



Brisk 

Steady - at 
about 4% p. a. 

Brisk 



2.5.1 



Windsor 



Situated across the river from the city of Detroit, Windsor is 
Canada's counterpart of Detroit as an automotive center. Ford, Chrysler, 
General Motors, and American Motors plants are in the city. Windsor is 
also a distilling center. There are light support industries and the 
University of Windsor within the City boundaries. 
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2.5.2 Burlington 

In Burlington, secondary industries predominate, along with 
service industries of the "clean" type. Slater Steel was expected to come 
into Burlington in late 1975. Population growth is expected to stay 
uniform. 

2.5.3 St. Catharines 

St. Catharines is strongly supported by the automotive 
industry. There are some electronic plants, a winery, and canning and 
agricultural industry sidelines. The construction industry is prominent. 
The western part of the city is devoted to agricultural use, fruit raising 
and vineyards, some sparse commercial and industrial uses, and some 
suburban-type residential development. There are no utilities other than 
electric and communication. The built up area contains mixed commercial, 
industrial, and residential developments. 

2.5.4 Kingston 

Kingston has little local industry. There is a sparse amount 
of light industry including a dairy. Kingston is a university center and 
includes Queen's University, The Royal Military College of Canada, and 
St. Lawrence Community College. It is also a tourist and historical 
center, serving the Thousand Island region in the St. Lawrence River. 

2.5.5 Thunder Bay 

Thunder Bay is one of Canada's main grain shipping ports. 
There are four large pulp and paper mills, a major steel and iron 
fabricating manufacturing plant, and some chemical manufacturers. 
Industrial support industries include saw mills, a malting plant, and an 
industrial grain product plant. 

The city is a tourist center and serves as gateway (eastern) to 
northwestern Ontario. It is well supplied with service industries and is 
the location of Lakehead University and Confederation College. It has 
shipyards. 

2.5.6 Sault Ste. Marie 

Sault Ste. Marie is virtually a one industry town with a steel 
company dominating. Smelted ore is shipped by rail and water to Sault 



Ste. Marie from Wawa to the North, and from the Mesabi and Atikokan 
ranges to the west by rail and water. There is a local iron foundry and 
numerous service industries. 

2.5.7 Toronto 

In Toronto, one of the most difficult things is predicting changes 
in land use because of the development dispute that is going on in the 
various segments of the city. Toronto has a yery wide diversity of land 
use. Residential density varies from a few persons per acre in such communi- 
ties as Rosedale to as high as 600 persons per acre (1,480 persons per 
hectare) in St. James Town. Toronto is also an important transportation 
centre and a port city. The City of Toronto is distinct from Metro Toronto, 
which is responsible for a variety of area-wide services. Major economic 
activities may be described as business, commerce, manufacturing, transporta- 
tion, which literally cover the entire field of production of modern goods 
and services, Toronto being the prime supplier for all of Canada. 

2.5.8 Kitchener 

Kitchener is a heavily industrialized city. Major economic 
activities include furniture manufacturing, automotive manufacturing, tire 
manufacturing, meat packing, shoe manufacturing, and other industries. 
Commercial activities are average, as are service activities. 

2.5.9 Guelph 

Guelph contains a wide range of commercial, industrial, and 
cultural activities. A sound economic base is apparent. Present 
indications suggest that the population should reach 130,000 by the year 
2000. The population of Guelph is increasing at an annual rate of about 
four percent. The northwest area of the city has not been developed while 
the remainder of the city has been developed. 

2.5.10 Milton 

Development pressure is wery heavy for the Town of Milton. The 
urban area of Milton still contains about 90 percent agricultural open 
space. Major activities in the area include a screw fastening company and 
a steel products company. There are many other smaller industries 
scattered throughout the town. 
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2.6 Summary 

Four main problems, namely, hydraulic overloading, system 
bypassing, solids deposition, and infiltration/inflow in the ten Canadian 
cities were confirmed. Undoubtedly there are many other problems such as 
flooding and runoff, but the problems commented upon in this report 
appeared to be most prevalent in almost all the cities of Ontario. 

Broadly characterized, the problems which appeared to exist are 
common to many cities in the United States. Inasmuch as there has not 
previously existed a major reason for attention to many of the problems of 
sewer overloading, bypassing, and such, little attention apparently has 
been given to the sewer system, the buried "conducts of civilization". 

The general opinion was that grants can act either as a restraint 
or as an encouragement. When a grant policy is in effect, programs are 
designed according to the availability of funds. A summary of Federal and 
Provincial Grant regulations is contained in Appendix IV. Local 
authorities have generally had favorable experience with Federal (CMHC) 
and Provincial (MOT&C) grants. 

Interest in runoff planning and management was at a high level 
and all communities seemed willing to cooperate with the technology 
transfer proposed and would provide input wherever possible. Most, 
however, were unable to provide computer capacities for the recently 
available storm water runoff models such as STORM [11]. 

Most communities did not have storm water quality /quantity 
management programs in operation although some work was being carried out 
by consultants. 



3. DESCRIPTION OF THE URBAN AREAS 

This section presents a summary and analysis of data on the 
following characteristics of urban areas within the portion of Ontario 
located in the Great Lakes Basin: 

1) population, land area, and location; and 

2) population density and land use distribution. 

3.1 Definition 

Urban areas in this study have been defined as: 

1) an incorporated city or urban core of 10,000 or more 
inhabitants; or 

2) an agglomeration of census tracts with population densities 
of one or more persons per acre. 

Using information found in the "1971 Census of Canada", Statistics 
Canada [12], a total of 56 urban areas were defined. The resulting 
population and land area for each urban area are shown in Table 6. In 
order to characterize all of the urban areas in the study area, ten test 
cities were chosen before the study was initiated: Burlington, Guelph, 
Kingston, Kitchener-Waterloo, Milton, St. Catharines, Sault Ste. Marie, 
Thunder Bay, Toronto, and Windsor. Milton was not used due to its small 
population. Toronto was characterized by a section of that municipality 
know as West Toronto, except for the determination of population density 
groups, in which the entire city was used. The location of these test 
cities is shown in Figure 1. 

3.2 PoP"lAt:i'_on DensT ty and Land Use Di stri bution 

The overall population density for an urban area may be 
obtained using the data in Table 6. In general, population densities 
have decreased during the past generation, reflecting the availability of 
improved transportation systems, the desire for individual home owner- 
ship, etc. No detailed data on urban land use for all of the urban 
areas could be found. For the nine urban areas (see Appendix II for 
maps), the area occupied by each of the following five types of land uses 
was determined: residential, institutional, industrial, commercial, and 
open space. Land use maps were unavailable. Aerial photographs were 
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employed if land use maps were unavailable. These photos presented a 
problem, in that they were taken at altitudes (1.6 miles or 2.6 km, 
approximately) that made differentiating among land uses difficult. 

The percentage of the land use in each of these five categories 
is shown in Table 7. Note that, with the exception of West Toronto which 
has a very high average population density (29.7 persons per acre), a 
large percent of the urban area falls in the open space category. This 
large amount of open space resulted from the definition of urban areas, 
which included population densities as low as one person per acre. Thus, 
significant acreages of land which were undeveloped and would not be 
served by sewerage systems were included. 

TABLE 7. LAND DISTRIBUTION IN NINE ONTARIO CITIES 

Population 
Density Percent of Urbanized Area in Indicated Land Use 

Urbanized Area persons/acre Residential Commercial Industrial Institutional Open Total 



3.0 


3.0 


5.0 


55.0 


100 


3.0 


8.0 


6.0 


54.0 


100 


5.0 


7.0 


12.0 


44.0 


100 


4.0 


10.0 


2.0 


57.0 


100 


4.0 


6.0 


6.0 


43.0 


100 


4.0 


19.0 


4.0 


33.0 


100 


13.0 


12.0 


6.0 


35.0 


100 


19.0 


14.0 


4.0 


5.0 


100 


6.0 


10.0 


3.0 


43.0 


100 



Burlington 5.84 34.0 

Guelph 5.01 29.0 

Kingston 10.39 32.0 

Kitchener-Waterloo 6.70 27.0 

St. Catharines 6.56 41.0 

SaultSte. Marie 5.56 40.0 

Thunder Bay 6.89 34.0 

Toronto, West 29.69 58.0 

Windsor 7.63 38.0 

Average, area weighted 10.47 38^^ 7^6 10^3 4.5 39.3 100 



Manvel et al [13] present data on land use in 106 United States 
cities. Analysis of these data indicates that the proportion of the 
urban area in each land use category was relatively similar after 
deducting the proportion of the urban area which is in the undeveloped 
category. This category is deducted from the total urban area to obtain 
the developed portion of the urban area. Figure 2, shows the percent 
undeveloped land as a function of population density for the U.S. cities 
and the percent open land for the Ontario cities. The open category 
would include undeveloped land and developed open space land, e.g., 
parks. As can be seen from Figure 2, the Ontario cities generally fall 
above the line of best fit for the U.S. undeveloped land. This result 
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FIGURE 2. PERCENT UNDEVELOPED LAND USE (U.S.) AND OPEN SPACE 
LAND USE (Ontario) VERSUS POPULATION DENSITY 

Note that best fit line is forced through 100 percent 
at PD=0 
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appears reasonable since the open space category would include the 
undeveloped land. Thus, it seemed reasonable to use the U.S. equation to 
estimate the percent of the urban area which is undeveloped, i.e., 

Z = i.0e-0-170(PD) (r = -0.57) (1) 

where: Z = proportion undeveloped land, 

PD = average gross (developed and undeveloped) population 
density, persons per acre, and 
r = correlation coefficient (-1.0 ± r <_ 1.0), 

Using this relationship, the developed population density, 

PD(i. can be expressed as a function of the gross population density 

using 

PO __1D_, PDe°-^^^" U) 

^^d (1-Z) " 0.17 PD_^ 

Equation (2) is shown in Figure 3. Note that the developed 
population density is about six persons per acre at the lowest level of 
urbanization (one person per acre). The developed population density 
approaches the gross population density as PD increases. Indeed, they 
are quite close at PD > 25 persons per acre. 

After correcting for the percent undeveloped, the proportion of 
the land in the developed uses was detennined as a percent of developed 
urban land only. After this transformation was made, the percent of land 
in the developed uses seemed to be independent of population density. 
The resultant distribution of developed land by use and undeveloped land 
is shown in Table 8. Note the similarity of the Ontario and U.S. land 
use distributions. 

The land distributions for all cities are determined using 
equation (1) and Table 8. The results are presented in Table 9. In 
determining the control costs, only the developed portion of the urban 
area is considered. Thus, it is important to check the validity of this 
assumption in future assessments. Actual field data need to be gathered 
and analyzed using a consistent set of assumptions regarding land use 
categories. 
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FIGURE 3. RELATIONSHIP BETWEEN GROSS AND DEVELOPED POPULATION DENSITY 

3.3 Population and Area Served by Type of Sewer System 

The area served by combined sewers was determined by on-site 
interviews In the nine cities, a survey by Waller [14], and questionnaire 
results fron a 1966 survey [15], These three sources provided estimates 
for 49 of the 56 cities. The remaining seven cities were assumed to have 
the same percentage of combined sewers as the other 49 cities, which is 
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TABLE 8. DISTRIBUTION OF DEVELOPED LAND USES IN 
ONTARIO TEST CITIES AND U.S. CITIES 

Percent of Total 
Land Use Ontario U.S.*^ 

Residential 52.5 58.4 

Commercial 10.3 8.6 

Industrial 14.0 14.8 

Other^ 23.2 18.2 

Recreational, schools and colleges, and cemeteries 

Heaney, J. P., W.C. Huber, et al . , Nationwide Evaluation of 
Combined Sewer Overflows and Stormwater Discharges: Volume II, 
Cost Assessment and Impacts, U.S. EPA 600/2-77-064, 1977. 

very close to 25 percent. Unfortunately, much of the data expressed the 
combined sewered areas as a percent of total sewered area^ which is 
unknown. Also, no data were available regarding population served by the 
combined sewer systems. Thus, indirect estimating procedures were used as 
discussed below. 

Information was obtained regarding the population density 
distributions of the nine test cities. Using these data, the census 
tracts were ranked by population density and grouped into five categories 
ranging from lowest density to highest density as shown in Table 10. A 
histogram for the city of Thunder Bay is shown in Figure 4. 

An equation of the form 

PD = axb (3) 

where: PD = gross population density, persons per acre (PD >^ u), 
X = percent of urban area (0 j^ x £ 100), and 
a,b = parameters, 

was fit to these data. The average population density in any interval, 
^1 to X2. is 

PD = . L_ / ^ ax^dx. ^^) 
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TABLE 9. LAND USE DISTRIBUTIONS FOR THE ONTARIO CITIES 
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TABLE 10. POPULATION DENSITY DISTRIBUTIONS IN NINE ONTARIO CITIES 



Urbanized Area 



% of Land Area 
Population Density Group (PDG) 
II III IV V 



Population Densities, 
persons/acre in PDG 
II III IV 
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Burlington 


36.2 


13.3 
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1.71 


3.82 


5.34 


8.98 


13.59 


Guelph 


35.4 


15.5 
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18.3 


15,5 


0.70 


2.63 


3.52 


8.81 


14.17 


Kingston 


37.1 
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13.2 


15.7 


4.28 


7,15 


11.60 


18.22 


25.92 


Kitchener-Waterloo 


37.6 
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17.4 


14.8 


12.5 


1,89 


4.65 


8,06 


11.67 


15.82 


St. Catharines 


32.6 


22.3 
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17.3 


3.33 


5.08 


6.53 


7.95 


13.37 


Sault Ste. Marie 
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1.34 


2.88 


5.32 


9.02 


12.59 


Thunder Bay 


33.8 
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2.45 


4.02 
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14.68 


Toronto 


32.9 
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8.56 
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36.95 


47.23 


74.21 


Windsor 


31,5 
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1.74 


4.71 


7.63 


11.13 


19.26 
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FIGURE 4. POPULATION DENSITY DISTRIBUTION OF THUMDER BAY, ONTARIO 
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To calibrate the overall average population density with the calculated 

popuTation density, an approximate value of x, was found such that 

PD lAn = PD 1 For example, for Thunder Bay, 

x, - 100 calc, 

PD = 86.4 x'^''"'^^ (5) 

and the actual average population density is 6,89 persons per acre (17.0 

persons/ha). 

Thus, 

^•*^^=T00^ x/'™ 84.6 x-"-^^ldx (6) 

This equation is solved to find X] , To calibrate the overall average 
population density with the calculated population density, for the other 
47 cities, values of b were assigned based on their similarity to the nine 
test cities. A value of x, =2 was used to avoid instability problems. 
Then, a was calculated such that the average population density, PD was 

PD. .. = 1 f^^^ ax'^dx (7) 



f 100 b 

'2 - 100 T7r;r= 77~ ^ J 



100 - 2 
_ 98 PD (1+b) (o^ 

^ = ^^^nTbTTTHbi- 

Thus, the final equation for gross population density is 

PD = ax with x^ < X < 100 (9) 

Given the equation in the form PD = ax , one can find the average 
population density, the proportion of the population within certain 
densities and so forth. 

The population density function, PD = ax , is given in tenms 
of the total urban area. Thus, it needs to be modified to integrate over 
only the developed portion of the urban area as shown in Figure 5. In 
order for the area under the two curves to be equal , one must have 

100 ^IOO(I-Z) 

a- xb dx (10) 



] ax^dx= ] 
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FIGURE 5. CHARACTERIZATION OF POPULATION DENSITY IN URBAN AREAS 



or 



Then, 



PD. = a'x'^ 112: 

d 



where: PD. = population density in developed portion of the urban area, 
a' = adjusted coefficient from equation (11), and 
x-j = calibrated lower limit on percent urbanized area. 

The percent of the urban area which is sewered is known for the 
nine test cities. Computing the corresponding PD . for seven of the 
cities resulted in the values shown in Table 11. Guelph and Toronto were 
considered extreme values and not entered into Table 11. Based on these 
data, a cutoff marginal developed population density of five persons per 
acre (12.4 persons/ha) was used to delineate the sewered part of the 
urban area. Solving equation (12) for x« yields: 

X2 = min [(5/a)^^,100(l-Z)] (13) 

where: X2 = percent of the urban area which is sewered. 

Knowing the percent of the urban area which is undeveloped, 
i.e., lOO(Z), the combined sewered area x , from the survey data, and 
the percent of the urban area which is sewered, x^, then the other 
sewered and unsewered developed areas can be calculated as residuals. 
The calculation procedure is summarized below: 

Sewered Areas As a Percentage of Total Urban Area 

1) Undeveloped Land = lOO(Z) = x^ 

2) Sewered Area = X2 

3) Combined Sewer Area = a(x2), where a is the proportion of the 
sewered area which is combined. Values of a are taken from survey 
data for 49 cities. For the remaining seven cities, it is assumed 
that a = 0.25. 

4) Storm Sewer Area = (1 - a)x2 

5) Unsewered Developed Area = 100 - Xy - X2. 
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TABLE 11. MARGINAL POPULATION DENSITY FOR SEWERED 
PORTION OF SEVEN URBAN AREAS IN CANADA 

Marginal Sewered Population 
Density, ?U^ 
City (persons/acre) (persons/ha 

Burlington 3.93 

Kingston 7.28 

Kitchener 4.49 

St. Catharines 5.79 

Sault St. Marie 4.67 

Thunder Bay 6.12 

Windsor 3.87 

Average of 7 cities 5.16 12.8 



The results of these calculations are shown in Table 12. 

The population served by the various types of sewerage systems was 
estimated indirectly since data were unavailable. Combined sewers are 
installed in the older and typically more densely populated portions of 
the urban areas. An earlier study by Waller of combined sewers in Canada 
indicated that the average population density in the combined sewered 
portion of the two largest cities was 30.7 persons per acre [14], The 
average density in the other 30 cities reporting combined sewers was 14.7 
persons per acre. These results suggested that a good approximation of 
the population served by combined sewers would be obtained by assuming 
that the high density areas were combined. The next highest density is 
served by storm sewers and the lowest densities are unsewered. Thus, 
the population served by type of sewer system is calculated as shown 
below: 

Population Served by Type of Sewerage System 

1 ) Combined Sewers: 

P = JL f '^ a' x^d, (14) 

^ 100 x^ -^ 
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TABLE 12. LAND USE BY TYPE OF SEWERAGE SYSTEM 
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2) Stonn Sewers: 



X, 



^. = A. f ^ a' x^ dx (15 

^ lUO X. -^ 



c 
3) Unsewered (no storm or combined sewers 

X 



P = A [ 

" Too xj 



^ a' X dx (16 



^2 

where: A = total urban area in acres 

The resulting population by type of sewerage system is shown in Table 
13. Lastly, the population densities by type of sewer system are shown 
in Table 14, 

3.4 List of Variables 

a constant 

a' adjusted coefficient 

■^■■- proportion of sewered area which is combined 

I? constant 

Pc population served by combined sewer system (persons) 

Ps population served by storm sewer system (persons) 

Pu unsewered population 

PD gross population density (persons per acre) 

PD average population density (persons per acre) 

P^calc calculated average population density (persons per acre) 

average population density in interval from x^ 

to X2 (persons per acre) 
PDd population density in developed portion of urban area 

(persons per acre) 
I* correlation coefficient 

U percent of urban area 

X"! calibrated lower limit on x such that average PD 

corresponds to the integrated average PD 
Xf calibrated upper limit on x such that average PD 

corresponds to the integrated average PD and percent of urban 

area which is sewered 



PD 



^1 " ^2 
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TABLE 14. POPULATION DENSITY BY lYPE OF SEWERAGE SYSTEM 
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Xc percent of urban area served by combined sewers 

^d percent of urban area which is developed 

\ percent of urban area which is undeveloped 

2 undeveloped portion of urban area 



4* RUNOFF AND POLLUTANT LOAD ESTIMATES 

The purpose of this Section is to estimate the quantity and 
quality of urban runoff from the 56 urban areas. The first part provides 
some background information regarding models which are used to assist in 
making decisions. Then, precipitation patterns are analyzed to form a 
basis for predicting the quantity of urban runoff. The relevant water 
quality parameters are discussed and the results of numerous attempts to 
estimate runoff quality are presented. Finally, a pollutant load 
prediction equation is developed which provides the basis for assessing 
pollutant loads. A summary of the methodology is presented elsewhere 
[62, 63], 

4.1 Modeninq of Urban Runoff 

4.1.1 Computer models 

The overall goal of urban runoff modelling is to aid in 
decision-making for the abatement of water quantity and quality problems. 
Thus, computer models do not provide "solutions" to problems. Rather, 
they serve as useful tools to those charged with devising such solutions. 
Within this context, subobjectives of the modelling process may be 
identified: planning, design, and operation. Models for the latter 
category are generally site-specific [16, 17] and were not considered 
during this study. However, numerous models are available for planning 
and design purposes, e.g. the U.S. Army Corps of Engineers' STORM and the 
U.S. EPA Storm Water Management Model (SWMM). However, they are not 
unique; several other urban runoff models are capable of similar tasks 
[18, 19. 20]. 

Computer models are merely mathematical abstractions of the 
physics of the urban runoff processes and do not necessarily produce 
accurate or even logical predictions without extensive calibration/ 
verification data. These are in addition to data required as model 
input, such as topography, land use, rainfall, antecedent conditions, 
description of drainage system and storage-treatment facilities. 

Among the principal 1965 findings of the ASCE Council on Urban 
Water Resources [21] was the serious need for field data on rainfall- 
runoff-quality for several catchments. Unfortunately^ in spite of efforts 
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by federal agencies in the United States and Canada, this need still 
exists, especially in the area of data for calibration/verification of 
urban hydrology models. For instance, few new discharge data have been 
acquired on sewered catchments, using flumes, weirs, or other accurate 
devices rather than conversion of stage gauge readings for determination 
of flow rates [22], The state-of-the-art in computation and simulation 
tools has, thus, outstripped its available calibration/verification data 
base. However, current analytical techniques will necessarily rely on 
computer models, especially for planning and design purposes. Hence, the 
ultimate goal of acquisition of salient field data remains worthwhile and 
necessary. Throughout this section, gaps in available data for input and 
calibration/verification will be apparent. But the useful analyses which 
can still be performed without these data should also be clear. 

The modelling procedures developed for the assessment wil 1 be 
discussed in detail. Two levels of sophistication are considered: use 
of STORM for the development of the parameters used in the assessment 
methodology described in Section 6. and use of a very simple runoff 
prediction technique for the 56 urban areas of the Ontario assessment 
itself. 

4.1.2 Runoff analysis using STORM 

The Storage, Treatment, Overflow and Runoff Model (STORM) was 
developed by Water Resources Engineers, Inc., (WRE) for the Hydrologic 
Engineering Center (HEC) of the U.S. Army Corps of Engineers [11, 23]. 
The model was designed for planning purposes, i.e., for long-term simula- 
tion of many storm events using an hourly time step. For instance, the 
model has been used to simulate runoff quality and simple storage-treat- 
ment options for a 63-year record of hourly rainfalls in San Francisco. 

When STORM is used for planning, the primary objective is to 
illustrate the effect of various storage-treatment combinations at the 
downstream end of an entire urban catchment that provide given levels of 
control, "Level of control" may refer to percent of runoff captured, 
percent BOD or other pollutant removed, number of overflows per year, 
quantity of overflow per year, etc. Use of the model for this primary 
objective is described in detail in Section 6 of this report, including a 
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discussion of the methodology employed. Thus, the use of STORM in this 
study is deferred to that section. 

4.2 Runoff Prediction for Ontario Assessment 

4.2.1 The hydrologic cycle 

The hydrologic cycle may be divided into three principal 
phases: 1) precipitation, 2) evaporation, and 3) surface and groundwater 
runoff. The hydrologic cycle has neither beginning nor end, as water 
evaporates to the atmosphere from land and water surfaces. The 
evaporated moisture eventually precipitates back to the earth where it 
may be intercepted or transpired by plants, may become surface runoff, or 
may infiltrate into the ground. Once in the ground, water may be stored 
as soil moisture and evapotranspired, or percolate to deeper zones to 
become part of groundwater flow. Surface and groundwater flow from the 
land eventually reaches streams, lakes, or oceans from which water 
evaporates to complete the cycle. 

4.2.2 Form of the equation 

Techniques for prediction of runoff quantities vary from ^/ery 
simple methods of the rational method type to sophisticated models of the 
nature of SWMM. The technique used in STORM is relatively simple, rely- 
ing on weighted average runoff coefficients and a simple loss function to 
predict hourly runoff volumes. Nonetheless, because of the nature of the 
continuous simulation involved, it is at a considerably higher level, and 
therefore more complex, than earlier, desk-top techniques. 

Due to the complexities and data requirements of STORM, it was 
not possible to run the model on all cities of the assessment, or even a 
majority. Rather it was run only on four test cities as discussed in 
Section 6. However, in its limited application, useful information was 
learned regarding formulation of a simple runoff prediction method for 
application to all the cities of the assessment. 

Runoff is a function of meteorologic, hydrologic, topographic, 
and demographic factors. On an annual basis, many of the factors may be 
considered constant, so that runoff is predicted on the basis of 
differences between areas rather than reflecting seasonal variations 
within a year. Hence, the prime meteorologic and hyrologic factor is 
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annual precipitation, and other factors are incorporated into a 
conversion to annual runoff. 

These considerations lead directly to the use of a simple runoff 
coefficient method in which runoff is merely a fraction of rainfall. 
This approach has been used successfully by Miller and Viessman [24] for 
runoff prediction on an individual storm basis in urban areas. This 
equation was: 

AR = 1.165 (I - 0.17) (P - la) (17) 

where: AR = runoff, inches; 

I = fraction imperviousness; 

P = precipitation, inches; and 

la - initial abstraction, inches. 

The recommended value of Ig, which accounts for initial losses such 
as depression storage, interception, etc., was between 0.10 and 0.15 in. 
(0.25 and 0.38 cm), and the equation was deemed valid for a range of 
imperviousness between 35 and 80 percent. Extrapolation for use on an 
annual average basis, however, may be questionable, particularly in the 
matter of how much water should be abstracted out of the cycle on an 
annual basis. Hence, an equation will be used that is similar in form 
to equation (17), but which is consistent with the STORM simulation 
runs, described in Section 6, on which the overall assessment was based. 

STORM computes a runoff coefficient, CR, weighted between 
pervious and impervious areas by: 

CR = 0.15 (1 - I) + 0.90 I 

= 0.15 + 0.75 I (18) 

where I is fraction imperviousness and the coefficients 0.15 and 0.90 are 
the default values used in STORM for runoff coefficients from pervious 
and impervious areas, respectively. Note that in both equations (17) and 
(18) the effect of demographic factors (e.g., land use, population 
density) is incorporated into the imperviousness, I. 

Imperviousness was estimated for the nine cities discussed in 
Section 4 using the same maps used to estimate land use. The average 
residence was assumed to have an impervious area of 1,500 ft^ (139 m^) 
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and the average street width was taken to be 34 ft (10.4 iii) regardless of 
land use. These assumptions were necessary due to the large scale of the 
aerial photos. Structures on institutional, industrial, and commercial 
lands were usually large enough to be measured. For each city, the 
procedure was to choose three or four representative areas for each land 
use and to determine a figure for percent imperviousness from these 
samples. The determination of curb length for each land use also 
required the use of aerial photographs. Again, representative areas of 
each land use were chosen. Street lengths were measured, and along with 
the assumption that curb length equals twice the street length, figures 
for curb length per acre were determined. The results of these studies 
are shown in Table 15. 

The American Public Works Association [3], Graham et al [25], 
and Stankowski [26] have developed equations to predict imperviousness as 
a function of population density. The imperviousness is to be estimated 
for the developed portion of the urbanized area only. The weighted 
average Imperviousness and population density were calculated for the 
nine Ontario cities. These results are plotted on Figure 6, along with 
the three estimating curves. Also, a tabulation was made of the 
imperviousness due to streets alone for various block sizes as shown in 
Table 16. These results are also plotted on Figure 6. A comparison of 
these various plots and the actual data indicates that the New Jersey 
[26] equation provides a suitable predictive equation with population 
density defined as developed population density. Thus, the equation used 
to estimate imperviousness is 

I = 0.96 PD, (0-"3-0.0391 log^^PD^) ,,,, 

where: I = imperviousness in percent, and 

PDd = population density in developed portion of the urbanized 
area, persons per acre^ 

The simplified equation for estimating annual runoff (AR) is now 

AR = {0.15 + 0.75I)P (20) 
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TABLE 15. IMPERVIOUSNESS AND CURB LENGTH DENSITY BY TYPE OF LAND USE IN NINE ONTARIO CITIES 





Urbanized Area 


Impe 
or 
De 


rviousness, I* 
Curb Length 
nsity. 6l** 


Residential 


Commercial 


Industrial 


Institutional 


Open 




Burlington 




I 


32 
365 (275) 


89 
178 (134) 


11 
136 (102) 


37 
125 ( 94) 


3 

57 ( 39) 




Guelph 




1 


30 
339 (255) 


89 
254 (191) 


43 
121 ( 91) 


36 
153 (115) 


1 
32 ( 24) 




Kingston 




1 


27 
332 (250) 


87 
314 (236) 


20 
82 ( 62) 


17 
56 ( 42) 


44 ( 33) 




Ki tchener-Waterl oo 




J 


29 
355 (267) 


82 
216 (163) 


52 
142 (107) 


36 
113 ( 85) 


8 
35 ( 26) 




St. Catharines 




1 


27 
331 (249) 


90 
238 (179) 


42 
146 (110) 


38 
150 (113) 


4 
60 ( 45) 




Sault Ste. Marie 




1 


22 
353 (266) 


n 

461 (347) 


57 
150 (113) 


22 
244 (184) 


3 

44 ( 33) 




Thunder Bay 




I 


29 
388 (292) 


78 
331 (249) 


44 
58 ( 44) 


32 
216 (163) 


2 
41 ( 31) 




Toronto 






44 
355 (267) 


52 
362 (273) 


44 
261 (197) 


31 
307 (231) 


14 
345 (260) 




Windsor 




I 
6l 


31 
382 (288) 


88 
337 (254) 


48 
121 ( 91) 


18 
133 (100) 


S 

71 ( 53) 




Average 




I 


30 
356 (268) 


81 
299 (225) 


40 
135 (102) 


30 
166 (125) 


81 ( 61) 



* % measured in percent of total 
** ft per acre (metres per hectare). 
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FIGURE 6. IMPERVIOUSNESS AS A FUNCTION OF POPULATION DENSITY 



54 



TABLE 16. EFFECT OF URBAN BLOCK SIZE ON CURB LENGTH 
DENSITY AND IMPERVIOUSNESS DUE TO STREETS 



Block Size Curb Length Density Imperviousness due 

ft X ft {m X n) Area, ac (ha) ft/ac (m/ha) to Street* 



660 X 330 5 { 2.02) 392.0 (298,0) 0.150 

{ 201 X 101) 

1,320 X 660 20 ( 8.09) 198.0 (148.0) 0.077 

( 4U2 X 2011 

2,640 X 1,320 80 { 32.40) 99.0 ( 74.6) 0.039 

( 805 X 402) 

5,280 X 2,640 320 (130.00) 49.5 ( 37.3) 0.019 

(1,609 X 807) 



* Assume 34 ft (10.4 m) wide street. 

where: AR = annual runoff, inches, 

I = fraction imperviousness from equation (19), and 
P = annual precipitation, inches. 

A comparison of STORM simulated runoff versus calculated runoff using 
equation (20) is shown in Table 17, The average difference is about 0.5 
Inches (1.27 cm) per year, A similar comparison in the U,S. assessment 
indicated a difference of 0.3 inches (0.76 cm) per year. Thus, a 
correction factor was added to equation (20) to reflect this difference. 
The final equation is 

AR = (0.15 + 0.751) P - 0,5 (21) 

Based on equation (21), wet-weather flow estimates were made for the 56 
urban areas for the combined, storm, and unsewered areas. The results 
are shown in Table 18. 

4.2.3 Dry-weather flow prediction 

Dry-weather flow is predicted, based on actual flow data for 
the test cities, indicating an average flow of 108 U.S. gallons per 
person-day (410 litres per person-day). Upon multiplication by population 
density and conversion to appropriate units: 
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TABLE 17. COMPARISON OF SIMULATED AND CALCULATED RUNOFF FOR FOUR TEST CITIES 







Annual 
Precipitation 
in/yr (cm/yr) 


Runoff'^ 
Coefficient 


Ri 


inoff: in/yr 


(cm/yr) 


City 


STORM^ 


Calculated 


Difference 


Burlington 




32.4 
(82.3) 


0.344 


10.62 
(26.97) 


11.15 
(28.31) 


0.53 
(1.34) 


Kingston 




37.8 
(96.0) 


0.306 


11.11 
(28.22) 


11.57 
(29.38) 


0.46 
(1.16) 


St. Cathari 


nes 


32.4 
(82.3) 


0.372 


11.41 
(28.98) 


12.05 
(30.61) 


0.64 
(1.63) 


Sault Ste. 


Marie 


36.7 
(93.2) 


0.413 


14.54 
(36.93) 


15.16 
(38.50) 


0.62 
(1.57) 



See Section 6.4. 

Equation (18) and Table 15. 

^Assumed depression storage = 0.01 inch, runoff coefficient from equation 
(18); see Section 6.4. 

^Equation (20) 

DWF = 1.45 PDjj (22) 

where: DWF = annual dry-weather flow, inches per year, and 
PD(j = developed population density, persons per acre. 

Results of these runoff calculations are shown in Table 19, Dry-weather 
flow is generated for entry into a sanitary sewer for storm or unsewered 
areas. 

Dry-weather flow and wet-weather flow for the developed portions of 
an urban area with a precipitation of 15, 30, or 45 inches per year are shown 
in Figure 7. Note that dry-weather flow predominates at higher population 
densities which have historically prevailed in cities. However, with the 
trend towards lower density urban living, wet-weather flows take on 
greater relative importance. Indeed, they are larger than dry-weather 
flows at the lower population densities. 
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TABLE 18. ANNUAL WET-WEATHER RUNOFF FOR COMBINED, 
STORM, AND UNSEWERED AREAS 



NO 



UPBANlZED A*<EA 



ANNL- 

HRfcrP 

IN/VK 



CJMB STURmIUnSFwI aver ' 



1 1 AjAy 

" AURPRA 

BARRIE 

S^LLE^LLF 
BRAMPTON 

BRAWTFORD 
BURLINGTON 
CHATHAM 
CWINGUACOUSV 

cnenufG 

DUNnAS 
ETOBICOKE 

fiEhSGETOHN 
GUELPH 

maRilt 

kingston 



ILTON 

gSt 
r 

LEAM 
CiNnSAY 



KftrM. -WATERLOO 
- MIN--- 

MARKAM 



ITCn.-WATI 
EAMINGTDN 
INHSA'- 
nNrtON 




_- NIAGRA FALLS 
27lNnRTh BAY 
dSIOAKvXLLE 

3piOSHAWA 
SilOWEW SOUN^ 



12 
33 

36 
39 
aO 
^1 

^£ 
as 
uu 
us 

46 

a? 

a? 

io 
>i 

53 
56 



PETFHBOROUO* 
PICKERING 
PORT CQLBOU''NE 
PORT ERIE 
PRESTHN 
RlC"MONr) HILL 
ST. CATHARINES 

st. thomas 
saRtna 

SLT.^STE. MARIE 
SCARBOROUGH 
STMCOE 
STRATFORD 

SMDnURY 

Twy^DEH 

TORnNTQ 

TRENTON 

WALLACEQURG 

WELL AND 

WHITBY 

WINDSOR 

HGQnSTOCK 

YHRK 



BAY 



'^'-'',m] 



YORK 



33 



n 



111 



19.1 
0.0 
0.0 

17.0 
0.0 

o.g 

0.0 
13.3 

0.0 
0.0 
0.0 
0.0 
0,0 
0.0 
0.0 
14.2 

i?:i 



u.i 

19. f 



iS:^ 



1 



4 



.8 

9 



Ot 



I 4 



WFIGHTEO AVE. 132.751 16. d 



12.B 



9,51 12. b 
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TABLE 19. ANNUAL DRY-WEATHER RUNOFF FOR COMBINED, 
STORM, AND UNSEWERED AREAS 



URBANIZED AREA 



ANNL. 
PRt'Cf 
IN/YK 



DRY HEaTHEH flow 
CiNCMgs PES YESR5 
CQ'^B STDRHJUN3EWI AVER 




AjAy 

AURORA 

BARPIE 

SELLEVILLE 
RAmPtQn 

rantforq 
urlingtdn 

HATHAM 
HING^ACOUSY 



OBDUR 
UNOAS 

tobJc 



OKE 



inqstqn 



HA 

K 

K 



GETOWN 
TON 



_ CH.-WATERtOQ 

eamington 
tndSay 



m°'- 



KAM 

NTAGRA FALLS 
NORTH BAY 
QAKVZLLE 

WEN^SOyNO 
fTERBQfiOL/tiH 



P 

PQRt 
RE8T0N 




OUHNE 



g 



st. thomas 
saAina 

'^t. ste. marie 
caSborough 



mn 



JOBURY 

HUNDER 
RDNT 



ORO 



BAY 

tqRdnto 

TREWTON 

.ACEBURG 
AND 
BY 
"OR 
*■ TOCK 
YORk 
YORK, EAST 
YORK.^^ORTM 



WALL 



l"l 



IWEISHTEn AVE. I3S.75 
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persons/ hectars 
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FIGURE 7 



10 20 30 40 50 

DEVELOPED POPULATION DENSITY, persons/acre 
COMPARATIVE MAGNITUDE OF ANNUAL WET- AND DRY-WEATHER FLOWS 



4.3 Quality Parameters 

4.3.1 Parameter definitions 

Quality analyses may be performed at several levels of detail, 
ranging from an explicit formulation of runoff quality for small 
subcatchments within a city, to a broad representation of pollutant loads 
for an entire urbanized area or province. It has been necessary to 
consider the entire spectrum during the course of this study. 

It is unfortunate that perhaps the only consistent remark about 
runoff quality analysis in general is that data and results of previous 
studies are so remarkably inconsistent. Few studies have been made of 
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characteristics of street litter, and they offer a wide range of values 
of concentrations and loads. Effluent data show a similar scatter. 
However, it is necessary that a decision be made regarding actual values 
for use in the analysis. This section will describe methods used for 
predicting runoff quality, data required for their use, and final results 
used in this study. 

Urban runoff quality may be characterized by a variety of 
parameters. However, the list is generally shortened for modelling 
purposes to those characteristic of solids, oxygen demand, health hazards 
and aquatic growth potential, as indicated in Table 20. It is 
discouraging that even at this juncture, a serious problem of definition 
of terms arises because of various possibilities for analyzing and 
reporting quality parameters. The assurance that analyses have been 
performed according to Standard Methods [27] is not enough information. 
For example, solids are sometimes reported as "residue" instead of 
solids, and "filterable residue" instead of "dissolved solids", because 
of the nature of the evaporation and filtration techniques utilized in 
the chemical analyses. Generally, "solids" and "residue" are synonymous, 
and "solids" will be used in this report. Another problem arises from 
the fact that pollutants may be in both soluble and insoluble forms. 
Some studies report concentrations of only the soluble portions of, say, 
BOD and PO4, leading to unrealistically low values if the reader 
mistakenly thinks of them as total (soluble plus insoluble) concentra- 
tions. On the other hand, it is important to know the relative soluble- 
insoluble fractions of pollutants since this has a major impact upon 
treatability. That is, pollutants that appear as suspended solids are 
relatively easy to remove (e.g., by sedimentation) compared to those that 
are soluble. 

To further complicate the picture, no clear relationship exists 
between data derived from studies of surface litter (gathered by 
sweeping, vacuuming, flushing) and those resulting from analysis of the 
runoff itself (e.g., samples of storm and combined sewerage effluent). 
Thus, a mixture of data exists, derived from both surface and effluent 
sources. However, there is no study in which samples of both types have 
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TABLE 20. TYPICAL QUALITY PARAMETERS OF URBAN RUNOFF MODELS 



Quality Characteristic 



Representative Quality 
Parameters 



Solids 



Oxygen Demand 



Health Hazards 



Aquatic Growth Potential 



Surface "Dust and Dirt" 
Surface "Solids" 
Total Solids 
Suspended Solids 
Dissolved Solids 
Volatile Solids 
Settleable Solids 

BOD. COD 

Total Organic Carbon 

Organic N, NH3 

Total Colifonns 
Fecal Colifonns 

0rtho-P04 
Total PO4 
NO2, NO3. Total N 



been gathered simultaneously. Hence, the relationship between the two 
is not well defined, and it Is difficult to draw conclusions from all data 
considered together. 

In this report, the solids relationship shown in Figure 8 was 
applied. Note that total solids (TS) is the sum of dissolved solids (DS) 
plus suspended solids (SS), and that total, dissolved, and suspended 
solids may be separated Into a volatile portion (generally considered the 
organic portion) and a fixed portion. Volatile solids (VS) will refer to 
a portion of total solids in this report, unless otherwise Indicated. 
Settleable solids are some fraction of suspended solids. Note, finally, 
that an upper limit on the size of total solids reported is imposed by the 
size of the openings in the sampling equipment {e.g., a quarter inch- 
screen). 

Similar diagrams may be prepared for nitrogen and phosphorus, as 
shown in Figure 9, and Figure 10. For these parameters, it is necessary 
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TOTAL SOLIDS {BY EVAPORATION) 



FILTER (TYPICALLY .45m) 



DISSOLVED SOLIDS 

/\ 



DRIED AT 60°C 



SUSPENDED SOLIDS - - - ^ 

I 
I 



/ I 



FIXED VOLATILE VOLATILE FIXED 

(DRIVEN OFF 
AT 600° C) 



SETTLEABLE SOLIDS 



FIGURE 8. RELATIONSHIPS AMONG SOLIDS PARAMETERS 



TOTAL N 



NO 2 NO 3 



KJELDAHL N 



j. 



ORGANIC N 



1 

NH 



FIGURE 9. RELATIONSHIPS AMONG NITROGEN 
PARAMETERS 



TOTAL PO. 



ORTHO-PO4 
(HYDROLIZED) 



OTHER 



FIGURE 10. RELATIONSHIPS AMONG PHOSPHORUS 
PARAMETERS 
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to know whether concentrations are being reported of the element itself 
(e.g., phosphorus) or the compound (e.g., PO4), although conversions 
can readily be made on the basis of the molecular weight of each. 
Regarding the nitrogen relationships, all concentrations should be 
reported in terms of N (i.e., NOg-N, NO3-N, NH3-N) in order for 
mass balances to be performed easily. 

4.3.2 Parameters for assessment 

For the Ontario assessment, five parameters were used that are 
representative to some degree of the quality characteristics indicated in 
Table 20. These are indicated in Table 21. 

TABLE 21. QUALITY PARAMETERS USED IN ONTARIO ASSESSMENT 
Parameter Abbreviation 

1. Five-Day Biochemical Oxygen Demand BOD5 or BOD 

2. Suspended Solids SS 

3. Volatile Solids VS 

4. Total Phosphate (as PO4) PO4 or TPO4 

5. Total Nitrogen (as N) N 

Note: All parameters (except suspended solids) are totals that include 
dissolved and insoluble portions, and are usually determined as 
in Standard Methods [27], All are usually reported in concentra- 
tion units of mg/L (equivalent to ppm) . 

Five-day BOD was used because of its broad acceptance and 
traditional role in water quality analysis. Its usefulness is severely 
impaired by the great difficulty in performing accurate and consistent 
laboratory analyses. For instance, there is no standard for laboratory 
comparison, and low-level values (e.g., 10 mg/L) are especially 
susceptible to errors of up to 100 percent. Moreover, studies have shown 
that results are affected by the percent dilution and are generally not 
reproducible [28]. In addition, samples are affected by amounts of heavy 
metals and other parameters present. Use of COD and/or TOC avoids some of 
these problems, but their relationship with traditional stream sanitation 
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analysis (i.e., prediction of dissolved oxygen) is unclear, and most 

people are used to thinking in terms of BOD. It is used in this study, 
realizing its limitations. 

The other four parameters were used because of general acceptance 

and availability of data. It should be borne In mind that many options 
are available for modelling purposes, and the choice of parameters is 
somewhat arbitrary. 

4.4 Quality Prediction Techniques 

4.4.1 Pollutant loads 

The quality prediction techniques found in most urban runoff 
models (e.g., SWMM, STORM) rely upon generation of an Initial surface load 
of pollutants. This load is usually expressed in units of lbs, lbs/acre, 
Ibs/curb-mile, Ibs/day-acre, or Ibs/day-curb-mile (or equivalent metric 
units). Normalized loads are, of course, multiplied by a unit of area, 
dry days, etc., to produce an initial mass of pollutants at the start of 
the storm. Pollutants are then "washed off" during a storm in an exponen- 
tial fashion in which the amount removed per time step is proportional to 
the amount present, the runoff rate, and other factors. SWMM [10, 29] and 
STORM [11, 23] documentation contain details of this methodology. The key 
factor in prediction of long-term (e.g., annual) pollutant loads from 
urban areas is, however, the surface loading rates themselves, and most of 
the following discussion will be devoted to them. 

Surface loadings are usually predicted by one of two means: 
estimates based on surface accumulation data, or estimates based on measure- 
ments of effluent concentrations and flows. As mentioned earlier, no one 
study has performed the analysis both ways, so comparisons are not easily 
accomplished. However, to obtain the study objective, normalization of 
loadings rates by some means that could be converted to total mass of 
pollutants upon multiplication by area, days, and/or other appropriate 
parameters was necessary. As a result, both methods were utilized in the 
developments that follow. 

4.4.2 Surface accumulation methods 

Both SWMM and STORM use this method for prediction of the total 
soluble mass of pollutants (except for solids) available at the beginning 
of a storm. For suspended or settleable solids calculations, the total 
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mass is given since there is no "soluble" portion. The method is based 
upon the following equation, given in representative English units: 

PiJ - ddi • Fi,j • Gl.I • Ai . Nd+ Po (23) 

where: P^- j = total soluble pounds of pollutant p on urban land use 

i at the beginning of the storm; 
dd^- = pounds of accumulated dust and dirt (or "surface 

solids") per curb mile per dry day; 
F^- J = total soluble pounds of pollutant p per pound of dust 

and dirt found on land use i; 
G|_ -j = number of curb-miles per acre of land use i; 
A^ = area of land use i, acres; 
Nq = number of dry days since last storm; and 
Pq = total soluble pounds of pollutant remaining on land 

use i at end of last storm. 

The dust and dirt accumulation rate is often given in terms of pounds/day 
per 100 feet of curb instead of curb-miles, but the latter units are used 
here for ease in cocnparison with other portions of the report. 

The parameter Nq is the number of dry days since the last 
storm, not the number of days since the last storm or street cleaning 
operation. This is because, in most cases the interval between storms is 
less than the street cleaning interval. The latter is generally in the 
order of several 10 's of days and the efficiency of street cleaning 
operations is uncertain in any event. 

The parameters ddi_, Fp^|_, and Gl are functions of land use, 
L. The dust and dirt loadings, ddL, and pollutant fractions, Fp [_, data 
are shown in Table 22 [31], These parameters may be updated to some 
degree, as will be shown. As indicated in Table 22, SWMM assumes that all 
dust and dirt will pass through a quarter-inch (6 mm) screen and is insolu- 
ble, thereby appearing as suspended solids (i.e., the SS fraction, F^ ■; 
is 1.0). STORM assumes that only from seven to 17 percent of dust and 
dirt meets these requirements. Both models assume that settleable solids 
are about ten percent of suspended. The SWMM assumptions imply that the 
total of all pollutants is slightly greater than 100 percent of dust and 





TABLE 22. 


PARAMETERS 


FOR SURFAC 


E POLLUTANT 










ACCUMULATION USED IN 


SWMM AND/OR 


STORM 












Land Use 










Single-family 


Multi-family 








Parameter 


Units 


Residential 


Residential 


Commercial 


Industrial 


Open" 


Dust and dirt 














loading, dd. 


lb/day -curb-mile 


40.0 


121.0 


174.0 


243.0 


79.2 




kg/day-curb-km 


11.4 


34.4 


49.4 


69.0 


22.5 


Pollutant frac- 














tions^, F,_ J 














SS^ (SWMM) 




1.0 


1.0 


1^ 


1.0 


1.0 


SS' (STORM) 




0.111 


0.08 


047 


0.067 


0.111 


Settleable Solids 














(SWMM) 




0.1 


0.1 


0.1 


0.1 


0.1 


Settleable Solids 














(STORM) 




0.011 


0.008 


0.017 


0.007 


0.011 


BODs 




0.005 


0.0036 


0.0077 


0.003 


0.005 


COD 




0.04 


0.04 


0.039 


0.04 


0.02 


Total PO4 




0.00005 


0.00005 


0.00007 


0.00003 


0.00001 


Total N 




0.00048 


0.00061 


0.00041 


0.00043 


0.00005 


Grease 




0.001 


0.001 


0.001 


0.001 


0.001 


Total Conforms 


MPN/g 


1.3 X 10^ 


2.7 X 10* 


1.7x10* 


1 .0 X 1 0* 


0.00 



Except as noted, values are for soluble portion and derived from the 
1969 APWA Chicago study [31]. 

^All values assumed. -, -^ \ 

Fraction refers only to soluble fraction of dust and dirt (except solids). 
*^A11 values assumed at 10% of value for SS. 

dirt, while the STORM assumptions imply that the total of all pollutants 
is only about 12 percent of dust and dirt. 

The insoluble portion of pollutants is accounted for {in SWMM and 
STORM) by addition of a fraction of the solids concentration to predicted 
effluent concentrations (of the soluble portion). For example, SWMM adds 
five percent of the SS concentration to the soluble BOD concentration to 
obtain total BOD, on the basis of calibration of the original SWMM in San 
Francisco. This is because of the reliance upon the 1969 APWA Chicago 
data in which only soluble fractions were reported. It is obvious that 
equation (23) could be used to predict the total (soluble plus insoluble) 
mass of surface pollutant accumulation simply by a redefinition of terms 
(and use of appropriate revised numbers). This would facilitate quality 
calibration of the models and probably be as accurate, considering the 
available data. Final surface pollutant loads derived subsequently will 
refer to total pounds of pollutants. 
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Starting with the Chicago study and followed subsequently by 
others, it has become customary to report data in terms of mass of 
pollutants per unit length of curb, on the assumption that the curbs and 
gutters represent the main source area for acquisition of pollutants by 
the storm runoff. In order to obtain loadings on a unit area basis, it 
is necessary to obtain the length of curb per area for each land use, 
thus defining the parameter G. . in equation (23), 

L , 1 

It is expected that G. . would be a function of land use, 
which in turn is a function of population density, PD. Curb length 
(taken as twice street length) was related to population density in the 
Washington, D.C., area by Graham et al [25]. Their data, augmented by 
data from other parts of the United States by APWA [3], resulting in: 

Gj_ = 0.0782 - 0.0668 • 0.839 ^^d (24) 

where: G = curb length per area, mile/acre, and 

PD = developed population density, persons/acre. 

Equation (24) seems to work well for residential areas, but the curb 
length concept is troublesome when one is evaluating commercial, indus- 
trial, or open areas. For example, what is the equivalent curb length of 
a shopping center? Data from other sources are compared in Table 23. 

An average of the Tulsa [32] and Ontario data was used in the 
analysis. Ontario data were not used by themselves because of possible 
differences in measurement techniques which may affect relationships with 
pollutant loadings that follow. In addition, it will be seen that ratios 
of curb lengths between different land uses are the important functions; 
these change little between cities. Specific data for residential areas 
were used in lieu of equation (24), since the equation was developed to 
predict curb length as a function of population density average over all 
land uses. However, the equation may be used when considering an overall 
urban area. 

To summarize, the surface accumulation methods are convenient 
for modelling purposes and illustrate the linkages between various 
causative factors. The key missing factor is a link between the surface 
loads and effluent loads that has been verified by measurements of both. 
Until this is accomplished, such a link must be hypothesized in its 



TABLE 23. MEASURED CURB LENGTHS FOR VARIOUS LAND USES 





mile/ 
acre 


Tuls^3] 

100 ft/ 
km/ha acre 


10 Ontario Cities^ 
mite/ 100 ft/ 
acre km/ha acre 


Average of Two 

Locations for Use In Study 

mile/ 100 ft/ 
acre km/ha acre 


Residential 


0.076 


0.30 4.01 


0.067 


0.27 


3.56 


0.072 


0.29 


3.78 


Commercial 


0.081 


0.32 4.28 


0.057 


0.23 


2.99 


0.069 


0.28 


3.64 


Industrial 


0.042 


0.17 2.22 


0.026 


0.10 


1.35 


0.034 


0.14 


1.78 


Park 


0.042 


0.17 2.22 


— 


— 


- 


— 


— 


— 


Open 
Institutional 


0.016 


0.063 0.84 


0.015 
0.030 


0.059 
0.12 


0.81 
1.66 


0.024^' 


0.097" 


1.29" 



^Average of daia collected by University of Florida, 1975. Guelph, Kingston. Kitchener- Waterloo, 
Milton, St. Catharines, Sault Ste. Marie, Thunder Bay, Toronto, Windsor. 

^A verage of open plus park 



mathematical formulation, as done in SWMM and STORM. However, equation 
(23) is used in developments that follow to relate loadings between 
different land uses and pollutants, hence the reason for the previous 
developments. The other side of the coin, that is, results derivable from 
effluent data alone, will be discussed next. 

4.4.3 Effluent concentration methods 

Many studies in recent years have reported measured concentra- 
tions of pollutants in stonn and combined sewer discharges. If the flow 
rate is also known, the mass flow pollutograph may be determined (e.g., 
Ibs/min of BOD) and integrated to produce the total mass emission for the 
storm discharge. When distributed over the area of the catchment and 
divided by the number of dry days preceding, normalized loadings (e.g., 
mass-BOD/area-day) may be determined. Some studies report these values 
directly, while others report a lesser amount of information. In general, 
the surface loading may be deduced from a measured flow-weighted average 
concentration and assumed runoff quantity: 



M = P . C . CR • p 

where: M = surface loading, mass/area-time; 
P = precipitation, depth/time; 

C = average concentration = mass pollutant per mass of 
of total sample; 



(25) 
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CR = runoff coefficient; and 
p = water density, mass/volume 

For an individual storm, preceded by N^ dry days, the total depth of 
precipitation, P5, may be given. Then, 

No 

For annual average computations it may be assumed that, on an 
average basis, 

Pg = P/n t^?) 

and 

Nd = 365/n ti^«) 

where: P = average annual precipitation, depth/year; 
Nq = average number of days between storms; and 
n = average number of storms per year. 

Equation (25) may thus be used to compute average annual values since it 
results from substitution of equations (27) and (28) into equation (26). 

Equation (25) may be converted to convenient units. For instance, 

M( 5-^—) = P(i^) . C(~^) ■ CR • 62.4(1^ 

day-acre yr 10 lb ft 

43.560 2 ft . yr 

acre ^ ' 12 in. 365 day 



ar 



M = 6.206 X 10-4 . P . C • CR 1 i^y) 



where: M = average surface loading, Ib/day-acre; 
P = annual precipitation, in/yr; 
C = pollutant concentration, mg/L or ppm; and 
CR = runoff coefficient, fraction. 

Use of equation (29) presents several difficulties. It is 
inherently an average, and is susceptible to the assumptions of equations 
(27) and (28). It requires the use of a flow-weighted average concentra- 
tion. Unfortunately, such values are seldom reported in the literature, 
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if indeed any specification is made as to the types of "average" concen- 
tration presented. Runoff is generated by the simplest of methods, that 
of a runoff coefficient, with all of its well-documented errors., 

On the other hand, measured concentrations do in fact represent 
the real amount of pollutants being discharged, and thus incorporate an 
of the unknown factors involved in trying to generate surface loads 
coupled with a wash-off and transport mechanism. These include such 
factors as dust fall, air pollution, and several others not specifically 
addressed in this study. Furthermore, for purposes of the assessment 
performed in this study, very simple methods of runoff and quality 
generation must be employed. Hence, equation (29) Is consistent with 
other levels of analysis used in this research. 

In the same manner that surface accumulations could be considered 
functions of population density and land use, so can surface loadings 
derived from effluent data. In particular, both the concentration and 
runoff coefficient are clearly such functions; the latter has been 
discussed previously. In order to ascertain the functional relationship 
between the surface loadings and population density, available data for 
residential areas for which population density is given have been tabulated 
Derived surface loadings are given in Table 24. The cities included in 
the table all had data for residential areas for which population density 
was specified and from which surface loadings could be derived. The list 
is not meant to be exclusive, but represents data that were readily 
available during the study. 

The vast disparity among all the data may be seen in Figure 11. 
Both separate and combined loadings vary by more than an order of magnitude. 
Unfortunately, the variation persists if normalized by dividing by annual 
precipitation (not shown). Three cities produced wery high results 
compared to the bulk of the data: Atlanta, Bucyrus, and Durham. The 
reason for this is primarily variation from strictly residential land use. 
In addition, the open channels sampled in Durham had characteristics of 
open sewers. The values were so high as to be inconsistent with the rest 
of the data and were omitted from subsequent analysis. The remaining data 
still showed considerable scatter, but were utilized to aid in deriving 
required relationships. 
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TABLE 24. SURFACE BOD LOADINGS FOR RESIDENTIAL AREAS AS DERIVED FROM EFFLUENT MEASUREMENTS 



City 



Site or 
Station^ 



Sewer 
System 



Catchmant 
Area 
ac (ha) 



Tulsa 



Bucyrus 



Atlanta 



Roanoke 



Milwaukee 



3 

f 

■T 

8 

9 

11 

ia 

IS 



17 



Confed. Ave. 
Blvd. 

McDan St. 
Harlan 
Caspian 
Fed. Pris. 

Trout Run 
Murray Run 
24 St. 

Hawley Rd. 



Separate 
Separate 
Separate 
Separate 
Separate 
Separate 
Separate 
Separate 

Combined 
Combined 
Combined 

Combined 

Combined 

Combined 

Separate 

Separate 

Separate 

Separate 
Separate 
Separate 



550 
507 
197 
211 

64 
815 
212 

74 



1,1?9 

2,421 

968 

954 

517 

1,498 



(223) 
(205) 
(80) 
(85) 
(26) 
(330) 
(86) 
(30) 



179 (72) 

614 (249) 
378 (153) 



(457) 
(980) 
(392) 
(386) 
(209) 
(606) 



997 (404) 

909 (368) 

1,034 (419) 



Annual 
Precip. 
in. (cm) 



Runoff 
Coef. 



BOD 
Cone. 

mg/L 



Surface Loading 
Ib/ac-day 

(kg/ha-day) 



Population 

Density 
Persons/ac 

(persons/ha) 



48 (122) 



35 (89) 



48 (122) 



34 (86) 



0.39 
0.41 
0.35 

0.31 
0.42 
0.42 
0.33 
0.56 
0.31 



120 
107 
108 

210 
84 

286 

7 

20 

26 



0.0381 
0.0901 
0.0417 
0.0899 
0.0544 
0.0963 
0.0679 
0.0688 

1.017 
0.953 
0.821 

1.94 

1 05 

3.58 

0.069 

0.334 

0.240 



(0.0428) 
(0.1012) 
(0.0468) 
(0.1009) 
(0.0611) 
(0.1081) 
(0.0762) 
(0.0772) 

(1.142) 
(1.070) 
(0.922) 

(2.178) 
(1.179) 
(4.019) 
(0.077) 
(0.375) 
(0.269) 



7.13 

8.93 

11.55 

11.37 

13.67 

9.57 

2.36 

11.22 

11.7 
9.1 
5.0 

10.9 
16.6 
13.2 

9.7 

7.3 

4.8 



0.0363 (0.0408) 11.0 
0.0428 (0.0481) 6.6 
0.0233 (0.0262) 9.7 



(17.61) 
(22.06) 
(28.53) 
(28.08) 
(33.76) 
(23.64) 
(5.83) 
(27.71) 

(28.9) 
(22.5) 
(12.4) 

(26.9) 
(41.0) 
(32.6) 
(24.0) 
(18.0) 
(11.9) 

(27.2) 
(16.3) 
(24.0) 



Combined 495 (200) 31 (79) 0.40 49 0.377 (0.423) 35.0 (86.5) 



Source 



32 



33 



34 



35 



36 



Note: Surface loadings are taken directly from the source if given, or derived from mass emission data. Otherwise, equation 29 is used (for cities for which runoff coefficient 
and BOD concentration are listed). 



TABLE 24. (CONT'D) 



' 
























Population 






• 






Catchment 


Annual 




BOD 


Surface Loading 


Density 








Site or 


Sewer 


Area 


Precip. 


Runoff 


Cone. 


Ib/ac-day 


Person s/ac 






City 


Station* 


System 


ac 


(ha) 


in. 


(em) 


Coef. 


mg/L 


(kg/ha-day) 


(persons /ha) 


Source 




Wash., D.C. 


Good Hope Run 


Separate 


265 


(107) 


41 


(104) 






0.063 


(0.0711 


37.6 


(92.3) 


37 






B4 


Combined 


105 


(43i 










0.247 


(0.277) 


43.6 


(107.7) 








G4 


Combined 


222 


(901 










0.381 


(0.428) 


52.6 


(129.9) 






Des (VIoines 


S-1 


Separate 


315 


(128) 


31 


(79) 


0.10 


48 


0.093 


(0.104) 


7.4 


(18.3) 


38 






m 


Separate 


356 


(144) 






0.10 


63 


0.121 


(0.136) 


5.3 


(13.1) 








0-3 


Combined 


4,050 


(1,640) 






0.15 


69 


0.199 


(0.223) 


7.5 


(18.5) 








0-6 


Combined 


5,600 


(2,267) 






0.15 


95 


0.275 


{0.309 


8.3 


(20.5) 








0-8 


Combined 


1,350 


(547) 






0.15 


68 


0.197 


{0.221) 


10.9 


(26.9) 








0-8A 


Combined 


927 


1375) 






0.15 


77 


0.222 


(0.249) 


10.9 


(26.9) 






Cincinnati 


Mt. Washington 


Separate 


27 


(11) 


40 


(102) 






0.0904 


(0.1015) 


9.0 


(22.2) 


39 


ro 


Durham 


E-1 


Separate 


56 


(23) 


45 


(114) 


0.29^ 


25 


0.202 


(0.227) 


14.9 


(36.8) 


40 






W-1 


Separate 


169 


(68) 






0.35'^ 


61 


0.596 


(0.669) 


2.6 


{6.4) 








W-2A 


Separate 


69 


(28) 






0.34^ 


38 


0.361 


(0.405) 


11.0 


(27.2) 








W-2B 


Separate 


138 


(56) 






0.36'^ 


51 


0.513 


(0.576) 


13.4 


(33.1) 








N-1 


Separate 


183 


(74) 






0.36'^ 


71 


0.714 


(0.802) 


4.2 


(10.4) 






Seattle 


Low Dens. 
Med. Dens. 
High Dens. 


Separate 
Separate 
Separate 


c 

c 
c 




36 


(91) 






0.04 
0.07 
0.13 


(0.045) 
(0.079) 
(0.146) 


11.0^ 
22.0*^ 
30.0^ 


(27.2) 
(54.3) 
(74.1) 


41 




Windsor 


Labadie Rd. 


Separate 


30 


{12) 


33 


(84) 






0.059 


(0.066) 


20.0 


(49.4) 


42 



Silc or statiein as listed in source documentation. 
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FIGURE 11. RESIDENTIAL BOD LOADINGS vs POPULATION DENSITY 
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The data of Droste [42] from Windsor were about the only data for 
the study region that provided both a BOD loading estimate and population 
density for a residential area. An earlier study in Windsor [43] also 
provided these required data but was conducted on a developing area that 
included construction activities and was considered somewhat atypical by 
Droste [42], Hence, it was not included in Table 24. 

Other studies of importance to the analysis included work in 
Halifax that had been published [44, 45] as well as work in Toronto, 
Burlington, and Aurora, sponsored by Environment Canada that had not yet 
been completed. Of these, only the latter Halifax study of Bhatia [45] 
provided BOD loads of the type required in this analysis, but data were 
taken there for only three months in 1971, and it was somewhat 
questionable as to whether they were representative of a whole year. 
(Almost all the data included in Table 24 were taken over a period of 
several months.) In addition, population densities were not given. 
However, the average surface runoff BOD load calculated for the 2.18 acre 
(0.9 ha) Cambridge St. residential area of 0.038 Ib-BOD/acre-day (0.042 
kg-BOD/ha-day) is within the range presented in Table 24 and Figure 11. 
Clearly a synthesis of data from current studies in Ontario and elsewhere 
across Canada will provide better estimates of parameters needed for this 
type of analysis. 

4.5 Pollutant Load Prediction for Ontario Assessment 

4.5.1 Form of equation 

Surface pollutant loads generated by the pollutant load 
estimating equation are assumed to "wash off" on an annual basis for 
purposes of the assessment. Hence, they must be representative of actual 
measured effluent loads. Moreover, they should be functionally related to 
causative factors in a reasonable manner. They are expected to be 
functions of land use and population density. In addition, there are 
apparent observed geographical variations in, say, dust and dirt loadings, 
although it is not immediately obvious as to why these loadings should 
differ in a commerical or industrial area from one point in the country to 
another, other than on the basis of climate. 



74 



The key climatic parameter is precipitation, since the more preci- 
pitation that occurs, the more likely it is that pollutants will be washed 
off the surface and appear as effluent loads instead of being removed by 
other means such as street sweeping or wind. Total annual pollutant loads 
from storm runoff are lower in arid regions for this reason [47], 
Precipitation includes both rain and snow, on an annual basis, incorporat- 
ing the assumption that pollutants accumulate during periods of snow cover 
and eventually are washed off during periods of melt. 

These considerations led to the selection of a prediction equa- 
tion, in which the loading is proportional to precipitation, for all land 
uses. It is also proportional to a function of population density for 
residential areas which is intended to account for many other implicit 
factors such as age of area, imperviousness, runoff coefficient, etc., all 
of which are functions of population density. This formulation was 
easily applied because precipitation and population density data were 
readily available. However, these parameters were about the only ones 
{other than areas) that were available, ruling out more complicated 
functions. The loading, M, is thus represented functionally as: 

M = a • filP) f2{PD) W) 

where the coefficient a and functions f-j and f2 are determined 
below. The procedure followed develops appropriate parameters for 
residential areas first, which are then extended to other land uses. 

4.5.2 Precipitation function 

If average BOD loadings for the cities in Table 24 {omitting 
Atlanta. Bucyrus, and Durham) are plotted versus annual precipitation (not 
shown), no clear relationship is indicated. Hence, the data are simply 
averaged to obtain the factor a and fi{P) of equation (30) for BOD. 
That is, it is assumed that the loadings are directly proportional to 
precipitation, such that zero precipitation generates zero storm water 
pollution. This is supported by equation (29). Hence, 

fl{P) = P (311 

and the parameter a is obtained as an average of the seven{ remaining 
cities in Table 24 for which separate data are available. Thus, for BOD 
for residential areas, 



I 7 loading- -.u, nnn 

a=4.2 — : ^=0.00219 '^'^°° 



7 1=1 P^- ' ac-day-in. 



= 0.799 



Ib-BOD 



. = 0.353 .^^3^^^ (32) 



ac-yr-in ' ha-yr-cm 

Annual average BOD loadings for residential areas are now predicted by 

M = 0.799 . P . f2(PDd) (33) 

where: M = annual average BOD loading for separate sewered, 
residential areas, Ib-BOD/ac-yr; 
P = annual precipitation, inches; and 
PDj = developed population density, persons/acre. 

For combined areas the equation is identical, except that a 
parameter b is employed instead of a to distinguish between combined and 
separate areas. For BOD for residential areas, the value of 6 is computed 
using average values for Des Moines, Milwaukee, and Washington, D.C. from 
Table 24. 

,-. 1 I '^^^ - 0.00902 l^-^O^ 



3 1=1 P. ' ac-day-in 

= 3.29 -iMOD_. 1.46 kMOD_ (34) 

ac-yr-in ha-yr-cm 

Annual BOD loadings for residential areas served by combined sewers are 
thus, 

M = 3.29 • P • fgtPDd* *^^* 

where parameters are as previously defined. 

It may be seen that, for the same population density and 
precipitation, combined BOD loadings are 3.29/0.799 = 4.12 times higher 
than separate loadings. This agrees with an independent survey of 
available data by Lager and Smith [46] in which average BOD concentrations 
in combined sewage of 115 mg/L are 3.83 times greater than the average BOD 
concentration of 30 mg/L in separate sewers. The difference in loadings 
is due mainly to residual matter left In conduits between storms, since 
simple mixing of storm water and dry-weather flow, or differences in 
population density between separate and combined sewer areas, will not 
explain the four-fold variation in concentrations and loadings. 
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4.5.3 Population function 

The data in Table 24 and Figure 11 incorporate all the avail- 
able information about the relationship of BOD loadings with population 
density implied by equations (29) and (30). In order to extend the data 
base slightly further, it is assumed that combined area loadings increase 
with population density, PD, in the same manner as do separate area 
loadings. The data base can then be extended slightly by normalizing by 
the average loading for separate and combined areas. Omitting the data 
from Atlanta, Bucyrus. and Durham, Table 25, was prepared. Finally, the 
data of Table 25 was plotted, as shown in Figure 12. A point has been 
added that represents the loading in open space of 0,00982 Ib-BOD/ac-day 
(0.0110 kg-BOD/ha-day) where presumably the population density is zero. 
(The derivation of this value is shown later.) 

Inspection of Figure 12 shows such scatter that no statistically 
significant relationship is likely to be derived from the data. Rather, 
an argument must be made upon the expected form of the functional 
relationship, and the data used only to obtain a calibration. This 
relationship is expected to be similar to those developed earlier for 
imperviousness and curb length, namely increasing rapidly at low 
population densities and leveling off at high ones. 

The concentration of storm water pollutants is M/AR, or 

a P fJPD ,) 
^/^^ = K[0.15'q.751 ] P (36) 

where: I = 0.096 Pd/^'"^ " ''''' ''ho ^°d, or 

I = 0.096 PD^^-^'* (37) 

and K is a conversion factor, for example, the value that appears in 
equation (29). Depression storage is omitted in the approximation of 
annual runoff. Thus, 

. [f PD ] 

M/AR . — ^—5 (^1 

K[0. 15+0.072 PD^O.54] 



m 



1 



TABLE 25. NORMALIZED BOD LOADING DATA 



Average 

Loading 

Ib-BOD/ 

ac-day 

(kg-BOD/ 

ha- day) 



City 



Loading/ 
Ave. Loading 



Population Density 
Persons/ ac Persons/ ha 



Separate 


0.0693 


Tulsa 


0.550 


7.13 


17.61 


Areas 


(0.0778) 




1.300 


8.93 


22.06 








0.0602 


11.55 


28.53 








1 .297 


11.37 


28.08 








0.785 


13.67 


33.76 








1.390 


9.57 


23.64 








0.980 


2.36 


5.83 








0.993 


11.22 


27.71 






Roanoke 


0.524 


11.0 


27.2 








0.617 


6.6 


16.3 








0.336 


9.7 


24.0 






Wash., D.C. 


0.909 


37.6 


92.9 






Des Moines 


1.342 


7.4 


18.3 








1.746 


SvB 


13.1 






Cincinnati 


1.305 


9.0 


22.2 






Seattle 


0.577 


11.0 


27.2 








1.010 


22.0 


54.3 








1.876 


30.0 


74.1 






Windsor 


0.851 


20.0 


49,4 


Combined 


0.271 


Wash.. D.C. 


0.911 


43.6 


107.7 


Areas 


(0.304) 




1.405 


52.6 


129.9 






Milwaukee 


1.391 


35.0 


86.5 






Des Moines 


0.734 


7.5 


18.5 








1.014 


8.3 


20.5 








0.727 


10.9 


26.9 








0.819 


10.9 


26.9 



Note: Values obtained from Table 24, omitting data from Atlanta, 
Bucyrus and Durham, 
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FIGURE 12. NORMALIZED BOD LOADINGS VS POPULATION DENSITY 
Data are from Table 25 
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It is assumed that f«(PDd) is 

f2(PDd) = a + b PDd'" (39) 

where: a = 0.142 = value at PDd = 0, 

and developed population density will be used for consistency. Note that, 

depending on the assumed value of m, the concentration of storm water 

pollution will vary accordingly. Since no firm arguments can be made on 

the nature of the concentration function, it will be assumed that m is 

equal to the approximate exponent in the runoff equation or m = 0.54. 

Thus, f2(PDj) = 0.142 + b PD^^'^*. Lastly, all data points with 

a PDd ranging from 5 to 15 persons per acre (12 to 37 persons per ha) 

are averaged to obtain a calibrated value of f2{PDd) = 0.895 at 

10 persons per acre (25 persons per ha). This range is chosen because 

data from most cities fall within it. Thus, the final equation is 

f^lPD^) = 0.142 + 0.218 PD^j (40) 

where: PD, = developed population density, persons per acre. 

The reasonableness of equation (40) can be checked by estimating the 

variation in concentration as a function of population density. From 

equations (33) and (40), the annual BOD loading is: 

54 
M = 0.799 • P • (0.142 + 0.218 PD^ ) (41) 

and annual runoff, AR, using the approximate New Jersey [26] equation for 
imperviousness is: 

54 
AR = [0.15 + 0.75 (0.096) PD^. ] • P (42) 

Thus, Q .. 
0.113 + 0.174 PD >^'^^ 
M/AR = ^ (43) 

K[0.15 + 0.072 PD^^*^'^] 

Using K = 0.227 for these units, this ratio, which is plotted on Figure 13, 
shows concentration increasing with population density, which does seem to 
be reasonable. The range of average annual concentrations is lower than 
values shown in Table 24 since it represents the average over the total 
residential area of a city. Unquestionably, the data base for estimating 
pollutant loads is ^ery weak, and the resulting estimating equation, 
supported by such a weak foundation, should be used with extreme caution. 
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FIGURE 13. ANNUAL AVERAGE RESIDENTIAL BOD CONCENTRATION 
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4.5.4 Conversion for alternate land uses and pollutants 

Different pollutants and land uses will generate different 
loadings for at least three reasons. First, the dust and dirt loadings 
for different land uses differ. Second, the conversion factor of curb 
length per area is different for different land uses. Third, the 
pollutant fractions (as a fraction of dust and dirt) are different for 
different land uses. These factors are used to extend the equation 
developed for BOD for residential areas to similar equations for 
conmercial . industrial and open land uses and for suspended solids, 
volatile solids, total PO^, and total N. 

It is assumed that fractions and ratios of pollutants as they 
appear in effluent will be the same as those determined from analysis of 
surface accumulation data. The parameters shown in Table 26 are used for 
conversion purposes. They are selected from the extensive survey material 
prepared by APWA [3], Where no data are available for pollutants as a 
fraction of surface dust and dirt, use is made (as a second choice) of 
similar data developed for pollutants as a fraction of total solids (TS). 

The BOD data are first converted to other land uses, using 
equation (23) as indicated below: 

a (i. BOD) = a (res, BOD) • -r~ • ^^^^ ' p^ -"" t^^J 

%es \,res res. BOD 

where: dd. = dust and dirt accumulation on land use i, 

lb /day-curb mile; 

G, . = curb miles per acre of land use i from Table 23, and 

F.'nnn = fraction of dust and dirt that is BOD on land use i. 

1 ,dUU 

For example, the parameter a for BOD for commercial land use for separate 
areas is: 

/ Dnn^ n 7qq 166 x 7190 0.069 

a (com, BOD) = 0.799 -3537^:65" • "0:^7? 

- 2 59 JMOD, . 1 14 k9-B0D ^45j 

^■^ ac-yr-in ' ha-yr-cm 

where the number of 353,465 is the average product of dd • F for BOD is 
equal to: 



Dust and Dirt (DD) 
Accumulation 



TABLE 26. SURFACE LOADING AND POLLUTANT FRACTION DATA 

Residential 
Single Multi 
Family Family Average Commercial Industrial 



Open 



tb 



62 



day-curb mite 

kg 



113 87.5 



166 



day-curb meter 
BOD - ppm of DD 

P|wn of total solids (TS)^ 
Total PO4 - ppm of DD 
ppm of TS 
Total N - ppm of DD 
ppm of TS^ 
Suspended Solids - ppm of TS^ 
Volatile Solids - ppm of TS^ 



17 32 
5,260 3,370 



24.8 



47 
7,190 
29,640 83,800 



609,200 582,300 
353,000 367,700 



319 



90 

2,920 

25,850 



14.2 

18,990 

1,670 



10,170*^ 

619,500 453,200 
306,100 437,500 



AH 
Date 



50 159 



45 



170 



664* 



Volatile Solids - ppm of TS" 353,000 367,j . . 

Except as noted, all data are from Table 82 [3], Missing entries are not 
given in original table or not used in analysis. 

^Values taken from Table 19 [3]. 
°Sum of K-N plus NO3-N. 
^Value for organic N only. 



62 X 5.260 -*- 113 X 3,370 
2 

After determination of BOD for each land use, i, other quality parameters, 
j, are computed on the basis of relative values of the fractions, F. Thus, 



46) 



47) 



a (i,j) = a{i,BOD) J '^ 

l,BOD 

For example, the parameter a for total PO^ (TPO^) in commercial 

areas is: 

170 ^^-™4 
a (com, TPO4) = 2.59 X j^^ = 0.0612 ^^^TFTTTi 

kg-TPO- 

= 0.027C 7-— — ~ 
ha-yr-cm 

For total nitrogen, N, in residential areas the calculation is similar but 
includes the average product of dd • F, 

(rp<; H\ - Q-799 - 664 ■ (62 + 113)/2 . „ ,^^ Ib-N (40) 

" *'^"»^' 353,465 - ^-^^^ dL-y r -in ^^^^ 



= 0.058 



kg-N 
ha-yr-cm 



83 



For open land use, and for suspended solids and volatile solids, no data 
are available for fractions of dust and dirt, so fractions of total solids 
are used for values of F in the ratios. For example, for suspended solids 
in commercial areas, 

, (CO..SS) = 2.59 X ^11^ = 13.0 ^I^ (49) 



= 7.95 



kq-SS 
ha-yr-cm 



Computations for combined areas are carried out in the same manner to 
calculate the b parameters. 

Results from the U.S. assessment show that there is a point after 
which the magnitude of street sweeping frequency has no effect on the 
computed values of average annual pollutant concentrations [2]. In Oes 
Moines, Iowa, if the streets are swept less frequently than every 20 days, 
then the STORM model, which accounts for street sweeping, does not show 
any significant reduction of pollutant load. For intervals up to 20 days, 
a linear buildup may be assumed. Thus, the final estimating equation 
includes a street sweeping factor v as a function of the sweeping interval, 
Ng, in days, i.e. 

f/20 if interval of street sweeping, f, is such that 

N /20 if < N < 20 days 



Y = 



50 



1.0 if N^ > 20 days 



No variation due to type of sewer system is included. For this assessment 
the street sweeping intervals exceeded 20 days so it was unnecessary to take 
explicit account of this factor. The final result is shown in Table 27. 

Use of the same adjustment factors for combined and separate areas 
leads to the same ratio e/a = 4.12 for all entries in the table. On the 
basis of measured concentration data [46], the assumption appears valid 
except for solids, wherein some studies have shown higher ratios of volatile 
solids to suspended solids, for example, in combined sewage than in storm 
runoff alone [44]. 

The BOD loadings are compared to dry-weather flow loadings in Table 
28 for residential land use. Stonn and combined runoff can be seen to be 
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TABLE 27. POLLUTANT LOADING FACTORS FOR ONTARIO ASSESSMENT 



The following equations say be used to predict annual average 
loading rates as a function of land use. precipitation and population 
density. 

Separate Areaa: M^ - a(l,J) • P • f2^*'°d* ' ^ acre-yr 
Corri?ined Areaer H^ - B(1,J) • P • ^2^^^^^ ' ^ ,cre-yr 

where M - pounds of pollutant J generated per acre of 

land use 1 per year, 

P - annual precipitation, inchea per year, 

PD - developed population density, persons per sere, 

a,B - factors given In table below, 

y ' Street sweeping effectiveness factor, and 

f-(PD.) - population density function. 
2 a 

Land Uaes : 1-1 Residential 

i - 2 Commercial 

1-3 Industrial 

1 ■■ A Other Developed, e.g., parks, cemeteries, schools 

(assume PD. - 0) 
d 



Pollutants: J - I BOD , Total 

- 2 Suspended Solids (SS) 

- 3 Volatile Solids, Total (VS) 

- 4 Total PO, (as PC.) 

- 5 Total N 

0.54 

Population Function: 1 - 1 ^j^^^a^ " ^-^^^ + 0.218 • PD 

i - 2,3 rA?V>) - 1.0 

I . 4 fU^Oj - 0.U2 

Z a 

Factore a and 8 for Equaii na : Separate factors, a, and combined factors, 
e. have units Ib/acre-in. To convert to kg/ha-cm, multiply 
by 0.4ii2. 









PC 


llutant , 


J. 








1 . 


Land Use, 1 


1. BOD^ 
0.799 


2. SS 


3. VS 


k. PO, 
0.03 36 


5. N 




Resident lal 


16.3 


9.4 S 


0.131 


Separate 
Areas, a 


?. 


Commerc ial 


2.59 


18.0 


11.4 


0.06 1 2 


0.239 


^ 


Industrial 


0.994 


23.8 


11.8 


0,0579 


0.226 


u. 


Other 


0.0969 


2.31 


2.23 


0,00852 


0,05 19 




1 


Residential 


3.29 


67.2 


38.9 


0.139 


0.540 




? 


Cotanerclal 


10.7 


74.2 


47.0 


0.252 


0.98S 


Areas, 6 


3 
U 


Industrial 
Other 


4.10 
0.399 


98.1 
9.52 


48.6 
9.19 


0,239 
0.0351 


0.931 
0.214 



Strpet Swppping: Fa 
N , (days): 



rtor r Is a function of Street sweeping interval, 



fN /20 If < N < 
1.0 If N * 20 d 



20 days 
days 



TABLE 28. COMPARISON OF BOD LOADINGS 



BOD Loading 
Ib/ac-yr kg/ha-yr 



Separate Areas 21 24 

Combined areas gg 9§ 

Dry-Weather Flow* M^ 697 

DWF at 85% Treatment* 93 



Assume residential land use; PD^^ = 10 persons/acre (24.7 persons/ 
ha) and P = 30 in/yr (76 cm/yr), and no influence of street sweeping 
(y = 1 ). 

* Assuming 0.17 Ib-BOD/persons-day (0.08 kg-BOD/persons-day) 

comparable to treatment plant effluent, although on a city-wide basis they 
would be greater because of higher loadings for commercial and industrial 
areas. Of course. BOD loads in both storm and combined sewage are in 
addition to the dry-weather flow loads since the usual BOD load for the 
latter of 0.17 Ib/person-day (0.08 kg/ person- day) is based upon measure- 
ment of flows actually received at treatment plants. The data from which the 
loadings shown in Table 24 were derived reflect discharges over and above 
those received by the plants. 

4.6 Tabulation of Ontario BOD Loads and Runoffs 

In order to minimize the volume of material presented for each 

city in the assessment, only BOD. PO^, and N loadings were tabulated. 

The equations indicated in Table 27 may easily be used to calculate 

loadings of any of the desired parameters, given the precipitation and 

population density of the area of interest. As described in Section 4, 

land use variations are determined by first computing the fraction of 

undeveloped land in the urbanized area. The remaining land has a constant 

distribution of land uses, and can be used to weight the pollutant 

loadings factors to give an average over-all land use as follows: 

- 4 

M = P z W. • a. • f2 (PD) • Y (51) 

i=l ' 1 

The land use distribution fractions for cities other than the nine test 
cities, w . , are given below. 
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Land Use Fraction, w. 
1 — , 

Residential 0,525 

Commercial 0.103 

Industrial O.UO 

Open Q.232 

1.000 

When equation (51) is applied to BOD loadings for separate areas, the 
result is: 

M = 0.419 - P . (0.142 + 0.218 PD^^'^^) + 0.409 P (52 

where: M = average annual BOD loadings over four land uses. 
Ib-BOD/ac-yr; 
P = annual precipitation, inches; and 
PD^ = developed population density, persons/acre. 

For application to combined areas, the result is: 



M = 1.726 . P . (0.142 + 0.218 PD^^"^^) + 1.685 P 



;53 



These composite equations may easily be applied over the non-test 
cities. Note that equation (53) is simply equation (52) multiplied by 
4.12. Using equation (51) and Table 27, similar equations can be 
developed for total phosphate and total nitrogen. These parameters are 
important as far as pollution of the Great Lakes is concerned. Thus, for 
T-PO^ in separate areas, 

M = 0.0176 • P • (0.142 + 0.218 PD^°'^^) + 0.0147 P (54) 

and for T-N in separate areas, 

M = 0.0688 • P . (0.142 + 0.218 PD^°'^^) + 0.0580 P (55) 

Equations (54) and (55) should be multiplied by 4,12 for combined area 
T-PO^ and T-N loads, respectively. 

Dry-weather flow loadings are computed simply on the basis of 
population density assuming average annual BOD generation of 0.17 
Ib/person-day (0.08 kg/ person- day). Thus, 

Mp = 62.1 . PD^ (561 



where: M^ = average annual dry-weather flow BOD loading, 
Ib-BOD/ac-yr. 

Dry-weather loadings of total phosphate and total nitrogen may be found 
using data of Lager and Smith [46], which indicate that the average 
concentration of T-PO. (as PO^) and T-N in dry-weather flow are 15 
and 20 percent, respectively, of the BOD concentration. Thus, equation 
(56) may be multiplied by these percentages for calculation of T-PO^ 
and T-N loadings. 

Results of the analysis may be seen for each city in Tables 29, 
30, 31, 32, 33 and 34. Area weighted averages for all areas are also 
shown. Owing to relatively low precipitation and relatively high 
population densities, dry-weather pollutant loads are generally higher 
than corresponding wet-weather values. However, as seen previously in 
Table 28, if 85 percent treatment is assumed for dry-weather loads, the 
resulting values are comparable to storm and combined sewered loads. 

4.7 List of Variables 

i Value of function f2(PD(i* ^^®" ^°d ^ ° 

A. Area of land use i 

M Wet-weather runoff, depth/time 

^ Normalized loading factor for separate sewered areas, 

mass/area- time- length 

b Coefficient in function fo(PD ,) 

BOD Biochemical oxygen demand 

BODc Biochemical oxygen demand at five days 

b 

B Normalized loading factor for combined sewered areas, 

mass/area- time-length 
C Concentration, mass pollutant/ total mass 

COD Chemical oxygen demand 

com Abbreviation for commercial 

CM Runoff coefficient 

DB Dust and dirt 

dd. Dust and dirt loading factor for land use i, 

mass/time-curb- length 
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DWF ADbreviation for dry-weather flow and dry-weather flow runoff, 

depth/time 
f,(P) Factor for adjustment of pollutant loads, a function of 

precipitation 
f (PD ) Factor for adjustment of pollutant loads, a function of 

population density 
F . Fraction of dust and dirt or land use L that is pollutant p 

G|^ ^ Length of curb per area of land use L 

T Street sweeping factor 

1 Imperviousness as a fraction or percent 

Ig Initial abstraction (loss) from precipitation, depth 

K Factor to convert runoff times concentration to pollutant 

loadings 

m Exponent in function f^(PD,) 

I a 

M Pollutant loading, mass/area-time 

M Pollutant loading averaged over different land uses, 

mass/area-time 

^ Pollutant loading in combined sewered areas, mass/area-time 

M^ Pollutant loading under dry-weather conditions, mass/area-time 

M^ Pollutant loading in separate sewered areas, mass/area-time 

n Number of storms per year 

N Nitrogen 

N Number of dry days preceding a storm 

N Street sweeping interval, days 

P Precipitation rate, depth/time 

pQ Mass of pollutant on surface at end of previous stom, mass 

^i j ^^^^ °^ pollutants on surface of land use L at beginning of 

storm, mass 

P^ Precipitation depth during one stonn, length 

PD Population density, persons/acre 

PD^ Developed population density, persons/acre 

PO^ Phosphate or total phosphate 

r Correlation coefficient 

res Abbreviation for residential 

p Water density, mass/volume 



SS Suspended solids 

TOC Total organic carbon 

TPO. Total phosphate 

TS Total solids 

VS Volatile solids 

m. Fraction of total area that is land use i 
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TABLE 29. DRY WEATHER BOD LOADINGS 



N(]| UPBANIZFD AREA 



AWNL. I DRY hEaThEK HOn 
PRtCPI TUBb/ACRE-YFA») 
IN/VRI CGMBjSTaRMlUNSFWI AVER 



1 I 




\ 

; a 

3 

4 
5 

lb 
17 

i? 

23 

2^ 
55 

2b 

27 
26 
29 

30 
31 

3^ 
33 
34 
35 
3b 
37 
■^6 
39 

ao 

43 
^4 
as 

4b 

47 
ak 

49 

55 
•^b 



AJAX 
AURH 
IbARR 

IBFLL 
I6RAM 

ibraw 

iburl 

(CHAT 
ICMINi 

icoan 

I DUNH 
lETOB 

Iueor 

GIIEL 

HAMT 

IKTNn 

iKirr 

ILFAM 
tLTNO 
I Lni-jn 
I f-l A R K 
IMTDI 

\m 

INIAP, 
INORT 
I OAKV 
lOOIL 
I OSHA 
I OWE '-J 
IPFT"^ 
IPICK 
IPDRT 
I PORT 
IPREfl 

1ST. 
SCAP 

isiir 

I3TRA 

'IgT 

TREM 
I NaCl 
I WELL 
'..HIT 

i rtilUF) 

! YHRK 
I YORK 
I YHRK 



RA 

IE 

EVILLE 
PTQK 

TFORD 
IMGTQn 
MAM 

GUArOMSY 
U^G ■ 

!cOKE 

GETOw^J 
pw 

LTQN 
STQW 

H.-wATEPLna 

SAY 
QH 

AM 

A^JQ 

ISS4UGA 

ARKET 

RA FALLS 

H BAY 

ILLF_ 

LIA 

WA 

S U U N n 
KBGRDIitiH 
ERING 

COLBnuRNF 

ERTE 
THN 

MONO HILL 
CATHARI^'ES 
THQMAS 
NA 

S T F. , 1 A H T t 
BORHuGH 
OE 
TFORD 

DER BAY 

NTO 

TON 

ACEBURG 

AND 

HY 

SOR 
itA^T 

,'iaRTM 



^3 

31 
32 

34 
31 
32 
32 
30 
^1 
34 
32 
31 
33 
32 
7 t 

32 
35 
3tt 
30 

}^ 

3 



^1 

33 

35 
31 
3fa 
34 
35 
32 
31 
34 
34 
34 
31 
32 
35 

\1 
]i 

38 

n 

30 
34 
34 

3U 

^1 



lfi40 



lUlfl 




n 










n 

990 

I39fe 



2098 



1794 

1061 

985 





?*79 

1992 

n 





1?6« 







e;?72 





927 

1322 

1779 

n 

1^23 

1958 





1370 

21/7 

n 

1637 

I'^bO 


ISdb 
f] 

IfiOft 

171'S 




73« 
742 
977 

^51 
bfiO 
66\ 
7dA 
497 
81? 
960 
I05n 
M04 
1033 
10S4 
928 
378 
550 
80b 
720 
9ifl 

8g« 

39b 
76'^ 

ao** 
37« 
749 
83fl 

1 004 
bet 
554 

1^73 

96W 
944 

709 
9b4 

«b4 
471 

40? 

543 

h58 
729 

104b 
813 
56S 
888 

109b 

^5? 
972 
?59 

i 

H4 



19P 

199 

274 
35b 
355 
259 

354 
t?a 

199 
197 
31B 
197 

198 

218 





299 

197 

m 

208 



1 

273 
327 
273 
19ft 
35^ 
204 
198 
200 
201 
19(1 
199 
209 
355 
198 

25? 

197 
197 
2b5 
317 


198 
198 
l9fl 
19Q 

?7fc 


ft 





518 
419 
498 
579 
629 
614 

5h 

hO?. 
790 
500 
541 
724 
529 
564 
543 
843 
778 
612 
550 
4H0 
610 
4§4 
475 
598 
434 

573 

762 
561 

513 
613 
468 
556 
494 

4A3 

576 

502 
451 
b07 
570 
884 
^15 
777 

532 

538 

495 

6PO 

1855 
573 
507 

4^6 

b52 

5lb 

lbl7 

1283 
834 



wFlGHTEn AVE. 



i?.75l 1545, I Tife 



255 



743, 
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TABLE 3U. WET -LEATHER 30D LOADINGS 



•^QIURBANIZED A^tL 



I ANNL. I 
iPRtCPt 
! IN/YH I 
I ' 






T hEA 
dS/Ac 



THEH 
IUN3 



EA"?) 
EWI AVER 



1? 



I 1 



I 



4 
5 
b 
7 
6 
IS 

23 



1^ 

26 
29 

32 
73 

34 
35 



38 

39 

ao 
«l 

43 
u4 
45 

51 



l^{ 



5^ 
55 

56 



I AJAX 
I AHRHrt 
IBar^RI 
IbELl £ 
IBRawP 
IBRant 

iburli 

ICHATH 

ICHINC 

IC^BOU 

IDUND* 

I&TOBI 

IGaLT 

IGFQRG 

tGUELr 

ihaRtl 

NGS 
TfH 
lUEAHl 

UTNOS 
LDNOO 

Rarka 

MIDI A 

NFWMA 
NlAr,3 
NriSTH 
OiKVl 
QRILL 
QSHAW 
gWEN 
i rrTtn 

m' 

pnRT 

PREST 

RICHM 
ST. C 
3T, T 
SARIN 
8LT. 

§caAb 

STRAT 

SMOnu 

THUND 
TDRON 
TRENT 
WALLA 
HfLLA 
WHITB 

YORK 
YORK* 

YORK, 



A 

e 

VILLE 

Ton 

FORD 
hJ G T K' 
AM 

uacousy 
s 

CQKF 

F r n ;j M 

H 

TjN 

TON 

.-waterlho 

NgTUN 

M 

ND 

SSANGA 

RkpT 

A 

SOUND 

SuRuuG" 



COLI 



BrnjRNF 

ONO HILL 

ATHARINES 

HOMAS 

gjE. MARIE 

E 

FORD 

On 

CERURG 

kin 
"U 

Y 

OR 
TOCK 



iUSH 



BAY 



EAST 

^DRTN 



131 
31! 
^21 
^41 
31 

li\ 
301 
311 
341 
32 1 

32l 
^3 

32 
35 
34 

iu 

}^ 



I 1 



l^\ .2 

0,0 
0.0 
135. 
0. 
0,0 
0,0 
106, b 
0.0 


0,0 
0,0 
0,0 
0.0 
15. fa 
43.9 
0,0 



15 



0.0 



1^4, 1 

0.0 

132.3 

152. b 

0.0 

0.0 

1^5.5 

iTa.3 

0,0 

\lt:\ 

0.0 
115.8 

0.0 
132.0 
129,9 

0,0 



29.7 

2a;5 



52; 



23.8 

22,6 

24.7 

?4,9 
22.3 
48,9 
58,2 
56.7 
30.3 
43,6 
37.2 
26,1 
23.6 
45.6 

142.5 
24,7 
35.3 
^3.2 
24,6 
36,6 
25.0 

114.7 
90.6 



WEIGHTED AVE, 132.751137.0 



25.71 21,3 



44,5 



92 



TABLE 31. DRY WEATHER T-PO. LOADINGS 







M M - L . 
KRLCP 


D R T w f 


"a^ 


rHEr^ 


pna 


'kjd 


UOB^NlZep ARLi 


Ct 


ds/acre-yfarj 






IM/YK 


C J'-ri 


STuRp-| 


UK'SEW 


AVER 


1 


A jAy 


53 


27fa.l 


lie. 


7 


^"^,7 


77.7 


2 


AMHnrtA 


31 


0,0 


in, 


3 


3a; 


9 


62p6 


3 


B AkRIF 


,?cf 


0.0 


l"b 


b 


29, 


9 


7a. 8 


a 


i5ELl£VlLLF 


3a 


212. fl 


'^'d 


7 


^i 


k 


86.9 


5 


d^AMPTG^j 


31 


0,0 


l<^d 





9a, a 


6 


dRA'^JTPORQ 


32 


O."; 


1^2 


2 


53 




-92.1 


7 


C HATH Ad 


32 


0,0 


lU 


5 


ti3 


3 


fl6.5 


8 


30 


ia9.fa 


7i4 


b 


53 


\l 


9Q-a 


9 


CHl'iG-'ACU'iaV 


31 


0.0 


I '"-l 


,9 


^9 


118.5 


1 10 


cnfiou^Ci 


34 


0,0 


ia7 


:i 


29 


,9 


75,0 


H 


DUNOAS 


32 


0,0 


157 


29 


• b 


81.1 


1 p 


ETtlolCUKE 


3 1 


0.0 


120 


.7 


^7 


,6 


ine.7 


13 


SALT 


3i 


0.0 


155 


,0 


29 


.6 


79.5 


U 


tiEURGETOwf' 


32 


0.0 


1*^2 


7 


?9 


.7 

,7 


Aa.b 


15 


CiiFi PN 


33 


0.0 


159 


.3 


32 


«l.a 


1^ 


HAMTLTON 


^ 


iS§:5 


ft2 


.5 







.0 
.0 


26.6 

ib.e 




KTTrH.-WATERLnD 


3a 


0.0 


U 


.0 


aa 


.a 




LPAMIN'GTQM 


30 


ii a, 9 


.0 


?9 


*6 


5^ 


LTi-^nSAY 

L u N I ) u '^ 
HARKAf'i 


33 

37 


0,0 


lag 


.8 


?S 


.2 

,0 


72,0 
66.1 


31 


159.3 


^^^ 


,3 


31 


,3 


23 


MIDLAND 

MtssissauriA 


36 


la7,9 


.5 


30 


,a 


B9.7 


?u 


31 


0,u 


117 


.P 


as 


• ^ 


25 


NfTi^HARKFT 


31 


0,0 


121 


.3 


32 


,7 


65.2 


26 


NTARrtA FALLS 


?3 


131.9 
29d,9 


5b 


.7 


ai 


»o 


Ab.l 


27 


NHf^TH fw. ¥ 


3!3 


112 


.a 


a9 


,0 


117,3 


?a 


OAKVILLF 


1^ 


O.u 


125 
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TABLE 32. UET WEATHLR T-PO,, LOADINGS 



^^qiurbanized area 



ANNL. I WET HEATHEH PO*! 
PRECPI CLH3/ACRE-YEAR) 
IN/YNI CQMblSTQRHlLlNSFWI AVER 



I 51 



lAJAX 
I AUROR* 
IBARRIE 
I 



m^'. 



LE 



I 

9 
JO 

n 

3 

a 



I 



ibrantfqrq 

iBu«LINCtON 
HATHAM 
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isiilRGETOWN 

iguelph 
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HA . 
KlTcH.-WATERLnU 



11 
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TABLE 33. DRY WEATHER T-N LOADINGS 
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TABLE 34. WET WEATHER T-N LOADINGS 
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3. OVERALL COST ASSESSMENT 

This section develops and applies a methodology to estimate the 
cost of controlling pollution from storm water discharges province wide. 
Costs of controlling combined sewer overflows, storm water runoff, and/or 
providiig tertiary treatment are compared. A general methodology for 
determining wet-weather pollution control costs is presented. Then, a 
procedure is described for determining the relationship between storage, 
treatment, and pollutant control for control of wet-weather flows. 
Generalized predictive equations are developed based on relatively 
intensive studies of four cities: Burlington, Kingston, St. Catharines, 
and Sault Ste. Marie. Knowing this "production function" one can 
determine the optimum combination of storage and treatment, by combining 
this information with data on the cost and performance of the available 
control options. This information is combined to produce the Ontario 
assessment. Results are presented for all cities >^ 10,000 persons. 
Related reports describing this methodology are available [62, 63], 

5.1 Methodology 

5,1 ,1 Principles 

There are several economic theories which, when applied to 
environmental resources management, assist in the decision-making process. 
One such theory is production theory, which provides techniques that aid 
in evaluating items such as the optimum size of a reservoir for water 
supply and flood control, or a wastewater treatment plant for pollution 
control. When the cost of inputs such as the reservoir or treatment plant 
is known, the cost of achieving a desired level of output (e.g., water 
Supply Or" puMuLiufi control) may be detennined. 

In storm water management, the inputs are usually in the form of 
a storage capacity and a treatment rate. Storage is expressed in terms of 
mil lion galluMi or inches over a certain area, typically the watershed 
being analyzed. The unit for treatment is either million gallons per day 
or inches per hour, using the same area as storage. 

When the degree of wet-weather control is considered as a single 
output, it can be expressed either in terms of the percent of the runoff 
treated or the number of overflows per year. This is with respect to 



quantity only and is, therefore, dependent upon the input storage capacity 
and treatment rate. 

When dealing with only two inputs it is feasible to use a 
graphical method to find the optimum combinations. Isoquants can be 
constructed which represent equal levels of output for different 
combinations of input (see Figure 14). For example, each isoquant could 
represent a specific percent of the runoff treated for different 
combinations of storage and treatment. Isoquants have the following 
properties [48]: 

1) Two isoquants cannot intersect. Intersecting isoquants would 
imply two different levels of output from the same input. 

2) Isoquants slope downward and to the right because as one 
input increases it takes less of the other input to achieve 
the same level of output 

3) Isoquants are convex to the origin because of the decreasing 
ability of one input to be substituted for another to obtain 
a given level of output. This is known as the principle of 
diminishing marginal rate of substitution. 

Also on Figure 14, a series of parallel lines have been 
constructed. These lines represent combinations of input 1 and input 2 
which may be achieved at the same total cost. The lines are known as 
isocost lines. The slope of the isocost lines is the relative unit cost 
between input 1 and input 2. The most economical combination of input 1 
and input 2 to produce a desired level of output is the point where the 
isocost lines become tangent to the isoquant representing the desired 
level of output. 

The line which joins the points of tangency among several 
isoquants and the isocost lines is called the expansion path. After the 
expansion path has boen detennined, the optimum combination of inputs car? 
be determined for any level of output by finding the intersection of th^ 
isoquant representing the desired level of output and the expansion path. 

The maximum output for a c^iven cost may be found by constructing 
the isocost line for the given total expenditure. The slope of the 
isocost line is the relative unit cost of the two inputs. The intercept 
of the axis depicting input 1 would be the allowed total cost divided by 



Cj = UNIT COST OF INPUT I 



--ISOQUANT C2 = UNIT COST OF INPUT 2 



EXPANSION PATH 



rj 



13 

a. 




S0C0S7 LINE 



INPUT I 



FIGURE 14. DETERMINATIONS OF LEAST COST COMBINATION OF INPUTS. 
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the unit cost of input 1, From this information, the isocost line may be 
drawn. The point where the isocost line intersects the expansion path 
gives the combination of inputs which produces the maximum output at the 
given cost. 

The storm water quality management problem can be expressed in 
the more compact mathematical form shown below: 



minimize 



Z = c^{S) + c^(T) (57) 

subject to 

f(R^;S.T) = 

R^;S,T 2 ^ 

where: Z = total control costs, 

c (S) = storage costs, 

c (T) = treatment costs, 

S = storage volume, 

T = treatment rate, 

R = percent pollutant control, and 

f(R -S.T) = production function relating the level of pollutant 

control attainable with specified availabilities of 

storage (S) and treatment (T). 

The next three sections describe: 

- the available storage/treatment options - their costs and 
effectiveness; 

- the production functions for evaluating tradeoffs between 
storage and treatment; and 

- the solution to the optimization problem yielding the optimum 
expansion path for any city. 

Given this information, the final assessment is presented. 

5.1.2 Control technology and associated costs 

A wide variety of control alternatives is available for improv- 
ing the quality of wet-weather flows [49, 50, 51], Rooftop and parking 
lot storage, surface and underground tanks and storage in treatment units 



are the flow attenuation control alternatives. Wet-weather quality 
control alternatives can be subdivided into two categories: primary 
devices and secondary devices. Primary devices take advantage of physical 
processes such as screening, settling, and flotation. Secondary devices 
take advantage of biological processes and physical-chemical processes. 
These control devices are suitable for treating storm water runoff as well 
as combined sewer overflows. However, the contact stabilization process 
is feasible only if the existing sewage treatment plant is of an activated 
sludge type. The quantities of wet-weather flows that can be treated by 
this process are limited by the amount of excess activated sludge 
available from the dry-weather plant. 

At the present time, there are several installations throughout 
the U.S. designed to evaluate the effectiveness of various primary and 
secondary devices. A summary of the design criteria and performance of 
these devices is presented in Table 35. Based on these data, the 
representative performance of primary devicp*; i^ assumed to be 40 percent 
BOD^ removal efficiency, and that of secondary devices to be 85 
percent BOD^^ removal efficiency. 

"Storage" devices will typically be used in conjunction with the 
above "treatment" devices. The two purposes are interrelated. Wastewater 
detained a sufficient time in a storage unit will undergo treatment. On 
the other hand, treatment units also function as storage units in that 
they equalize fluctuations in influent flow and concentration. DiToro 
presents approaches for evaluating the equalization and treatment which 
occurs in both of these units [52], The STORM model, which was used in 
this assessment, assumes the configuration tor storage and treatment shown 
in Figure 15. No treatment is assumed to occur in storage, and "treat- 
ment" is assumed to be complete removal of all pollutants routed through 
treatment. Thus, for the purpose of this assessment, no treatment is 
assumed to occur in storage and control costs are assigned accordingly. 
This assumption tends to underestimate the costs of storage since all 
provisions for solias handling are Included in treatment. 

5.2 Cost of treatment and storage 

Cost data for installea -wet-weather treatment devices Are 
listed in Table 35. Since wet-weather control facilities operate 
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TABLE 35. WET-WEATHER TREATMENT PLANT PERFORMANCE DATA 



Device 



Reported 
Design Criteria BOD5 Removal 

Control Alternatives gpm/ft^ (l/min-m^| Efficiency, r} 



Primary Swirl Concentrator^''^ 60.0 

Microstrainer= 20.0 

Dissolved Air 
Flotation w/ Chemical 

Addition*^ 2.5 

Sedimentation* 0.5 

Representative Performance 

Secondary Contact Stabilization* Cont 



Physical-Chemical^ 
Representative Performance 



(2,448.0) 
( 816.0} 



0.25-0.50 
0.40 - 0.60 



( 102.0) 


0.50 - 0.60 


{ 20.4) 


0.25 - 0.40 




0.40 


0.25 hours 


0.75 - 0.88 


3.00 hours 




3.00 hours 


0.85 - 0.95 




0.85 



""Field and Moffa [53] 

^APWA [54] 

■^Maher [55] 

^Lager and Smith [46] 



Performance data based on domestic 
wastewater treatment. 

Agnew et a1 [56] 

^Estimate based on performance of these 
units for domestic wastewater. 



TABLE 36. INSTALLED COSTS FOR WET-WEATHER TREATMENT DEVICES 

Annual Cost per mgd:'^ (m^/day) S/yr 
Operation and 





Capacity 


Amortized Capital^-^ 


Maintc 


jnance 


Total 




Control Device 


mgd 


(m^/day) 


per mgd 


(m^/day) 


per mgd 


im^/day) 


mgd 


(lO^m^/day) 


Swirl Concentrator^ 


8,9 


(34,500) 


5,600 


{ 1.48} 


2,100 


(0.55) 


7,700 




(29.6) 


Microstrainer^ 


7.4 


(28,700) 


14,230 


( 3.76) 


3,895 


(1.03) 


18,125 




( 8.6) 


Dissolved Air Flotation^ 


25.0 


(96,900) 


71,706 


(18.94) 


16,700 


(4.41) 


88,406 




(33.5} 


Contact Stabilization^ 


20.0 


(77,500) 


120,000 


(31.70) 


24,000 


(6.34) 


144,000 




(54.4) 



Based on 8% interest for 20 years. 

Construction cost; does not include sludge handling costs. 
^Field and Moffa [53] 
^Maher [55] 

Lager and Smith [46] for Racine, Wisconsin adjusted to ENR = 2,200. 

Operation and maintenance costs based on 480 hours of operation (9 50.341/ 
1000 gallons ($0.0126 per 1000 litres). 

^Agnew et al [56]. Operation and maintenance costs based on 960 hours of 
operation. 

All gallons are U.S. All costs are U.S. dollars. 
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Runoff 




Receiving Water 



FIGURE 15. STORM MODEL SIMULATION OF STORAGE AND 

TREATMENT FOR WET-WEATHER QUALITY CONTROL 

intermittently, annual operations and maintenance costs are greatly affected 
by the number of hours the facility is utilized. As a general rule, a 
facility will operate a greater amount of the time if it incorporates 
storage. An examination of Table 36 reveals that annual operation and 
maintenance costs are 16.7 percent of the total annual costs for the 
contact stabilization unit. In the case of the swirl concentrator, the 
percentage is 27.3. Annual operation and maintenance costs for other 
units fall in between these two values. Based on this analysis, it was 
decided to assume annual operation and maintenance costs as 20 percent of 
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the total annual costs for all treatment devices. Cost functions 
developed for various wet -weather quality control devices are presented in 
Table 37. These costs include provisions for sludge handling, engineer- 
ing, contingencies, and land costs. 

All treatment units exhibit economies of scale, i.e., unit cost 
decreases as plant size increases. Thus, there is an incentive Lu build 
larger units. The optimum size treatment unit can be found by comparing 
the savings in treatment cost of going to a larger unit with the increased 
piping costs. For example, if one is comparing building two 10 myu 
(37,850 m^/day) plants with building one 20 mgd {75.700 m"^/day} 
plant and a pipeline, the breakeven pipe length, L is found using 

Two plants One plant + pipeline 

s{10)^ + s(10)^ = s(20)^ + K{10)^{L) (58 

where: s, x, z, K and y = coefficients. 

Unfortunately, data on the number and flow rates of storm water discharges 
in urban areas could not be found. Thus, it is not possible to determine 
the optimum mix of treatment plants and pipelines. Therefore, representa- 
tive treatment costs were used as shown in Table 36. A 30 mgd plant size 
was selected since it represents a reasonable upper limit on the range of 
strong economies of scale. The average costs are based on a microstrainer 
for primary treatment and contact stabilization for secondary treatment. 

Cost data on detention basins built in the Chicago area for 
temporary storage of runoff are listed in Table 38. Costs of storage 
tanks built for the purpose of wet-weather quantity and quality control as 
well as for dry-weather quantity control are also included in this table. 
Due to the wide variations in these figures, an attempt was made to verify 
these costs using excavation costs as the basis. Storage costs based on 
unit excavation costs are listed in Table 38. The unit cost of equaliza- 
tion and the estimated costs of rooftop and parking lot storage basins for 
sewage treatment plants are also shown in Table 38. Lastly, analysis of 
recent estimates of storage costs developed by Benjes et al indicate the 
following unamortized capital cost C ($ x 10 } as a function of 
storage volume, S (million gallons) [58], 
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TABLE 37. ANNUAL COST FUNCTIONS FOR WET WEATHER CONTROL DEVICES ^'^'^" 



Device 



Primary 









Annual Cost 


$/yr 






Amortized Capital 


Operation 


and 


Total 




CA = ITm 




Maintenance 


TC = wT' 




orISm 




OM = pT'' 


or wS' 


Control Alternative 


1 


m 


P 


q 


w z 


Swirl Concentrator'^''*'^ 


1,971.0 


0.70 


584.0 


0.70 


2.555.0 0.70 


Microstrainer^'* 


7,343.8 


0.76 


1,836.0 


0.76 


9,179.8 0.76 


Dissolved Air Flotation® 


8.161.4 


0.81 


2,036.7 


0.84 


10,198.1 0.84 


Sedimentation" 


32,634.7 


0.70 


8,157.8 


0.70 


40,792.5 0.70 



Secondary 



Representative Primary Device — Total Annual Cost = $4,000 per mgd ($1.06/m'/day) 



Contact Stabilization^ 


19,585.7 


0.85 


4,894.7 


0.85 


24,480.4 


0.85 


Physical-Chemical^ 


32,634.7 


0.85 


8,157.8 


0.85 


40,792.5 


0.85 



Representative Secondary Device Total Annual Cost = $15,000 per mgd ($3.93/m^/day) 



Storage 



High Density {15/ac) 


51,000.0 


1.00 


Low/ Density (5/ac) 


10,200.0 


1.00 


Parking Lot^ 


10,200.0 


1.00 


Rooftop*^ 


5,100.0 


1.00 



51,000.0 


1.00 


10,200.0 


1.00 


10,200.0 


1.00 


5,100.0 


1.00 



Representative Annual Storage Cost' ($ per ac-in} - $122 e°-^^'^'^d' 



T = Wet-weather treatment rate in million gallons per day (mgd) 
S = Storage volume in million gallons 
one million gallons = 3,785 m3 

ENR = 2,200. Includes land costs, chlorination, sludge handling, 
engineering and contingencies. 

Sludge handling costs based on data from Batelle Northwest [57] 
^Field and Moffa [53] 
^Benjes et al [58] 

Lager and Smith [46] 
^Maher [55] 
^Agnew et al [56] 
"^Wiswall and Robbins [59] 
Vor T < 100 mgd {378,500 m-^/day). 

PO^ = developed population density, persons/acre. 
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TABLE 38. INITIAL CAPITAL COST OF STORAGE FACILITIES' 



Capacity 
mil ga' (1,000 m^t 



Capital Cost 
S/gal (S/liter) 



Stuiage heservoirs" 

Hillside Park 
Heritage Park 
Oak Lawn 

Middle Fork North Branch 
Wilke-Kirchoff 
Melvina Dutch 
Oak Hill Park 
Dolphine Park 
Average 

Storage Tanks® 
Cottage Farm, Boston*^ 
Spring Creek, New York'= 
Chippewa Falls, Wisconsin'^ 
Humboldt Avenue, Milwaukee^ 
Seattle, Washington 
Whittier Narrow, Columbus*^ 
Average 

Based on Excavation Costs^ 

$2/yd^ {S2.62/m^) 
$5/Yd^ (S6.54/m^) 

Equalization Basins for Dry 
Weather Sewage Treatment 
Plants^ 



Other*' 

Parking Lots 
Rooftops 



11.4 

36.5 
7.8 
195.5 
32.6 
53.8 
25.1 
53.8 
52.1 



1.3 

10.0 
2.8 
4.0 

32.0 
4.0 



9.0 



( 43.1) 
(138.0) 
( 29.5) 
(740.0) 
(123.0) 
(204.0) 
( 95.0) 
(204.0) 
(197.0) 



: 4.9) 
: 37.8) 
; 10.6) 
[ 15.1) 
1121.0) 
( 15.1) 



34. i; 



w. 


( 3.8) 


3.0 


(11.4) 


10.0 


(37.8) 


IM 


( 3.8) 


3.0 


(11.4) 


10.0 


(37.8) 



(0.003} 
(0.003) 
(0.005) 
(0.003) 
(0.008) 
(0.008) 
(0.005) 

(0-0Q3) 

0.019 (0.005) 



0.01 
0,01 
0.02 
0.02 
0.03 
0.03 
0.02 
0.01 



5.21'' 


(1.38) 


2.33 


(0.62) 


0.29 


(0.08) 


0.55 


(0.14) 


0.25 


(0.07) 


1.70 


(0.45) 


1.72 


(0.45) 


0.01 


(0.003) 


0.025 


(0.007) 



0.22 


(0.06) 


0.10 


(0.03) 


0.06 


(0.02) 


0.39 


(0.10) 


0.28 


(0.07) 


0.25 


(0.07) 


0.10 


(0.03) 


0.05 


(0.02) 



Earthen Basin 
Earthen Basin 

Earthen Basin 
Earthen Basin 
Earthen Basin 
Earthen Basin 
Earthen Basin 
Earthen Basin 



Covered Cone. Tanks 
Covered Cone. Tanks 
Asphalt Paved Basin 
Covered Cone. Tanks 
In line 
Open Concrete Tanks 



Earthen Basin 
Earthen Basin in Rock 



Earthen Basin 
Earthen Basin 
Earthen Basin 
Concrete Basin 
Concrete Basin 
ConcTete Basin 



^Based on ENR = 2,200 

^Metropolitan Sanitary District of Greater Chicaqo 

^Also used for storm water treatment 

^Includes pumping station, chlorination and outfall facilities. 

\ager and Smith [46] 

^Soil Conservation Service, Gainesville, FL. 

^Flow Equalization - Plus for Wastewater Treatment Plants, Civil Eng. . 9/75. 



h 



'..'is-^Qll dna Robbins [59] 



Unit Cost @ S = 10 mg 
Type Equation $/gal {$/litre) 

Earthen C = 0.02b s'^'^^ $0,013 ($0.0034) 

Concrete w/o Cover C = 0.350 S°'^^ $0,133 ($0.0350) 

Concrete w Cover C = 0.400 S^*^^ $0,250 ($0.0660) 



The data indicate wide variation in the costs of storage. Thus, the 
relatively simple relationship shown in Table 38 was used. Annual storage 
costs are estimated as a function of gross population density. At low 
population densities, land values are relatively low. Thus land-intensive 
storage facilities, e.g., shallow ponds, can be used. At higher popula- 
tion densities, land values increase to the point where storage tanks 
become more economical. The data presented in Table 38 are based on 
differing assumptions regarding land values. In some cases the land is 
free (part of an easement) whereas in others it is valued highly. Thus, a 
simplified approximation was used. The curve was derived using an 
unamortized capital cost of $0.10 per gallon ($0,026 per litre) for 
PD. = 5 persons per acre (12.4 per ha) and $0.50 per gallon ($0,132 
per litre) for PD = 15 persons per acre (37.1 per ha). 

5.3 Kelationsmp petween btorage/Treatment and Percent Pollution 

Control 

5.3.1 Use of STORM 

STORM was used to evaluate various storage/treatment options 
for controlling storm water runoff pollution. This model assumes that the 
study area can be characterized as a single catchiTient from which hourly 
runoff is directed to storage and treatrnent. 

STORM uses a simplified rainfal 1 /runoff relationship, neglects 
the transport of water through thp city and assume^; a s/e^ry simple 
relationship between storage and treatment. However, these simplifica- 
tions are essential if one hopes to do a continuous simulation. The 
continuous simulation approach was used because o.o general concurrence 
exists regarding an appropriate single event that one should analyze. 
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The degree of control can be expresseo in t2nns of the percent of the 
runoff treated, the annual number of overflows, or the amount of 
pollutants discharged to the receiving waters. 

As described in the User's Manual, STORM computes the runoff 
based on the composite runoff coefficient and the effective precipitation 
[23], The depression storage must be satisfied before the runoff 
coefficient is applied to the precipitation. The amount of depression 
storage available in ditches, depressions, and other surfaces is a 
function of the past precipitation and the evaporation rates. Each hour 
that runoff occurs, the model compares it to the treatment rate. As long 
as the runoff rate is less than or equal to the treatment rate, all the 
runoff passes directly through the treatment plant and storage is not 
utilized. When the runoff rate exceeds the treatment rate, the excess 
runoff is sent to storage. If excess runoff occurs frequently enough to 
exceed the storage capacity then overflow occurs. When runoff falls below 
the treatment rate then storage is depleted at the excess treatment rate. 
The hourly occurrence of treated runoff, stored runoff, and runoff that 
has overflowed is taoulated for the entire record of rainfall. Included 
in the output is the annual number of overflow events and the percentage 
of the runoff that overflowed to the receiving waters. This type of 
analysis was carried out for different storage capacities and treatment 
rates. 

5.3.2 STORM input data 

STORM requires several input parameters that characterize the 
urban area under study. These include hourly precipitation, total area, 
land use types and percentages, percent imperviousness, and curb length 
per area for each land use [23]. In order to apply STORM, these data were 
collected or derived for nine urban areas within the Ontario portion of 
the Great Lakes Basin, The local data were collected by on-site 
interviews. 

Precipitation data were acquired from the Department of the 
Environment. The records were of varying length, from a few years to 
several decades. Additionally, the records were supplied in two parts. 
One part included only rainfall and was usually restricted to the months 
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from April to October. This portion was previously formatted on magnetic 
tape, which facilitated a frequency analysis of this restricted rainfall 
record. This allowed the selection of one year (limited) of rainfall to 
characterize the entire record of each city. After the year was selected, 
the accompanying months of precipitation (primarily snowfall) were taken 
from written records and added to the restricted rainfall record. The 
precipitation for the selected year for each city was totalled and 
compared to the mean annual average. These values are shown in Figure 16. 
In some cases, it was apparent that although the rainfall frequency for 
the selected year was typical, the snowfall was not. Therefore, four of 
the nine cities, Burlington, St, Catharines, Kingston, and Sault Ste. 
Marie, were chosen to make the final STORM runs. It should be noted that 
precipitation patterns over Ontario show two distinct trends. In Figure 
17, it can be seen that below 44°N, the monthly precipitation is very 
stable. Above 44°N the trend is toward a peak in the summer months. 
Therefore the selection of the four cities to represent the entire region 
seemed adequate. Sault Ste. Marie provides a good representation of the 
northern section of Ontario, while the three other cities adequately 
represent the southern section of the study area. 

Snowfall is a significant component of Ontario's precipitation 
total, ranging from a low of 4.0 inches (10.2 cm) water equivalent near 
Windsor, to a high of 14.0 inches (36 cm) north of Sault Ste. Marie. 
STORM uses the degree-day for the computation of snowmelt with the option 
of applying this formulation or simply ignoring snowmelt and allowing snow 
to act as rainfall. A runoff frequency analysis was performed on the 
cities of Thunder Bay and Windsor to determine the effect of this routine. 
Runs were made with and without the method shown above. The results are 
shown in Figures 18 and 19. The figures show a difference in frequency; 
therefore, the snowmelt routine in STORM was used in the final STORM run. 

Dally evaporation rates for each month were estimated from a 
report by Phillips and McCulloch [60], These were available for only a 
few stations in Ontario and the evaporation rates for the four test cities 
were assigned on the basis of proximity to these stations. The depression 
storage is assumed to be 0.01 inches for all cities. The input data used 
to run STORM for the developed areas of the four selected test cities is 
summarized in Table 39. 
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FIGURE 16. MEAN ANIMUy^L PRECIPITATION FOR THE GREAT LAKES BASIN. 

{After Phillips and McCutloch,^^ Chart 21.) 




FrCURE 17. MEAN MONTHLY PRECIPITATION (After Brown. McKay and Chapman,^! p. 35) 
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TABLE 39. STORM INPUT DATA FOR TEST CITIES 



Study Area: 
Developed Area: 
Test Year: 
Precipitation: 



Burlington 

8687 ac (3516 ha) 

1973 

32.38 in. (82.25 cm) 



Daily evaporation rates for each month, Jan-Dec, in-/day (cm/day) 
0.00 0.00 0.05 0.05 0.13 0.15 0.17 0.14 0.10 0.06 0.05 0.00 
(0.00) (0.00) (0.13) (0.13) (0.33) (0.38) (0.43) (0.36) (0.25) (0.15) (0.13) (0.00) 



Study Area: 
Developed Area: 
Test Year: 
Precipitation: 



Kingston 

4706 ac (1905 ha) 

1965 

37.81 in. (96.04 cm) 



Daily evaporation rates for each month, Jan-Dec, in/day (cm/day) 
0.00 0.00 0.05 0.05 0.13 0.15 0.17 0.14 0.10 0.06 0.05 0.00 
(0.00) (0.00) (0.13) (0.13) (0.33) (0.38) (0.43) (0.36) (0.25) (0.15) (0.13) (0.00) 



Study Area: 
Developed Area: 
Test Year: 
Precipitation: 



St. Catharines 

10976 ac (4442 ha) 

1973 

32.37 in. (82.22 cm) 



Daily evaporation rates for each month, Jan-Dec, in/day (cm/day) 
0.00 0.00 0.05 0.05 0.13 0.15 0.17 0.14 0.10 0.06 0.05 0.00 
(0.00) (0.00X0. 13) (0.13) (0.33) (0.38) (0.43) (0.36) (0.25) (0.15) (0.13) (0.00) 



Study Area: 
Developed Area: 
Test Year: 
Precipitation: 



Sault Ste. Marie 

8516 ac (3446 ha) 

1969 

36.69 in. (93.19 cm) 



Daily evaporation rates for each month, Jan-Dec, in-Zday (cm/day) 
0.00 0.00 0.04 0.04 0.09 0.11 0.12 0.10 0.06 0.04 0.04 0.00 
(0.00) (0.00) (0.10) (0.10) (0.23) (0.2^) (0.30) (0.25) (0.15) (0.10) (0.10) (0.00) 



Burlington, Kingston, St. Catharines - estimated from Guelph "pan" evaporation 
land evaporation = 0.70 "pan" evaporation. 

Sault Ste. Marie - estimated from Seney, Michigan "pan" evaporation: land 
evaporation = 0.70 "pan" evaporation. 
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5.4 Results 

For each storage/treatment rate combination there is a value 
for the percent of the runoff and pollutants which are "treated". 
Preliminary analysis of STORM runs made for the U.S. assessment indicated 
little year to year variation in results. Thus, only one year of 
precipitation was used to derive the isoquant curves. By making several 
runs at different combinations of treatment and storage, points were 
generated representing different levels of control. Then isoquants were 
drawn connecting the points that represent combinations of storage 
capacities and treatment rates which give equivalent percent runoff and/or 
pollutant "treated". If the concentration of pollutants is constant and 
"treatment" efficiency, n, is 1.0, then percent runoff control is 
synonymous with percent pollutant control. Obviously, this is not the 
case. Thus, account needs to be taken of: 

1) treatment efficiency, and 

2) variable concentration due to first flush effects. 

5.4.1 Adjustment for treatment efficiency 

Let R denote the percent runoff control and equal treatment 
plant efficiency. If R, denotes the percent pollutant control, then 
to realize R, , one needs to process R,/n of the runoff. Note that 
R^ may be percent BOD removal, percent SS removal, etc. In Table 39, 
representative treatment efficiencies, in terms of BOD removal, were 
derived for primary and secondary devices. These values are as follows: 

Assumed Efficiency, n 
Treatment Device (BODc Removal ) 



Primary 0.40 

Secondary 0,85 

Thus, if one desires 25 percent BOD^ removal with a primary device, 
then 62.5 percent of the runoff volume must be processed, whereas only 
29,4 percent of the runoff needs to be processed if a secondary device is 
selected. Thus, to convert percent runoff control isoquants to percent 
pollutant control isoquants, one uses: 

R = R^/n (59) 
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5.4.2 Adjustment for first flush 

STORM estimates the percent pollutant control a^ well as 
percent runoff control. The STORM model runs incorporated the standard 
first flusti assumption which is used in the model, i.e., that the amount 
of pollutant removal at any titi^.e, t, is proportional tc the amount 
remaining and that a unifonn rainfall of one-half inch per hour would wash 
away 90 percent of the pollutant in one hour. If a first flush is 
assumed, then storage and treatment can be operated more effectively 
because of the greater relative importance of capturing the initial 
runoff. The first flush is accounted for by defining the output in terms 
of pollutant control directly. 

5.4.3 Mathematical representation of isoquants 

The storage/treatment isoquants are of the form: 

T =T^ + (T^ - T^)e"^^ (60) 

where: T = wet-weather treatment rate. Inches per hour, 

T = treatment rate at which isoquant becomes asymptotic to 

the ordinate, inches per hour, 
T„ = treatment rate at which isoquant Intersects the 

abscissa, Inches per hour, 
S = storage volume, inches, and 
K = constant, inch' . 

A relatively large storage reservoir is required to operate the 
treatment unit continuously. Thus, first flush effects would be dampened 
out and the effluent concentration from the reservoir should be relatively 
uniform. Thus, if storm water entering the treatment plant was a 
relatively uniform concentration, then T-. can be found as follows: 

^1 - 87760 'ilo' = *« '^^' 

where: a = coefficient, 

AR = annual runoff, inches per year, and 
R = percent runoff control. 
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By relating the parameters ^-i » ^o - T, and K to the level of 
control R, one equation was developed for each of the four cities, 
The ly - T and K terms versus R were found to be of the 
following general form: 

. cR 



T T 

'2 '1 



(62) 

K = de'^^ (63) 

Based on this analysis the following general equation for the isoquants is 



obtained: 



T = aR + be 



cR-(de"^^)S 



(64) 



The values of paraineters a, b, c, u, and f for various cities are 
presented in Table 40, for percent pollutant control. The correlation 
coefficients for the equations for the four cities are also shown in this 



+-aKl 



Tn 



■i-ho *;■; + 



The results for the four cities are shown in Figures 20, 21, 22 
and 23. Each figure shows the isoquants calculated by the isoquant 
equation. Alsc shown are some actual data points for a treatment rate of 
U.Ol inches per hour and varying amounts of storage, 

TABLE 40. VALUES OF PARAMETERS AND CORRELATION 
COEFFICIENTS FOR ISOQUANT FACTORS 

Percent BOD Control with First Flush, n = 1.0 



Test City 



a 


b 


c 


in hr ' (% R) ' 


in hr ' 


{% K) ' 


(cm hf' ) 


(cm hr"^ ) 





Burlington 0.0000121 0.0017093 0.0414918 

(0.0000310) (0.0043400} 

Kingston 0.0000127 0.0013611 0.0384055 

(0.0000320) (a0034G00) 

St. Catharines 0.0000130 0.0016126 0.0434050 

(0.0000330) (0.0041000; 

Sault Stc. Marie 0.0000166 0.0018704 0.0449201 

(0.0000421) (0.0074900) 



d 

{cm-* ) 



m^ 



f Correlation Coefficients 

[%Rr* T^-Ti- K=de^" 

be= 



,cR 



210.4327 0.0298024 0.994 -0.93R 
(32.8000) 

241.9431 0.0306992 0.995 -0.992 

(95.3000) 

240.4267 0.0298348 0.995 -0.983 

(94.7000) 

191.309 0.0334145 0.992 -0.994 

(75.300) 
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FIGURE 20. STORAGE-TREATMENT ISOQUANTS FOR PERCENT 

BOD REMOVAL WITH FIRST FLUSH FOR BURLINGTON 
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FIGURE 21. STORAGE-TREATMENT ISOQUANTS FOR PERCENT BOD 
REMOVAL WITH FIRST FLUSH FOR KINGSTON 
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FIGURE 22. STORAGE-TREATMENT ISOQUANTS FOR PERCENT BOD 
REMOVAL WITH FIRST FLUSH FOR ST. CATHARINES 
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FIGURE 23. STORAGE-TREATMENT ISOQUANTS FOR PERCENT BOD 
REMOVAL WITH FIRST FLUSH FOR SAULTSTE. MARIE 



The optimum expansion path can be found using 

where: c^ = unit cost of storage, 

c = unit cost of treatment, and 

MRS = marginal rate of substitution of storage for treatment. 

The value of c^. and c^ are presented in Table 34. 

Analysis of the figures indicates that if c^/c^ < 25, 
then treatment alone should be used. From Table 34: 



C3 ,,,^0A6(P0) 

J 4. * 4. - CO ^-..^ , acre-inch 

For primary treatment, c = $2,610/ — ^^^ — . 

Thus, even at zero population density, c-j-Zc = 21.4 so that the 
optimum policy is to use treatment only. 

So that, 

for secondary treatment, letting c /c^ = 25, and knowing that 

c^ = S9.800/i^IS^, yields 

PD . = 7.29 persons per acre. 

If PD , is higher than about 7.5, then the relative cost of storage is 
such that it is again optimal to use treatment only. Using 7.5 persons 
per acre as the cutoff, then 12 of the 56 cities would use treatment only 
for the secondary control level. The remaining 44 cities would select a 
mix of storage and treatment. 

It is simple to find the optimum expansion path graphically for 
the four test cities. Unfortunately, these results need to be extrapolated 
to all urbanized areas. It appeared that an analytical approach would 
provide a more general and consistent procedure. Thus, the isoquant 
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parameters were adjusted based on the runoff in the city under 
consideration relative to the reference city, i.e., let 

AR. = annual runoff in city i; i = 1,2, ..., 56 

AR . = annual runoff for test year in test city for region j 
(see Figures 20-23); j = 7.17,38,41. 

Then, the isoquant coefficients are 

a. . = AR./(8,76 X 10^) (66) 

AR 
b._. =^bj. (67) 

=i.j =^j' (68) 

AR. 

where a. ., b. ., c. ., d. ., and f. . are parameters for city i 

in region j; and b., c., and f. are the parameters for the test city in 

region j. The test cities are denoted as follows: 

C City 

7 Burlington 

t7 Kingston 

Si St. Catharines 

ft Sault Ste. Marie 

5.4,4 Wet-weather quality control optimization 

The wet-weather optimization problem, assuming linear costs, 
may be stated as follows: 

minimize 

Z = CgS + c^T (71) 

subject to 

T = T^ + (T^ - T^)e"'^^ 
T,S > 
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Solving this constrained optimization problem yields: 

S* = max[^ In^LKd^ - T^)]. 0] <72) 

where: S* = optimal amount of storage, inches, 

and T* = T^ + (T2 - T^)e''^^* ^73) 

where: T* = optimal amount of treatment, inches per hour. 

Note that T* is expressed as a function of S*. so it is necessary to find 
S* first. Knowing S* and T*, the optimal solution is: 

Z* = c S + c^T* (7*1 

where: Z * = total annual cost for optimal solution, dollars per acre. 

Data needed to estimate T , T and K have already been presented 
in the previous section. 

For a primary device, 

c^ = S4,000/mgd = $2,510/^^^^|^ ($1 ,05/mVday) 

and ^ = 0.40. 

For a secondary device, 

c = $15,000/mgd = $9,800/ ^'^'^^^^" {$2.32/m /day). 

For storage cost, 

C3 ($/acre-in.) - 122 e°-^^<'"^' <") 

where: PD = gross population density in persons per acre. 

The above optimization procedure was programmed to generate 
curves, e.g.. Figure 24, showing percent pollutant removed versus total 
annual costs for primary and secondary treatment in conjunction with 
storage. Note that for wet-weather control, marginal costs are increasing 
because of the disproportionately larger sized control units needed to 
capture the less frequent larger runoff volumes. 
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FIGURE 24. CONTROL COSTS FOR PRIMARY AND SECONDARY 
UNITS IN STORM SEWERED AREAS IN BURLINGTON 



These results also permit one to decide whether a primary or a 
secondary level is more cost-effective in controlling smaller percentages 
of pollution. As seen in Figure 24, a primary control device is less 
expensive for low removals (say < 20 percent), but it loses effectiveness 
at higher levels because of the disproportionately large storage require- 
ments. Costs will be reported for 25, 50, and 75 percent control levels. 
Thus, the secondary cost curve can be used in this range. The primary 
curve will not be discussed further. 

The curves shown in Figure 24 can be approximated by functions of 
the form: 

Z = ke ^'^l (76) 

where: Z = total annual cost, dollars per acre, 
k B = parameters, 

R, = percent pollutant removal, < R, < R", , and 
R, = maximum percent pollutant removal. 

The resulting costs for 25, 50, and 75 percent pollutant control for 
combined, storm, and unsewered areas are shown- in Tables 41, 42 and 43, 
respectively. Note that the reference city and values of the cost 
equation parameters are also shown, 

5.4.5 Estimating number of overflow events 

Some urban areas have used the number of overflow events per 
year as an indication of level of control due to different storage/treat- 
ment combinations. The objective in this case would be to find the most 
economical combination of storage and treatment which would not allow the 
annual number of overflows to exceed a predetermined value. It would not 
seem logical to increase the treatment rate or storage capacity if the 
number of overflows did not decrease. 

The number of overflow events is affected by the definition of an 
"event" used in the STORM model wherein an event is defined as starting 
when storage is utilized and ending when storage is depleted. Even though 
overflow may take place in two separate time frames, the two occurrences 
are considered to be parts of the same event if storage is utilized 
throughout the time frame. 
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TABLE 41. ANNUAL CONTROL COSTS - COMBINED AREAS 



Mai URBANIZED ARE.A 
1 



KEF 

CTYI 

---I 



fleps. 



CUNTRQL rOST 
t C$/aCh£Ii 
I ?t>X I «iOX I JbX 



1 I AJAy 
21 AtlROKA 
3IBARRIE 

BRANTFQRD 
BURLt^''GTON 
CHATHAM 
WiNG'JACUUbY 

nanuRG 

DUNOAS 

ETDnlCQKE 

GALT 

GEDRGETCIt^N 
15IGUELHH 
IblHAMTLTQN 
ITIKINGSTOM 

iBUTTCri,-WATERLnO 

t^lUEAMlNGTDN 

gOlLlNnSAY 

ai unNnON 

' If 

36INIAGRA FALLS 
PTINDRTH 3AY 
PfllQAKVlLLE 

29IDfiILLlA 
30I0SHAWA 

31 IDWEW SOUNO 

32IPFTFK3QP0iiGH 

|3IP!CK£RING 

32 SSS"^ colbourne 

35IP0RT ERIE 
|36|PRE9T0N 

i3ri«rcMM0ND hill 

3fllST. CATHARINES I 

39I3T. THOMAS | 

aolSARlNA 

4liaLT. 5TF. MARltl 

asiarAflBORounH i 

a3tSTMcOE I 

aAlSTRATFORD I 

aSlSMDPuRY 

i!6|THUMDER 

«7 ITHRHNTU 

aeiTRENTHN 

491 WALLACEBuRG 




b4Y 



NFLl AND 
jilwt^ITttY 
52tKT^nSnR 
WOUOSTQCK 



53» 

n 



YORK 

YnHK,EAST 
YnRK,NORTM 



7130,34 
71 0,0 
71 0.0 
17122. ft£ 
0.0 
0.0 
0.0 



7 

36 
7 



3fil 13, B6 



7 

'? 

7 
7 
7 
7 
7 

17 
7 

36 



0.0 

o.n 

0.0 
0,0 
0,0 

olo 

0,0 
lfe.97 
23,^2 

0,0 

25, a3 



8:§ 



i 



,Oa«ii 
lo I 
.0 I 
.0a3bl 

•2 ' 

.0 I 
.0396 

•2 
'I 

: 

.0 

-" 

.0394 

• 043ci 

tO 

,0504 

.0 

,0b04 

.0i99l 

.03941 

_ ;8 I 

3a| 12.001 .03901 

41 126.751 .0521 I 

71 0.0 1,0 I 

7 I 0,0 I ,0 I 

710.0 1,0 I 

7128, fl7l ,04201 

17 I 0,0 I ,0 I 

71 0.0 I ,0 I 

361 0,0 1.0 I 

36129. 7JI ,0bO4l 

71 0,0 I ,0 I 

71 0.0 I .0 I 

36112. 8bl ,03^41 

3612^.831 ,04321 

38127,401 ,05031 

411 oin \\o 

7I2«.3/ ' 

36129,09 

71 0.0 
an 0.0 
41120, og 

7 I 30,161 ;0463l 
171 0,0 I ,0 I 
36124,051 ,049i 
38127.56 



,04331 
.05041 

•2 ' 
,0 I 

.04761 



0,0 1,0 



,04821 



I I 



7 
3o 

361 

I2fl,2ei ;04ft3 
127,791 04821 
7 I 0.0 1 .0 I 



24,16) ;P4H6 



10 



k 



66. 
0. 

0. 

0. j 

45, t 
69.1 

90,1 



¥\ 
8:1 

3?. I 

98, I 

2*' 
0,1 

n.i 

62, I 

0.1 

0, I 

0,1 

1U5, I 

0, 
34,1 
79.1 
«**.! 

0, I 

103,1 
0,1 
0.1 

62.1 
9 



61 ,1 

0,1 

94,1 

"•I 



339 



?01, I 
0.1 

102, 

I- 
I' 

0, 

122; 

204. 

0, 
316, 

0, 
371 
136. 

144. 

2- 

85. 
362. 

0. 



11132. 

0, 
597, 

275. 

2- 
2* 
2* 
2- 
2* 

3d6, 
599, I 
0, 
1116. 

11307; 

I 374, 
'85. 

I 0, 
I 224, 
I 1334, 







0.1 
236.1 

0.1 
0, I 
0.1 



370.1 1304, 

0. I 0. 

0.1 0, 

92.1 246, 

233.1 685, 

339.IU93. 

0.1 0, 

317.11056, 

362, f 1276, 

2* 

0, 

713, 

1125. 

0, 

966, 

1021. 

0, 

0, 
J055. 
1035; 

Of 





o; 

2p; 

0, 
282. 

306. 

o; 

"^ 

0. 

316. 

3io; 

0. 



0. 
0,1 

o;t 

673, I 

2-' 

Pi 



WEIGHTED AVE, I 123.221.04681 75,1 240. 



775. 
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TABLE 42. ANNUAL CONTROL COSTS - STORM AREAS 



\ 



TY 



rfiEFS. 



ra.TPQ^^r^ 



25% 



CS/Ai 
1 50X 



ST 1 
7bX I 



1 ! AJAV 

3ldARniF 

aioELLEVlLLF 
5iaRAMHTQN 

blrtRA^JTFQPQ 

7lriijRLlWGTQW 

SICMATHAM 

9iCHlrJGHArQn3y 

11 tOMNOAS 
laiETOfllCOKE 

HlSEoRGFTOhN 

tSlGUELHH 
IfalriAMTUTUN 

17 IKINGSTQN 
UIKITrh.-WATbRunQ 
19 LPAMJMGTON 



LINDSAY 

nNno^" 




ablNIARHA FALLS 
P/INORTH BAY 
28liJ«KVll.LF 
qiQRII.LlA 
OKINAWA 

'WEN ^Q"ND^ 
FTFRBaPOuGH 

iPTCKEfilNG ^ 
IPORT CULbnURNF 
[PORT HRTE 
I PRESTON 
IRTC^MClNn HILL 
131. CATHARIWES 
lai, THUJAS 

ISaPtna 
Hi ISI,T. STF. MARIE 
aaiSCAPbDHnuGH 
a3 STMrOF. 
«a STRATFORD 



30 
3i 
3S 
33 

3a 
35 

^^ 

3fl 
^9 

ao 



asisuoBuMY 



BAY 



abl TMUNOFI 

ja? iTnnnNTQ 

labiTREWT'^i^ 

1^91 HALLArEHURG 
1501 nFLLAND 
|bl 1 WHITdV 
l^aihlNntjOK 

I^JlHOOnSTOrK 
1 5a ! YORK 
isSl YHRK^rAST 
i5ol yHHK, NORTH 

I *li*FIRHT£n AVe. 



71 lO.^ai .03801 
7 1 9 . *i 9 I , i S I 

7 1 ib.sai .n493i 

I? I 5.ai ! ,03691 

71 8, an .03771 
_ ^^ 



3tt 
7 



7 
17 



0379 
\l I .03SSI 
, 03^591 
, "3*^31 

D3fl9" 



7 

I u 

I a, 591 

111.341 

171 ia.36l , 

7 I 19.001 ,039o! 

7111.141 .03^31 

7! l*i.99| .0397 1 

I I5i.a«| .03901 

I a.pii .03*)£>i 

I 5,971,03^91 
7! I2.^<el .03*^31 
3fll 7.67l,03flll 
171 12.901 ,n3«7l 
iBI 7,aft|.0375| 

71 U.?4I .03831 

7ltl.24i;n5»iSi 
38! 3. flat .03bai 

an 9.?oi.oai8i 

7 I 15.001 .03flbl 
71 19, 691 .039b| 
9;?9I. 03771 
7.101 .0371 I 
15.901 .039St 
15.721 . 03*^31 
14,191 ,03951 
«.S3l,'^3«l- 



7 1 

71 
171 

71 
38 1 
38! 

71 

71 
TBI 
18 1 
361 
411 

71 
381 

7 I 

ai I 

an 

71 



17.851 ,n39ai 

12.681 .OjHfal 
4.621 .03681 
4.31 I .03641 
5.831 .03721 

10.891 .04^11 
6,^1 I .03TI t 
9,?0I .03«2I 

22.491 .03981 
9.631 .04221 
4.801 .0407 1 

I3.3bl .Ojafti 

l7.04i;0597l 
6.051 ,03731 
5.9bl ,03721 

17,421 .03^)3 



I 16. lb 
1 7.?l 
1 6.14 

in -fle 



.03^91 

.03731 

.03701 
.03851 



29' 
27' 

141 
34i 



1 



I 9.951.03HUI 26. 



bfl 



180. 
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TABLE 43. ANNUAL CONTROL COSTS - UNSEWERED AREAS 



I 

NU I U^'bftNlZED AREA 



KEP 
CTY 



no 



C763 I rUNT^UL COST 
EFS. 1 fi/ACHF) 

?5% I «SOX I 75% 



Oibal 
035il 

03551 
03fell 
OShi J 

03f)l I 
03591 
0351 I 
0?(53I 
03591 
03531 
03531 
035ai 
I 
I 
3581 
035a| 
0351 I 
035b| 
O353I 
03531 
03581 
03541 

03571 
0391 I 

03561 

03531 

0361 I 

O35J1 

0351 I 

03531 

035il 

03521 

03531 

03531 

03^1 I 

035^1 

t 

03861 

03591 

03551 

0353! 

03fl7l 

039UI 

1 

0351 

03521 

03521 

03531 

03571 

03521 

I 

I 

n I 



1 1 
I 



a 

51 
61 
71 
SI 
91 

lot 
u 

la 
115 

1^ 

18 
19 

ao 

li, 

341 
251 
261 
271 
2AI 
291 
301 
111 

331 

3a I 

351 

1^ 

36 
39 

ao 

42 

43 
aa 
45 

2^ 

481 
491 
501 
51 
52 
531 
5m 
551 
56 1 
I--I 



AJAV 

AUHHRA 

BARfllE 

BPLLtVlLLF 

BPAMPTOm 

BRawTFQRD 

BliRLlNGTON 

CHATHAM 

CHIWGUACOUSV 

COf^nuRG 

DUNDAS 

ETQBICQKE 

GALT 

GEOfySETOwN 

GUELPH 

hartlton 

kinrston 

kttch.-uaterl 
lfamingtqm 

ClNnSAY 

DMOQ'^ 

AHKAM 
rilDLAND 

mts9i9sauga 
newmapket 
ntagra falls 

NORTH dAY 
OAKVILLF 
ORILLIA 
QSH&rs'A 

OWEW Sl!i!Nn 
PFTFHBQRDllGH 
PltKERlNG 
PORT COLPnuRN 
PORT EKIE 
PRESTON 

Richmond hill 

ST. CATHARINE 
ST. THUWAft 

saAtna 

slt, stf. mak 

scarborough 

SlMcDE 

STRATFORD 
3MDRURV 
THUWOER BAY 
TORONTO 
TRENTON 
WALLACEPIJRG 
WELl AND 
WHITBV 

WTNDSOK 

wnpnsTQCK 

YORK 

YORK, East 

YORK, NORTH 



00 



3 

TE 



2,69 
2.7b 
?.'S2 



17 


3.01 


7 


3.P5 


76 


3.50 


7 


3,74 


76 


3,26 


7 


3.58 


'? 


2.48 


2,fr0 


7 


3.49 


7 


2.69 


7 


a.fti 


7 


2,85 


1^ 


0.0 


0.0 


d 


3,66 


2.1H 




?.a2 


2.60 


7.17 


7 


3.75 


7 


2.f»7 


38 


2.<»5 


41 


3.45 


7 


3.10 


7 


2.^5 


7 


4.22 


f 


2. "52 


'? 


?.3l 


2.53 


38 


2. Si 


5B 


2.49 


78 


2.79 


2,&1 
3:55 


76 


2.57 


7ft 


0,0 


41 


3.03 


^l 


5,56 


2.57 


7 


3.12 


41 


2.89 


"1 


2,64 


0.0 


17 


2. 59 


^8 


2.19 


76 


2.49 


7 


2.79 


78 


2.88 


381 


2,49 


7! 


0,0 


7i 


0,0 


7 1 


O.n 



7 
7 
6 
7 

9 

9 
8 
8 

Q 

6 
9 
6 
6 
7 


9 
5 
6 
8 
k 
8 
8 

8 
7 
9 

8 

7 
10 
7 
6 
6 
8 
6 
7 
h 
9 

6 


fa 

Q 
6 
8 
7 
7 

fe 
5 
6 
7 
7 
6 

n 





16, 
16, 
15, 

1*; 

23, 

21. 

20. 
20, 

22. 

14. 
15, 
21, 

II: 

'I' 
ai: 

13. 

\h 
IS: 

24. 
IS. 
17. 
26. 

^Z- 

15. 

15, 
16, 

il: 

21. 
21 . 

'|: 
1?: 

0, 
14, 
15, 
14, 

16. 
16. 

0. 

0, 



3B. 
40, 
i7. 
47, 
58. 
5?. 

49. 
53, 

i^: 

52, 
S», 

37, 

40. 

0, 
54. 
31, 
34, 

46. 

45, 
49. I 
38,1 
43.1 
65. I 
45.1 
42.1 
63.1 

41. i 

32.1 
36. I 
35.1 
35.1 
39,1 
37.1 
53.1 
36.1 

2' 
55, 

53. 

36, 

aQ. 

0. 

33, 

31. 
35. 
39, 
39; 
35. 
0. 

0. 



li»^Fir,HTEO AVE. 



2.961 .03601 



7. 



1#». 



44. 



Because of this definition, the number of overflows may increase 
with an increase in treatment rate as shown in Figure 25. If the treatment 
rate Is high enough to deplete storage after the first overflow, then the 
event is over. When storage is utilized later a new event starts, and any 
overflow occurring in this event is considered separate from the first 
overflow. Thus, the number of overflow events was increased from one to 
two events, even though the treatment rate was increased. 

The number of overflow events appears to provide a irore meaning- 
ful parameter if the event is defined differently than the definition used 
by STORM. The overflows shown in case 2 in Figure 25, should be considered 
as a single event since they occur so closely together. Based on the U.S. 
assessment, a storm event terminates after 12 hours of no precipitation 
[2], Using this definition, the approximate relationship between overflow 
events and percent volume control Is derived as shown in Figure 26, for 
Burlington and Kingston. Thus, a rough approximation of the relationship 
between percent R and overflow events (OE) is d(OE)/dR = -1. 




TREATMENT- 0.03 In/hr {.076 cm/hO 
-•—TREATMENT » 0.01 In/hr (.025 cm/hr) 



CASE l'(T* 0.01 in/hr)[ 025cm/hr) 
NO OF OVERFLOW • ! 
VOL. OF OVERFLOW • 

.45in.[i.l43ciT,) 

■'*^ — STORAGE • OlOln. 
(.224 cm) 

^ OVERFLOW 



2S 30 

CASE 2 (T'O 05 ifi /ht)l 076cm/liO 
NO. OF OVERFLOW • 2 
VOL OF OVERFLOW '.IBin 

^STORAGE . O.IOIn * '^"''^"^ 

(.254 cm) 



hO'jrs 



,o 



-o r 



-o C 



FIGURE 25. EFFECT OF STORAGE AND TREATMENT CAPACITY 
ON NUMBER OF OVERFLOW EVENTS. 
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40 60 80 

% CONTROL OF VOLUME, R 



FIGURE 26. NUMBER OF OVERFLOW EVENTS vs PERCENT CONTROL. 
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5,5 The Overall Cost Assessment 

5,5,1 Overan results 

General results thus far are summarized in Tables 44, 45, 46, 
47 and 48, 

The only remaining problem is to estimate the province wide costs 
for 25, 50, and 75 percent control. As a first approximation, assume that 
an overall 25 percent control level is achieved by 25 percent control on 
the combined (A^), storm {A ) and unsewered (A^) areas. Thus, the 
approximate total annual costs, TAG, are: 

(TAG) 25 = 75A^ + 26A2 + 7f<^ 

= 75 (71,000) + 26 (210,000) + 7 (114,000) (77) 
= $ll,583,000/yr 

Likewise, 

(TAC)^Q = $33,372.000/yr (78) 

(TAC)^g = $97,841 ,000/yr (79) 

Recall that the cost of wet-weather control using secondary 
facilities is: 

1 - ke^^l (80) 

where: Z = annual cost, dollars per acre, 
k 8 = constants, and 
R = percent BOD removal (0 < R < 85). 

The cost of wet-weather control in tenns of pounds of pollutant removed, 
w, is: 

^s ' ^^^ *"^* ^ *^^* 

The marginal cost of BOD removal is: 

S=MLAJ^^554^ (82) 

However, the optimum mix of control of storm runoff from combined 
storm and unsewered areas is found by equating marginal costs. Using 
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TABLE 44. GENERAL INFORMATION 



Total Urbanized Area 
Total Population Area 
Average Population Density 
Average Precipitation 



586,000 acres (237,000 ha) 
4,720,000 persons 

8.07 persons/acre (19,9 persons/ha) 
32.75 inches/year (83.2 cm/hr) 



TABLE 45. LAND USE BY TYPE OF USE 



Use 



1,000 Acres 



(ha) 



Undeveloped 
Residential 
Commercial 
Industrial 

Other 
Total 



191 
» 

m 
Jl 

586 



( 77) 
( 85) 
( 16) 
( 22) 
( 37) 
(237) 



TABLE 46. LAND USE AND POPULATION BY TYPE OF SEWERAGE SYSTEM 



1,000 Acres (ha) 



1,000 Persons 



Developed 
Population 
Density 
Persons/ 
Acre (ha) 



Undeveloped 
Combined 
Storm 
Unsewered 



191 
71 

m 

TT4 



77) 
29) 
85) 
46) 



586 



237 





1,773 

2,485 

468 

4.726 







24.9 

11.9 

4.1 

12.0 



(61.5) 
(29.4) 
(10.1) 
(29.7) 



TABLE 47. QUANTITY AND QUALITY OF SEWAGE AND STORM WATER RUNOFF 



Sewerage System 



Flow: in/yr (cm/yr) 
Sewage Storm Runoff 



Quality: Annual Pounds of 

BCD/Acre {kg/ha) 
Sewage Storm Runoff 



Combined 
Storm 

Un sewered 



33.4 (84.8) 16.2 

15.9 (40.4) 12.8 

5.5 (14. U) 9,5 



41.1) 1.545 
32.5) 736 
24.1) 255 



16.1 (40.9) 12.5 (31.8) 



743 



1,733 
826 
286 
834 



137.0 (153.7) 
25.7 ( 28.8) 
21.3 ( 23.9) 



44.5 



49.9) 



TABLE 48, ANNUAL CONTROL COSTS PER UNIT OF DEVELOPED URBAN AREA: $/ACRE ($/ha) 



Sewerage System 



Level of Control, % of Total 
25 50 75 



Coefficients 
K B 



Combined 


75 


(185) 


240 


(593) 


775 


(1915) 


23.2 


0.047 


Storm 


26 


( 64) 


68 


(168) 


180 


{ 445) 


9.9 


0.039 


Un sewered 


7 


( 17) 


18 


( 44) 


44 


( 109) 


3.0 


0.036 
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equation (82) with the subscript (1) denoting combined, (2) denoting 
storm, and (3) denoting unsewered, yields: 

100 6 w- 100 6„ w- 

100 ^ K sr-^ 1006 „k„ /-^ 

L_L e M^ ^ ^_2 M2 



"V "2 

100 e^ w- 

100 k_ -^-^ 

P rr^-^ e "3 (83) 

"3 

If the above approximation is used, and marginal costs are 
compared for say, 50 percent BOD removal, one obtains: 

MP _ 100(0.047) (23.2) 100(0,047) (0.5) ioa) 

"^1 13770 ® ^ ' 

= $8,35/lb BOD ($18.39/kg BOD) 

„(. ^ 100(0.039) (9.9) 100(0.039) (0.5) /qs) 

2 25.7 

= $10.56/lb BOD ($23.26/kg BOD) 



MP _ 100(0.036) (3.0) 100(0.036) (0.5) 
"^3 2TT3 ^ 



:86: 



= $3.07/lb BOD ($6.76/kg BOD) 

This result indicates that, to achieve 50 percent control, 
unsewered and combined sewer areas should be controlled more intensively 
due to their relatively low marginal costs. Storm sewer areas should be 
the least intensively controlled due to their relatively high marginal 
cost. 

The correct solution can be found by solving for w, and 
Wj as functions of Wg. i.e., Sc 

W] = a^ 2 * ^ 2 ^2 *^^^ 

**2 " ^3 2 ^ '^3 2 *^2 <^^* 

Ml Bo k- M- 



where: 



h "1 

"r "2' "3' h' ^2' ^3' ''r ^2' h^ ^' ^2' "3 ^^^ ^^ defined 
earlier. 

The total wet-weather pollution load, WP, is 

WP = \ M/i *89) 

where: M. = annual pounds per acre from i area, and 

«"' * 4th 

A. = area of i area. 

Let p denote the proportion of WP that one wishes to control. Then, the 
optimum solution for a given p is found by substituting equations (87) and 
(88) into (89) or 

p(WP) = w^ A^ + Wg A^ + W3 Aj (901 

P(WP) Ma^ 2 ^ b^ 2 ''2*^ ^ ^2 ^2 ^ *^3 2 ^ ^3 2 ^2^'^3 *^^* 

^ * ^(WP) - a, 2 A, - ^3 2 ^3 (92) 

^2 - b^ 2 ^ ^ ^2 ^ ^3 2 ^^3 

Knowing w* the optimum pounds of pollutant removal for area 2, one 
can find the optimum levels of pollutant removed for areas 1 and 3, 
w* and W3, by substituting into equations (87) and (88), 

The results of the Ontario assessment indicate the values for the 
parameters shown earlier. Using these data, one obtains 

a) WP = M^ A^ + M^ A^ + M^ A3 

= 137.0(71.000) + 25.7(210,000) + 21.3(114,000) 

WP = 17,552,000 lbs BOD/yr (7.969,000 kg BOD/yr) 
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- 137.0 , p. 0.039 . , 9.9. J 37.0 ,, 
100(0.047) '"'-^0.047' ^23. 2^ ^"ZSTl^^ 



Similarly, 



a^ 2 = IS-^^ 



'3 2 = ^-^3 



c) b, „ and b^ „ 
62 M 



f O.039 . . 137.0 . 
^0.047' 45.7 ' 



Similarly, 
Thus, 



b^ 2 = 4.42 



b^ 2 = 0-898 



0r 



w* = 0(17,552,000) - (18.52)(71,000) - (6.43)(n4,000) /q. 

2 4.42171,000) + 210.000 + 0.898(114.000) 

w* = 28.03p - 3.27 (94! 



Then, substituting into equations (87) and (88) to find w* and w* 
yields: 



w* = 123. 89p + 4.07 (95) 

w* = 25.17P + 3.49 (96) 

Let (R*)p = 100 (w*)p/M^. denote the optimum percent control of 
the i^h source for control level, p. Then, 

(R*^ - 10Q(123.89p + 4,07) ^^ ^ ^ ^ ,„7. 

Vq 1"37~0 ^ 90. 4p + 3.0 (97) 



. 100(28.03p - 3.27) ._ ^^^^^^ . ^.^ (98) 

(R,)p = 1Q0(25.17p . 3.49) ^ ^^q^^p . 16.4 (99) 

3 '^'•■^ 

Let X . = A.M./WP for 1 = 1. 2, 3 denote the proportion of 
total pollutant load from the three areas. 

^ 137.0(71,000) ^ POO) ^ 

^1 17.552.000 "*^^^ 

, = 25.7(210,000) ^ (101) 

^2 17.552,000 ^' 

, , 21.3(114.000) , 0.138 (102) 

3 17.552.000 1.000 

The optimum percent control of 25. 50 and 75 percent is shown in Table 49, 

TABLE 49. OPTIMUM PERCENT CONTROL FOR SPECIFIED OVERALL PERCENT CONTROL 



Optimum Level of Control 3 

I 

Level of Control R^ ^ ^ i=l 



I Ri X. 



0,25 25.6 14.6 46.0 25.0 

0.50 48.2 41,9 75.5 50.0 

0.75 73.8 72.6 85.0 75,0 



Knowing (R.)p. one can find the cost per acre by simply substituting 
into equation (80), i.e., 

Z^ = ke^^ (103) 

to obtain the optimum annual cost per acre as shown in Table 50. Thus, 
the optimum annual control costs are shown in Table 51. 

5.5.2 Potential savings due to multipurpose planning 

The cost of wet-weather quality control can be reduced by 
integrating this purpose with dry-weather treatment plants and/or storage 
reservoirs for storm water quantity control. Dry-weather sewage treatment 
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TABLE 50. OPTIMUM ANNUAL COST PER ACRE FOR SPECIFIED PERCENT CONTROL 

Type of Optimal Annual Cost/Acre (ha) for Specified % Control 

Sewerage System 25 50 75 



Combined 


77(190) 


224(553) 


745(1,841) 


Storm 


17{ 42) 


51(126) 


168( 415) 


Unsewered 


16( 40) 


45(111) 


64( 158) 



TABLE 51. OPTIMUM ANNUAL CONTROL COSTS 
Level of Control Optimal Annual Cost: $/yr 



25 10.861,000 

|g 31,744,000 

1^ 95.471,000 



plants are designed to handle the peak flow anticipated 10 to 15 years 
after construction. The full capacity of these plants is seldom utilized 
because peak flows occur infrequently and also because additional capacity 
is frequently added before the actual flow approaches the design capacity 
of the plant. Provision of storage to equalize peak flows can greatly 
enhance the effective capacity of the existing treatment units. 
Utilization of this excess capacity can reduce the treatment capacity 
needed for wet-weather quality control. Similarly, utilization of storage 
available for wet-weather quantity control can result in reducing the 
storage and treatment requirements for wet-weather quality control. 

Preliminary results from the U.S. assessment, wherein excess 
capacity of 10 mgd in a 20 mgd plant was assumed and the city needed to 
store the excess (over natural) runoff for a two-year, 24-hour storm, 
indicated that it was possible to achieve significant savings due to 
multipurpose planning [2]. For these assumed conditions, the savings as 
compared to the single purpose venture are as shown below: 



Savings as % of 
Level of Control Singl e Purpose Control Costs 

25% 54% 

50% 50% 

75% 30% 

These results are suggestive only. Specific studies are needed to refine 
these rough estimates. 

5.5.3 Tertiary treatment versus wet-weather treatment 

The optimum mix of tertiary treatment and wet-weather control 
can be found by equating the marginal cost of tertiary treatment with the 
marginal cost of wet-weather pollution control. The estimated total 
annual incremental cost of a tertiary treatment plant is [2]: 

C^g^^ = 87,000 D°*^^^ 1104) 

where: C^ ^ = total annual incremental cost of tertiary treatment 
tert 

plant, dollars per year, and 
D = plant size, mgd. 

Assume a 25 mgd plant. Then, C^ ^ = $1,096,000 per year. The plant 
is assumed to increase the BOD removal from 85 percent to 95 percent or 
about 0.017 pounds (7.71 g) per capita-day or 1,550,000 pounds (704.000 
kg) per year for the city of 250,000 people. Thus, the unit cost of 
tertiary treatment, c. ., is $0.71 per pound ($1.56 per kg) of BOD 
removed. 

Equating the marginal cost of wet-weather control to the unit 
cost of tertiary treatment yields: 

c^ ^ = la^e'l^ (105) 



tert M 



ai* 



M c ^(M) 

"* 100b '" hooe (K)-" 

where: w* = optimum pounds of wet-weather pollution to control prior 
to using tertiary treatment. 
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. -n-aiu. . _..^ 



The optimum percent control in terms of R, is: 

1 c^ (M) 

The overall average BOD loading per acre, M. is 

- MP _ 17,552,000 

^ (A^ + A^ + A3) - 395.000 



Then, 



M = 44.5 lb BOD/acre (49.9 kg/ha) 

«1 = -X r^^ in io°o(-o!o4li(3.60)'°] = ''■' '^°«) 

Thus, for these assumed conditions, approximately 16 percent of the 
wet-weather pollution should be controlled prior to initiating tertiary 
treatment. While these results are for one specific set of assumptions, 
they do suggest that it is highly desirable to do this tradeoff analysis 
before committing a community to tertiary treatment. 

1*^ Summary 

The purpose of this assessment was to evaluate the cost of 
controlling varying levels of wet-weather pollution emanating from the 
4,720,000 people in 56 cities in Ontario. Reliable procedures for assess^ 
ing storm water pollution were not yet available. Thus, a considerable 
amount of developmental effort was expended in devising such procedures. 
Major results are presented, by item, in the next paragraphs. 

1) Land Use - Using a definition of urbanized areas which includes 
population densities as low as one person per acre leads to 
inclusion of the relatively large amount of land which is 
undeveloped (about one-third of the total land). Residential 
development utilizes the majority of the developed land. 

2) Type of Sewerage System - ^ery limited data exist on the population 
and area served by various types of storm drainage systems. 
Population served by combined, storm, and unsewered areas was 
derived by assuming that the highest density areas were served by 
combined sewers, the intermediate level by storm sewers, and the 



lowest density areas were unsewered. The transition points were 
identified using available data on areas served by the three 
systems. This method would tend to overestimate the population 
served by combined sewers. 

3) Quantity of Storm Water - An average of 28.6 inches of water per 
year leaves Ontario cities, of which 12.5 inches comprise storm 
water runoff. Annual storm water runoff volumes exceed sewage flows 
in low density urban areas. Urban runoff is significant relative to 
sewage flows in low density urban development. Thus, on a volumetric 
basis, urban runoff is significant relative to sewage flows. 

4) Quality of Storm Water - Storm water pollution loads approach waste- 
water effluent loads after secondary treatment has been installed. 
The exact quantity of storm water pollution remains unknown due to 
lack of sufficient data. Numerous assumptions were needed to 
develop a general pollution loading equation. There seems to be 
general agreement that combined sewer overflows are much more 
serious than storm water runoff. However, the results of this 
assessment indicate relatively low loadings of storm water pollu- 
tants (about one-half of the load coming from secondary treatment 
plant effluent). Only through carefully conducted sampling programs 
can these estimates be refined. 

5) Total Single-Purpose Control Costs - Relatively detailed studies in 
four Ontario cities provided the basis for evaluating storage-treat- 
ment alternatives for wet-weather control. One year of hourly data 
was simulated to predict the performance of various storage-treatment 
configurations. These results were put in analytical form to 
expedite extrapolation to all urbanized areas. These results are 
combined with data on the cost of storage and treatment to derive 
the optimum mix of storage and treatment to use to obtain a given 
level of control. The final result is shown in Figure 27. A 
striking feature of the curve is that it bends upward (convex) 
indicating increasing incremental costs (particularly at higher 
levels of control). The primary reason the curve has this shape is 
due to the disproportionately larger amounts of storage and 
treatment required to control the larger storms. 
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*Assumes management plan integrating dry-weather quality control, 
wet-weather quality control, and wet-weather quantity control. 

FIGURE 27. TOTAL ANNUAL AND INITIAL CAPITAL COSTS FOR VARIOUS 
LEVELS OF WET-WEATHER POLLUTION CONTROL IN ONTARIO 

6) Total Multiple-Purpose Control Costs - Significant savings can be 
realized if one integrates dry- and wet-weather treatment and 
storage for quality as well as quantity control. The lower curve in 
Figure 27 indicates the cost of storm water quality control in an 
integrated system. This result suggests that the potential savings 
are significant enough to warrant further study in evaluating storm 
water systems. 

7) Tertiary Treatment versus Storm Water Quality Management - A compari- 
son of the marginal costs of tertiary treatment of sewage for 
further BOD control with initiating control of wet-weather quality 



143 



indicates that one should initiate some level of wet-weather quality 
control prior to using tertiary treatment. Of course, a different 
result would occur if nutrient control is used instead of BOD 
control. Nevertheless, the relative low marginal costs of wet- 
weather control at low levels of control indicate that it should be 
given serious consideration as an alternative to tertiary treatment. 

5,7 List of Variables 

a coefficient (inches per hour) 

A:. combined sewer area 

^w storm sewer area 

A- unsewered area 

A . total developed area 

tot 

AR annual runoff (inches per year) 

fe coefficient (inches per hour) 

t coefficient in cost equation 

€ coefficient (percent R ) 

c unit cost of tertiary treatment (dollars per pound) 

tert 
c unit cost of storage (annual dollars per acre-inch) 

G_ unit cost of treatment (annual dollars per inch per hour) 

C. ^ Total annual incremental cost of tertiary treatment 

tert 

plant (dollars per year) 

(1 coefficient (inch ) 

D plant size (mgd) 

ENR Engineering News Record Cost Index 

ri treatment plant efficiency 

f coefficient (percent R)~ 

f(R,:S,T) Production function relating percent pollutant control 

(R ) to storage (S) and treatment (T) 

I coefficient 

K coefficient 

I breakeven pipe length 

X proportion of total pollutant load from i area 

§ annual pounds of pollutant (pounds per acre-year) 
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M average annual pollutant loading {pounds per acre-year) 

K marginal cost 

MRS^^ marginal rate of substitution of storage for treatment 

QE number of overflow events per year 

^. gross population density 

PD^ population density in developed portion of urban area 

q coefficient 

I percent runoff control 

R^ percent pollutant control 

^ optimum percent of pollutant control prior to using tertiary 

treatment 

t maximum percent pollutant control 

f proportion of WP which is controlled 

i coefficient 

S storage volume, inches 

S* optimum storage volume (inches) 

t treatment rate (inches per hour) 

y* optimum treatment rate (inches per hour) 

T^ treatment rate at which Isoquant is parallel to the ordinate 

(inches per hour) 

T^ treatment rate at which isoquant intersects the abscissa 

(Inches per hour) 

TAG total annual cost: $/year 

t annual pounds of pollutant removed 

#* optimum pounds of wet-weather pollutants to control prior to 

using tertiary treatment 

W total wet-weather pollutant load (lbs/year) 

^ coefficient 

z coefficient 

Z total annual cost (dollars per acre) 

I* optimum total annual cost (dollars per acre) 

Zg total annual cost for secondary control unit (dollars per 
acre 
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GLOSSARY 

Antecedent conditions - Initial conditions in catchment as determined from 
hydro! ogic events prior to stonm. 

Biological treatment processes - Means of treatment in which bacterial 
or biochemical action is intensified to stabilize, oxidize, and nitrify 
the unstable organic matter present. Trickling filters, activated sludge 
processes, and lagoons are examples. 

Catchment - Surface drainage area. 

Combined sewage - Sewage containing both domestic sewage and surface 
water or stonn water, with or without industrial wastes. Includes flow in 
heavily infiltrated sanitary sewer systems as well as combined sewer 
systems. 

Combined sewer - a sewer receiving both intercepted surface runoff and 
municipal sewage. 

Combined sewer overflow - Flow from a combined sewer in excess of the 
interceptor capacity that is discharged into a receiving water. 

Conservative - Non-interacting substance, undergoing no kinetic reac- 
tion; examples dre salinity, total dissolved solids, total nitrogen, total 
phosphorus. 

Convective precipitation - Precipitation caused by lifting due to 
convective currents, as in thunderstonns. 

Cyclonic precipitation - Precipitation caused by lifting associated with 
junctions of different air masses, as for instance, with most warm and 
cold fronts. 

Depression storage - Amount of precipitation which can fall on an area 
without causing runoff. 

Detention - The slowing, dampening, or attenuating of flows either 
entering the sewer system or within the sewer system by temporarily 
holding the water on a surface area, in a storage basin, or within the 
sewer itself. 
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Domestic sewage - Sewage derived principally from dwellings, business 
buildings, institutions, and the like. It may or may not contain 
groundwater. 

Economies of scale - Unit costs decrease as output increases. 

Equalization - The averaging {or method for averaging) of variations in 
flow and composition of a liquid. 

Expansion path - Locus of points connecting numerous isoquants indicat- 
ing the optimum combination of inputs. 

First flush - The condition, often occurring in storm sewer discharges 
and combined sewer overflows, in which a disproportionately high pollution 
load is carried in the first portion of the discharge or overflow. 

Frequency diagram - Curve which relates the number of occurrences of 
events to their magnitude. 

Initial abstraction - Initial precipitation loss including interception 
and depression storage. 

In-system - Within the physical confines of the sewer pipe network. 

Interception - Initial loss of precipitation due to vegetation. 

Isoquant Lines - Lines of equal cost. 

Isoquants - Curves representing combinations of the inputs yielding the 
same amount of output. 

Non-conservative - Substance undergoing kinetic interaction, assumed to 
be a first-order reaction; examples are biochemical oxygen demand (BOD), 
col i form bacteria, dissolved oxygen (DO). 

Orographic precipitation - Precipitation caused by lifting of an air 
mass over mountains. 

Orthophosphate - Phosphate that appears as PO^, HPO| or H^PO^, 
i.e. is hydrolizable. Creates a growth response in algae. 

Physical-chemical treatment process - Means of treatment in which the 
removal of pollutants is brought about primarily by chemical clarification 
in conjunction with physical processes. The process string generally, 
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includes preliminary treatment, chemical clarification, filtration, carbon 
adsorption, and disinfection. 

Pollutant - Any harmful or objectionable material in, or change in 
physical characteristics of water or sewage. 

Precipitation event - A precipitation event terminates if zero rainfall 
has been recorded for the previous specified time interval. 

Primary treatment - Process which removes about 35 percent of the 
biochemical oxygen demand of the waste. 

Retention - The prevention of runoff from entering the sewer system by 
storing on a surface area or in a storage basin. 

Runoff coefficient - Fraction of rainfall that appears as runoff after 
subtracting depression storage and interception. Typically accounts for 
infiltration into ground and evaporation. 

Sanitary sewer - A sewer that carries liquid and water-carried wastes 
from residences, commercial buildings, industrial plants, and institu- 
tions, together with relatively low quantities of ground, storm, and 
surface waters that are not admitted intentionally. 

Secondary treatment - Process which removes about 90 percent of the 
biochemical oxygen demand of the waste. 

Sewer - A pipe or conduit generally closed, but normally not flowing 
full, for carrying sewage or other waste liquids. 

Sewerage - System of piping, with appurtenances, for collecting and 
conveying wastewaters from source to discharge. 

Storm flow - Overland flow, sewer flow, or receiving stream flow caused 
totally or partially by surface runoff or snowmelt. 

Storm sewer - A sewer that carries intercepted surface runoff, street 
wash and other wash waters, or drainage, but excludes domestic sewage and 
industrial wastes. 

Storm sewer discharge - Flow from a storm sewer that is discharged into 
a receiving water. 
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storm water - Water resulting from precipitation which either percolates 
into the soil, runs off freely from the surface, or is captured by storm 
sewer, combined sewer, and to a limited degree sanitary sewer facilities. 

Surface runoff - Precipitation that falls onto the surfaces of roofs, 
streets, ground, etc., and is not absorbed or retained by that surface, 
thereby collecting and running off. 

Tertiary treatment - Process which removes about 96 percent of the 
biochemical oxygen demand of the waste. 

Urbanized area - Central city, or cities, and surrounding closely 
settled territory. Central city (cities) have population of 50,000 or 
more. Peripheral areas with population density of one person per acre or 
more are included. (United States city definition) 

Urban runoff - Surface runoff from an urban drainage area that reaches a 
stream or other body of water or a sewer. 

Wastewater - The spent water of a community. 
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APPENDIX II 
MAPS OF AREAS OF DATA TABULATION 
FOR SEVEN TEST CITIES IN ONTARIO 



Appendix II 
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FIGURE 28. AREA OF DATA TABULATIOM IN BURLINGTON 

(inside hatched line) 
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FIGURE 29. AREA OF DATA TABULATION IN GUELPH 

(inside hatched line) 
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KINGSTON 













FIGURE 30. AREA OF DATA TABULATION IN KINGSTON 

(inside hatched line) 



Tea 



KITCHENER 







v;l;:s 



.5 .5 1 





^^^^-^ r 



KiLO%'£rp.rs. d 



.5 



FIGURE 31. AREA OF DATA TABULATION IN KITCHENER 

(inside hatched line) 
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MILTON 
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FIGURE 32. AREA OF DATA TABULATION IN MILTON 

(inside hatched line) 
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FIGURE 33. AREA OF DATA TABULATION IN ST. CATHARINES 

{inside hatched line) 
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FIGURE 34. AREA OF DATA TABULATION IN SAULT STE. MARIE 

(inside hatched line) 
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APPENDIX III 
FIELD INTERVIEW OUTLINE 

A. General Information 

1. Community Information 

a. Demographic - population, land uses, etc. 

b. Community service area 

c. Major economic activities 

2. Physical Infonnation 

a. Topographic characteristics 

b. Prevalent soils - permeability, etc. 

c. General patterns of drainage - slopes, location in 
catchment, etc. 

3. Annual Climatic Information 

a. Precipitation- rainfall, snowfall 

b. Temperature distribution, etc, 

4. Local Governmental. Infonnation 

a. Organization of responsibility for: 

- design and construction of 

- storm runoff works 

- sanitary wastewater works 

- operation and maintenance of 

- storm runoff works 

- sanitary wastewater works 

5. Financing for storm and sanitary wastewater works 

a. Local financing methods 

b. Provincial grants 

- administrative vehicle 

- determination of provincial share 

- detenmination of local share 

- expense attributable to province 

- expense attributable to local jurisdiction 

c. National grants 

- administrative vehicles or agencies involved 
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- determination of national, provincial, and local shares 

- expenses attributable to national, provincial, and local 
juristictions 

d. Past grants experience 

B. Col 1 ection Systems Information 

1. Service area coverage 

a. Combined sewer 

b. Sanitary sewers 

£> Separate storm sewers 
d. Unsewered 

2. Mileage of Sewer Systems 

a. Combined sewerage 

b. Sanitary sewerage 

c. Separate storm sewerage 

3. Problems Experienced 

a. Hydraulic overloading 

b. System bypasses (and their purpose) 

- number 

- location 

- hydraulic capacity limitation 

- receiving water 

C« System solids deposition problems 

- general location 

- cleaning frequencies 

d. Infiltration and/or inflow 

- roof leader connections 

- other direct inflow locations - manholes, etc. 

- infiltration problems 

- estimated infiltration rates 

- general infiltration problem locations 

e. Storm runoff inundation areas 

- locations 

- cause {no collection systems, hydraulic overloading, etc.) 

- threshold runoff rate that results in inundation problems 
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4. Local Construction Practice 

a. Use of catch basins or inlets 

b. Construction of system relief bypasses 

c. Control of bypasses or overflows 

- regulators 

- detention facilities 

d. Use of detention facilities 

- Type - on-site. In-system, off-system 

- preferable locations for each type 

- approximate design criteria 

e. Use of retention facilities 

f. Practices concerning extension of combined sewer systems 

g. Practices concerning extension of separate sanitary sewer 

systems 
h. Practices concerning extension of separate storm sewer systems 

C. Dry-Weather Treatment Works 

a. Mean daily flows 

b. Treatment unit processes 

c. Equalization of influent flows 

d. Influent strength 

- BOD 

- COD 

- SS 

- Nutrients 

- Col i form 

e. Treatment removal efficiencies 

f. Effluent strength 

g. Hydraulic treatment plant capacity 
h. Bypass capability 

- how often? 

- bypassing due to? 

D. Receiving Waters For: 

1 . Runoff discharges 

2. Overland flow 

3. Combined sewer overflow 



4. System and plant bypasses 

5. Treatment plant effluents 

E. Wet Weather 

1. Control Activities 

a. Regulator operation and maintenance 

2. Sampling Activities 

a. Combined sewer overflows 

b. Urban runoff 

c. Bypasses 

3. Planning Activities for: 

a. Control of quantity 

b. Control of quality 

c. Abatement of runoff effects 

- storage/detention 

- types 

- likely locations for each 

- treatment 

- treatment types 

standard - sedimentation, disinfection 
exotic - air flotation, microstrainer, swirl 
technologies 

F. Interest In Runoff Planning 

1. Current local activities for investigating runoff-induced problems 

a. Quantity management (flood control, etc.) 

b. Quality management 

- effluents 

- bypasses 

- combined sewer overflows 

- runoff discharges 

- overland runoff flow 

2. Does interest exist in locally applied runoff planning tools? 

a. Would the local jurisdiction be able to produce inputs? 

b. Do they have computer capacity? 

I c. Would they use Ministry level assistance? 

t 
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G. Verification of Existing Reduced Data 

H. Miscellaneous 

1. Land acquisition costs 

a. Central city land costs for 

- residential 

- commercial 

- Industrial 

- open space 

b. Non-central city land costs for 

- residential 

- commercial 

- industrial 

- open space 

- undeveloped 
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APPENDIX IV 

PROVINCIAL GRANTS 

All municipalities in the Province of Ontario are eligible under 
provincial legislation for grants. 

I. General grants - These grants are available for general 
purposes and are not tied to specific uses. Most municipalities use them 
for relief from taxation by reduction of the mill rate. None of the 
municipalities interviewed applied for grants for construction and/or 
maintenance of urban runoff facilities. 

II. Subsidy under the Public Transportation and Highway 
Improvement Act for urban municipalities - The municipalities are eligible 
for a subsidy allotted annually by the Minister of Transportation and 
Communications for eligible expenditures made on all roads and streets 
subject to the limitations of MTC policy on subdivision roads and streets; 
which limitations are set out in DTC Circular 72-010. 

Part I - Construction 

A4 Drainage is covered by the subsidy and includes: 

(1) Open ditching, including off-take ditches to nearest outlet 

(2) Underdrains 

(3) Storm sewers, including pumping stations where required, subject 
to the limitations of MTC policy currently outlined in DTC 
Circular 71-040 

(4) Catch basins and connections 

(5) Curbs 

(6) Gutters 

(7) Municipal drainage assessment on roads 

(8) Stream improvement, if required beyond those limits defined in 
Section B 

B. Bridges, culverts, and grade separations subsidizable at bridge and 
culvert rate. 

3. Outlet sewers for underpasses, including pumping stations when 
required, subject to limitations outlined in DTC Circular 71-040 
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5. Stream improvement, if required, for 100 feet along bed of stream 
measured from the outer extremities of a new structure 

6. Stream diversion in lieu of structures only if approved by Head 
Office 

Part II - Maintenance 

A Roads and Streets 
3. Roadside 

(1) Entrance culverts - cleaning, repairs 

(2) Ditches, including off-take-cl eaning, repairs, relocation 

(3) Erosion control 

(4) Catch basins - cleaning, repairs, replacement 

(5) Storm sewers, subject to limitations of MTC policy currently 
outlined in DTC Circular 71-040 

(6) Underdrain cleaning, repairs, replacement 
(8) Curb and gutter repairs and replacement 
(12) Municipal drainage assessment on roads 

B4 Replacement of pipe culvert 

Part III - Overhead 

Overhead can be compensated for at 7 percent of the cost of items 
eligible for subsidy under Parts I and II. Details are set forth in Part 
III of published guidelines. 

In DHO OB-M-69 a schedule is set out giving allowable rates for 
equipment rental . 

In Circular No. 71-038 a list of approved protective clothing and 
safety equipment is given. 

Circular 71-040 sets out the subsidy on storm sewer construction 
and maintenance. 

Circular 71-020 outlines the conditions under which the Minister 
will authorize the payment of subsidy on municipal expenditures made on 
roads or streets opened or constructed by private interest, or opened or 
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constructed by municipality acting as a subdivider, or opened or constructed 
by the municipality acting as a constructor for the subdivider. Three cases 
are set out and eligibility can be determined by a decision-tree concept 
graph included therein. 

In addition there are certain municipalities, as of April 1, 
1975, that are eligible for provincial grants for water and sewage 
facilities through the Ministry of the Environment: 

- Municipality of Metropolitan Toronto 

- The Regional Municipalities of 
Durham 
Haldimand-Norfolk 

Hal ton 

Hami 1 ton-Wentworth 

Niagara 

Ottawa-Carl eton 

Peel 

Sudbury 

Waterloo 

York 

- The District Municipality of Muskoka 

- The cities of Thunder Bay and Timmins 

Sewage facilities include all treatment works operated by the 
municipality or on behalf of an area municipality or any local branch 
thereof, and intercepting and trunk sanitary sewers including ancillary 
structures which make practicable the matter of adequate but efficient 
pollution control (the trunk sewer must terminate at a sewage treatment 
plant, may serve as a local collector, and must have a theoretical 
capacity in excess of 6 cfs). 

Eligible costs include: 

(a) Cost of design and supervision 

(b) Actual construction costs incurred by the municipality 

(c) Interim financing charges during actual construction 

(d) Expenditures incurred in associated land purchases 

(e) Other approved costs (savings clause). 
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FEDERAL GRANTS 

Federal grants of interest to municipalities are (1) those 
administered by the Central Mortgage and Housing Corporation under the 
terms of the National Housing Act, 1973, and (2) grants administered in 
the past by different federal departments under the category of Winter 
Works or Winter Capital Projects Funds for New Employment Opportunities in 
the Construction Industry. This was to be succeeded by a new Local 
Initiative Program to reduce seasonally high unemployment in 1976 and 1977, 
with emphasis being given to municipal works and projects. (Announcement 
by Minister of Finance in budget speech of June 21, 1975.) 

(1) The sections of the National Housing Act noted below are of most 
interest to municipal and regional public works departments: 

Part III - Urban Renewal, Sections 22 through 27 

Part III.l - Neighborhood Improvement Program, Sections 27.1 
through 27.6 

Part VI, 1 - New Communities Federal-Provincial Agreements, 

Sections 45.1 through 45.3 

Part VIII - Loans for Municipal Sewage Treatment Projects, 
Sections 50 through 54 

(2) Federal-Provincial Winter Capital Project Funds 

This has been a federally financed program to provide essentially 
new employment opportunities for people largely in the construction industry, 

Capital works programs are those which would not ordinarily have 
been undertaken. The labor cost forgiveness is 100 percent for employment 
created between December 1 and May 31 in any year, and 50 percent in the 
summer season. The program is said to generate up to four times the 
amount of work (in dollars) of the allotment granted. 

The above program has been in effect from time to time since the 
1950's. On June 21, 1975, the Federal Minister of Finance announced in 
his budget speech, that a new Local Initiative Program with available 
funds up to $285 million was to be set up. 

The program would provide for projects over two eight-month 
segments, November /June 1976 and November /June 1977, to generate 400,000 
man-months of employment. 
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The program would be administered by the Department of Manpower 
and Immigration, and was being created to reduce high seasonal unemploy- 
ment in specific areas. It is a labour-intensive community betterment 
program. 

Federal funds will be available for approved projects sponsored 
by private groups or by municipalities with the emphasis on municipal 
works and projects. 
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